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1.0  SUMMARY 


The  technology  of  high  optical  quality  gain  elements  based  on  Co,  Ni,  and  Fe  doped  ZnSe  and 
ZnS  crystals  promising  for  room-temperature  lasing  was  developed.  Multiple  crystals  have  been 
fabricated  and  studied  via  absorption  and  photoluminescence  spectroscopy.  For  the  first  time 
lasing  was  achieved  on  Fe  ions  in  ZnS  and  ZnSe  crystals  co-doped  with  Co  via  effective  Co-Fe 
energy  transfer.  For  the  first  time  Co:ZnS  crystals  were  successfully  used  for  passive  Q- 
switching  of  continuous  wave  (CW)  pumped  Er:YAG  1645  nm  lasers. 

The  matrices  for  calculating  the  energy  and  intensity  of  electric  dipole  transitions  of  3 cf 
electronic  configurations  were  created.  These  matrices  enabled  calculation  of  the  energies, 
oscillator  strengths,  and  lifetimes  for  any  transition  of  the  3 dn  (n  =  2,  3,  4,  6,  7,  8)  electronic 
configuration.  This  thereby  enabled  calculation  of  the  absorption,  luminescence,  luminescence 
excitation  spectra,  as  well  as  the  absorption  spectra  for  transitions  from  the  excited  states.  Within 
the  energy  and  intensity  matrices,  the  configuration  curves  of  the  energy  levels  are  developed  by 
approximating  the  zero-  and  fourth-order  crystal  field  parameters  with  parabolic  dependences  on 
the  vibration  amplitudes.  The  configuration  diagrams  of  Co2+,  Fe2+,  and  Ni2+  ions  are  calculated 
for  different  cases  of  anhannonicity.  The  proposed  theory  is  applied  to  analyze  the  experimental 
spectra  of  Co2+,  Fe2+,  and  Ni2+  ions  in  ZnSe  available  in  the  literature  and  yields  good  agreement 
between  the  experimental  and  modeled  spectra. 
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2.0  INTRODUCTION 


Middle-infrared  (mid-IR)  laser  sources  are  in  great  demand  for  a  variety  of  applications 
including  molecular  spectroscopy,  non-invasive  medical  diagnostics,  laser  medical  scalpels  that 
make  use  of  the  ability  of  the  laser  to  tune  in  and  out  of  the  strong  absorption  band  of  liquid 
water,  industrial  process  control,  environmental  monitoring,  atmospheric  sensing  and  free  space 
communication,  oil  prospecting,  and  numerous  defense  related  applications  such  as  infrared 
aircraft  countermeasures,  monitoring  of  munitions  disposal,  stand-off  detection  of  explosion 
hazards,  eye  safe  seekers  for  smart  munitions  and  cruise  missiles,  covert  communications 
systems,  and  injection  seeding  of  high  power  chemical  (e.g.  HF)  lasers. 

Mid-IR  wavelengths  are  typically  generated  using  relatively  complex  nonlinear  optical 
conversion  techniques  or  by  means  of  direct  generation  in  hetero -junction  lead-salt,  antimonide, 
or  quantum  cascade  semiconductor  lasers  featuring  limited  output  power  and  tuning  range. 
Impurity  doped  crystalline  vibronic  lasers,  having  a  gain  bandwidth  up  to  50%  of  central 
wavelength,  constitute  another  viable  route  for  broadly  tunable  mid-IR  coherent  sources.  Among 
them  are  transition  metal  (TM2+,  e.g.,  Cr2+  or  Fe2+)  doped  binary  (e.g.,  ZnSe,  ZnS,  CdSe,  CdS, 
ZnTe)  and  ternary  (e.g.,  CdMnTe,  CdZnTe,  ZnSSe)  chalcogenide  crystals  that  represent  a 
relatively  new  class  of  solid  state  gain  media  with  strong,  ultra-broad  absorption  and  emission 
bands  in  the  mid-IR  range  of  the  optical  spectrum.  In  1996,  Scientists  from  the  Lawrence 
Livermore  National  Laboratory  [1]  were  the  first  to  show  that  among  the  different  types  of 
crystalline  gain  materials,  TM”  doped  wide  bandgap  II- VI  semiconductor  crystals  could  be  very 
special  for  mid-IR  lasing: 

•  An  important  feature  of  II- VI  compounds  is  their  tendency  to  crystallize  as  tetrahedral 
coordinated  structures,  as  opposed  to  the  typical  octahedral  coordination  at  the  dopant  site. 
Tetrahedral  coordination  results  in  a  smaller  crystal  field  splitting  which  has  the  effect  of 
placing  the  dopant  transitions  further  into  the  infrared  (IR). 

•  Another  key  feature  of  these  materials  is  that  the  heavy  anions  in  the  crystal  enable  a  very  low 
energy  optical  phonon  cutoff  which  makes  them  transparent  in  a  wide  spectral  region  while 
decreasing  the  efficiency  of  non-radioactive  decay.  This  gives  the  promise  of  a  high 
fluorescence  yield  at  room  temperature  (RT). 

Active  interest  in  transition  metal  (TM)  doped  II- VI  compounds  inspired  by  [1]  was  explained  by 
the  fact  that  these  media  are  close  mid-IR  analogs  of  titanium-doped  sapphire  (Ti-S)  in  terms  of 
spectroscopic  and  laser  characteristics.  It  was  anticipated  that  TM“  doped  chalcogenides  would 
lase  in  the  mid-IR  with  a  great  variety  of  possible  modes  of  oscillation,  similar  to  what  has  been 
found  with  the  Ti-S  laser.  Indeed,  optically  pumped  RT  lasers  based  on  Cr:ZnS,  Cr:ZnSe, 
Cr:Cdi_xMnxTe,  Cr:CdSe,  and  Fe:ZnSe  crystals  providing  access  to  the  2-3  and  3.8-5  pm  spectral 
regions  with  high  (up  to  70%)  efficiency,  multi- Watt-level  output  powers,  tunability  in  excess  of 
1000  mn,  and  narrow  spectral  linewidth  (<20  MHz),  have  been  reported  by  several  groups  (see 
e.g.  reviews  [2,  3,  4,  5,  6]  and  references  therein).  Arguably  these  lasers  represent  nowadays  the 
simplest  and  the  most  cost-effective  route  for  high  power,  broadly  tunable  lasing  in  the  2-3  and 
3.8-5  pm  wavelength  ranges. 

There  are  still  several  challenges  to  be  overcome  in  the  development  of  affordable  high  power 
tunable  TMTI-VI  lasers  continuously  covering  the  2-5  pm  spectral  range.  Among  them  are  the: 
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1)  Lack  of  developed  gain  materials  capable  of  effectively  operating  over  the  currently 
unavailable  3-4  pm  spectral  range. 

2)  Lack  of  convenient  pump  sources  for  3.8-5  pm  Fe:ZnSe  lasers 

3)  Lack  of  TM:II-VI  gain  materials  pumped  by  the  radiation  of  readily  available 
neodymium  or  ytterbium  lasers. 

To  address  the  enumerated  challenges,  the  following  research  objectives  are  proposed: 

1 .  Development  of  the  technology  associated  with  high  optical  quality  gain  elements  based  on 
Co2+  doped  ZnSe  crystals  that  are  promising  for  room-temperature  lasing  over  the  3-4  pm 
spectral  range. 

2.  Fabrication  and  characterization  of  Co"  :ZnSe  experimental  samples  that  lase  over  the  3-4  pm 
spectral  range  under  Er- liber,  808  mn  diode,  and/or  Alexandrite  laser  pumping. 

3.  Development  of  the  technology  associated  with  high  optical  quality  gain  elements  based  on 
Ni"  co-doped  Co"  :ZnSe  crystals  which  are  promising  for  room-temperature  lasing  over  the 
3-4  pm  spectral  range  via  energy  transfer  from  Ni  to  Co. 

4.  Fabrication  and  characterization  of  Ni"  co-doped  Co"  :ZnSe  experimental  samples  for  3-4 
pm  lasing  under  Yb-fiber  and/or  Nd  laser  pumping. 

5.  Development  of  the  technology  associated  with  high  optical  quality  gain  elements  based  on 
Co2+  co-doped  Fe:ZnSe  crystals  which  are  promising  for  room-temperature  lasing  over  the 
3. 8-5. 2  pm  spectral  range  via  energy  transfer  from  Co  to  Fe. 

2d”  2d” 

6.  Fabrication  and  characterization  of  Co  co-doped  Fe"  :ZnSe  experimental  samples  for  3.8- 
5.2  pm  lasing  under  Er-fiber  and/or  808  mn  diode  laser  pumping  via  energy  transfer  from  Co 
to  Fe. 

The  specific  tasks  to  be  performed  at  the  University  of  Alabama  Birmingham  (UAB)  in 
collaboration  with  the  A.M.  Prokhorov  General  Physics  Institute  of  the  Russian  Academy  of 
Science  business  representative  towards  realization  of  the  proposed  objectives  are  as  follows: 

Task  1.  Analyze  and  evaluate  the  basic  parameters  relevant  for  fabrication  of  TM  doped  ZnSe 
crystals  by  post-growth  thennal  diffusion  of  Ni,  Co,  and  Fe  ions  in  high  optical  quality  IR 
window-grade  polycrystalline  ZnSe.  (UAB) 

Task  2.  Theoretically  analyze  energy  level  diagrams,  radiative  and  nonradiative  transitions  of  3d 
group  ions  (Fe,  Ni,  Co)  in  II- VI  crystals.  (Subcontractor  -  Prokhorov  General  Physics 
Institute  (GPI),  RAS  business  representative). 

Task  3.  Fabricate  25  ZnSe  samples  having  a  concentration  series  of  Ni,  Co,  Fe,  Ni-Co,  and  Co- 
Fe  concentrations  ranging  from  10  to  5 x  1 0  cm"  (5  samples  for  each  dopant).  (UAB) 

Task  4.  Perform  temperature  dependent  absorption,  photo  luminescence  (PL),  time  resolved  PL, 
and  microscopic  characterization  of  fabricated  crystals.  (UAB) 

Task  5.  Perform  decay  kinetic  analysis  of  optical  energy  migration  and  concentration  quenching 
in  fabricated  samples.  (UAB) 
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Task  6.  Theoretically  analyze  the  optical  excitation  energy  transfer  schemes,  sensitization  as 
well  as  up  and  down  energy  conversion.  Evaluate  the  multipolarity,  dimensionality,  macro¬ 
efficiency  and  micro-efficiency.  Derive  an  analytical  solution  and  perfonn  a  Monte  Carlo 
simulation.  (GPI) 

Task  7.  Using  information  from  tasks  1-6  fabricate  and  optically  characterize  the  following: 

a)  Co:ZnSe  gain  elements  optimal  for  3-4  pm  lasing  under  Er- liber  laser  excitation.  (UAB) 

b)  Co:ZnSe  gain  elements  optimal  for  3-4  pm  lasing  under  808  mn  diode  laser  excitation. 
(UAB) 

c)  Co:ZnSe  gain  elements  optimal  for  3-4  pm  lasing  under  gain-switched  alexandrite  laser 
excitation.  (UAB) 

Task  8.  Using  information  from  tasks  1-6  fabricate  and  optically  characterize: 

•2d-  2d- 

a)  Ni"  co-doped  Co"  :ZnSe  gain  elements  optimal  for  3-4  pm  lasing  under  Yb-fiber  laser 
excitation.  (UAB) 

•2d-  2d- 

b)  Ni"  co-doped  Co"  :ZnSe  gain  elements  optimal  for  3-4  pm  lasing  under  gain  switched 
Nd:YAG  laser  excitation.  (UAB) 

Task  9.  Using  information  from  tasks  1-6  fabricate  and  optically  characterize: 

2d- 

a)  Co  co-doped  Fe:ZnSe  gain  elements  optimal  for  3. 8-5. 2  pm  lasing  under  Er-fiber  laser 
excitation.  (UAB) 

2d_ 

b)  Co  co-doped  Fe:ZnSe  gain  elements  optimal  for  3. 8-5. 2  pm  lasing  under  808  mn  diode 
laser  excitation.  (UAB) 

2d_ 

c)  Co  co-doped  Fe:ZnSe  gain  elements  optimal  for  3. 8-5. 2  pm  lasing  under  alexandrite 
laser  excitation.  (UAB) 

Task  10.  Laser  characterization  of  gain  elements  fabricated  in  tasks  7-9.  (UAB  and  Kirtland 

AFB) 

In  the  final  report,  the  results  of  our  research  are  summarized  in  the  following  chapters: 

■  Chapter  2:  Fabrication  of  II- VI  crystals  with  a  concentration  series  of  dopants  (tasks: 

1,3, 7, 8, 9) 

■  Chapter  3:  Theoretical  analysis  of  energy  level  diagrams,  radiative  and  nonradiative 
transitions  of  3d  group  ions  (Fe,  Ni,  Co)  in  II-VI  crystals  (tasks:  2,6) 

■  Chapter  4:  Spectroscopic  characterization  of  TM  doped  samples  (tasks:  4, 5, 7, 8, 9) 

■  Chapter  5:  Laser  characterization  of  TM  doped  samples  (tasks:  10)  Fabrication  of  II-VI 
crystals  with  concentration  series  of  dopants 

2.1  Fabrication  of  Ni2+;  Fe2+;  and  Co2+  doped  ZnSe/ZnS  samples 

Post-growth  thennal  diffusion  from  a  metal  film  on  the  surface  of  a  crystal  was  used  to  fabricate 
Co,  Ni,  Fe  and  Cr  co-doped  ZnSe/ZnS  samples  for  this  project.  This  technique  utilizes  thermally 
activated  diffusion  of  transition  metal  ions  into  II-VI  crystals.  The  typical  experimental 

procedure  includes  the  following  technological  operations: 
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•  Blank  crystal  preparation 

•  Deposition  of  the  metal  film  on  the  crystal  surface 

•  Sealing  of  a  sample  into  a  quartz  ample 

•  Annealing  of  the  sample  in  an  oven  at  high  temperature  (800-1200°C) 

The  average  TM  concentration  in  a  sample  can  be  controlled  by  the  specific  thickness  of  the  TM 
thin  film  deposited  on  the  crystal  surface.  The  thicknesses  of  a  TM  thin  film  required  for  a 
concentration  of  Njm  =  10 19  cm'3  in  a  1  mm  thick  crystal  are  presented  in  the  Table  1.  The 
density  associated  with  bulk  metals  was  used  for  these  calculations. 


Table  1.  Thickness,  h,  of  the  transition  metal  film  required  for  fabrication  of  a  1  mm 
sample  with  an  average  concentration  of  NTm  =  1019  cm'3 


Metal 

ft(nm) 

Cr 

120 

Fe 

118 

Co 

110 

Ni 

110 

Values  for  the  diffusion  coefficients  of  TM  ions  in  II- VI  hosts  reported  by  several  research 
groups  were  not  similar  with  differences  up  to  several  orders  of  magnitude.  These  differences 
may  have  resulted  from  several  parameters  which  are  difficult  to  control.  These  parameters 
include:  the  polycrystalline  structure  of  the  II-VI  host;  the  surface  crystal  preparation,  as  well  as 
sublimation  of  the  metal  and  host  during  thermal  diffusion.  However,  based  on  data  reported  on 
the  thermal  diffusion  coefficients  for  different  transition  metal  ions  under  similar  experimental 
conditions,  one  could  consider  the  following  to  be  the  rank  of  the  TM  diffusion  coefficients  in  a 
ZnSe  host  at  T  =  946°C  from  lowest  to  highest: 

D(Fe )  <  D(Co)  <  D(Ni )  <  D(Cr) 

IPG  Photonics  Corporation  has  developed  a  commercial  fabrication  procedure  for  quantitative 
(accuracy  of  pre-assigned  concentration  of  dopant  is  better  than  3%),  fast  thermo-diffusion  of 
chromium  ions  in  ZnSe  crystals  with  suppressed  sublimation  of  the  anion  and  cation  sublattices. 
The  resulting  fabricated  crystals  are  uniformly  doped  through  a  thickness  of  up  to  7  mm  and 
feature  a  low  scattering  loss  of  1-2  %  per  cm  for  samples  with  a  Cr  concentration  of  5  x  10  cm'  . 
The  diffusion  procedures  for  other  TM  ions  were  based  on  the  one  developed  for  chromium. 

Metal  targets  for  deposition  of  thin  films  having  purity  better  than  99.95%  were  purchased  from 
the  Kurt  J.  Lesker  Company.  Transition  metal  films  were  deposited  by  either  thermal 
evaporation  or  magnetron  spattering  in  a  vacuum  chamber.  Thennal  diffusion  was  then  carried 
out  in  sealed  vacuumed  (~10'5  Torr)  ampoules  at  temperatures  between  900-1 100°C  over  a 
period  of  7-20  days.  The  most  reliable  and  controllable  results  for  deposition  of  the  metal  film 
were  obtained  with  the  use  of  magnetron  sputtering  (see  Figure.  1). 
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Figure  1.  Plasma  Sputtering  System  at  IPG  Photonics  Corporation 


The  thin  film  deposition  rate  and  the  substrate  temperature  have  been  optimized  by  choosing 
parameters  which  provide  the  most  uniform  thin  film  thickness  according  to  Atomic  Force 
Microscopy  (AFM)  data,  and  the  best  film  adhesion  determined  by  the  stability  of  the  film 
during  rapid  heating  of  the  crystal  under  vacuum  conditions.  The  thickness  deposition  monitor 
has  been  calibrated  with  Atomic  Force  Microscopy  (AFM,  Park  Systems)  to  enable  a  precise 
determination  of  the  thickness  of  the  TM  film  on  ZnSe.  For  the  deposition  monitor  calibration, 
several  microscopic  slides  having  different  thicknesses  of  deposited  TM  films  were  prepared. 

Their  thicknesses  were  measured  by  the  AFM.  As  shown  in  Figure  2,  the  thickness  of  the  iron 
films  in  the  left  and  the  right  samples  were  estimated  to  be  1020  and  1305  nm,  respectively. 


Figure  2.  Topography  of  iron  films  having  different  thicknesses  as  measured  by  AFM  (each 

grid  represents  500  nm) 
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In  addition  to  the  calibration  of  thin  film  deposition  by  AFM,  the  percentage  transmission  of 
deposited  metal  films  was  also  measured.  Figure  3  shows  the  percentage  transmission  of  a 
cobalt  thin  film  deposited  on  a  glass  slide  measured  at  a  wavelength  of  1000  mn  as  a  function  of 
the  film  thickness.  As  one  can  see  from  Figure  3,  this  simple  method  can  be  used  to  determine 
the  thickness  of  a  film  lying  between  5-125  nm,  post  deposition. 


Cobalt  calibration  Cobalt  calibration 


Expected  film  thickness,  nm  Film  thickness,  nm 


Figure  3.  A)  Calibration  curve  for  the  deposition  monitor  used  for  cobalt  deposition 
during  magnetron  spattering  at  1000  nm;  B)  percentage  transmission  of  a  cobalt  thin  film 
deposited  on  a  microscopic  glass  slide  vs  thin-film  thickness. 


The  fabricated  samples  were  then  sealed  in  separate  quartz  ampoules  at  10'5  Torr  and  annealed  at 
900-1000  C  for  7-24  days.  Figure  4  shows  a  photo  of  cobalt  doped  ZnSe  samples  annealed  for  7 
days  at  950°C.  For  spectroscopic  characterization,  a  number  of  ZnSe  samples  having  a 
concentration  series  of  Ni,  Co,  and/or  Fe,  ranging  from  10  to  5><10  cm'  were  prepared. 
Several  samples  having  different  iron  concentrations  were  fabricated  from  the  initial 
Co:ZnS/Co:ZnSe  crystals.  Study  of  the  luminescence  quenching  of  cobalt  (nickel) 
photoluminescence  in  iron  co-doped  crystals  should  provide  information  on  the  energy  transfer 
rate  from  the  Co  (Ni)  co-dopants  to  the  Fe  ions. 


Figure  4.  Cobalt  doped  ZnSe  (7x3.1x50  mm3)  samples  after  annealing  for  7  days  at  950°C. 
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2.2  Fabrication  of  Ni2+  and  Co2+  doped  ZnSe/ZnS  samples 

Initially,  ZnS  and  ZnSe  samples  were  doped  with  the  cobalt  impurity  by  using  a  thermo- 
diffusion  method  which  is  described  in  detail  in  a  previous  report  [7].  The  Co:ZnSe/ZnS  samples 
were  then  co-doped  with  nickel.  The  nickel  co-doped  samples  were  also  fabricated  by  thermal 
diffusion  from  a  metal  film  deposited  by  the  Plasma  Sputtering  System  at  IPG  Photonics 
Corporation.  A  summary  of  all  the  fabricated  samples  is  shown  in  Table  2  and  pictures  of  the 
samples  811-60/62/63/64/65  are  shown  in  Figure  5  below. 


Table  2.  Summary  of  fabricated  Ni:Co:Fe:ZnS/Se  samples 


Crystal 

# 

Mat. 

Thicknes 
s  [cm] 

Width 

[cm] 

Length 

[cm] 

Concentration 
calculated  (Ni) 
[cm'3] 

Ann.time 
(at  950°C 
for  Ni) 

Concentration 
calculated  (Co) 
[cm'3] 

811-59 

ZnS 

0.33 

0.64 

5.7 

3.00X1018 

10  days 

l.OOxlO18 

811-60 

ZnS 

0.33 

0.64 

5.74 

3-OOxlO18 

10  days 

3.00xl018 

811-61 

ZnS 

0.33 

0.64 

3.15 

3-OOxlO18 

10  days 

3.00xl0iy 

811-62 

ZnS 

0.32 

0.64 

5.73 

3-OOxlO18 

10  days 

l.OOxlO20 

811-63 

ZnS 

0.33 

0.64 

2.74 

3.00xl019 

10  days 

3.00xl019 

811-64 

ZnSe 

0.33 

0.75 

4.3 

3.00xl018 

10  days 

l.OOxlO18 

811-65 

ZnSe 

0.3 

0.75 

4.2 

3-OOxlO18 

10  days 

3.00xl018 

811-66 

ZnSe 

0.35 

0.75 

2.34 

3.00xl018 

10  days 

3.00xl019 

811-67 

ZnSe 

0.33 

0.75 

4.33 

3.00xl018 

10  days 

l.OOxlO20 

811-68 

ZnSe 

0.33 

0.74 

2.26 

3.00xl0iy 

10  days 

3.00xl019 

Crystal 

# 

Mat. 

Thicknes 
s  [cm] 

Width 

[cm] 

Length 

[cm] 

Concentration 
calculated 
(Co)  [cm'3] 

Concentration 
calculated  (Fe) 
[cm'3] 

811-la 

ZnS 

0.25 

0.55 

0.96 

5.560xl017 

811-lb 

ZnS 

0.25 

0.55 

0.96 

6.007xlOr/ 

811-lc 

ZnS 

0.25 

0.55 

0.96 

8.261xl017 

811-ld 

ZnS 

0.25 

0.55 

0.96 

7.824xl017 

811-le 

ZnS 

0.28 

0.24 

0.55 

2.908xl017 

8 1 1  -If 

ZnS 

0.28 

0.24 

0.55 

1.596xl017 

811-lg 

ZnS 

0.28 

0.24 

0.55 

2.674xl017 

811-lh 

ZnS 

0.28 

0.24 

0.55 

1.846xl017 
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81 1-lk 


8.826x10 


811  -44a 
811  -44b 
811  -44c 
81 l-14d 
81 l-14e 
8 1 1  -1 5a 
81 l-15b 
81 l-43a 
81 l-43b 
811 -43c 
8 1 1  -1 7a 
81 l-17b 
8 1 1  -1 7c 
81 1 -1 7d 
811-29 
811-28 
81  l-58a 
8 1  1  -58b 
81 1 -57a 
81  l-57b 


Crystal  Mat.  Thick 

#  [cm] 

FI  109  ZnSe  L0 

A2109  ZnSe  082 

A1309  ZnSe  071 

J3110  ZnS  0.84 

J3210  ZnS  0.84 

J3310  ZnS  0.84 

J3510  ZnS  0.84 

J2110  ZnSe  0.77 


Width  Length 
[cm]  [cm] 


1.949xlOiy 
1.885xlOiy 
1.943xlOiy 
Brewster  1.781xl019 
Brewster  1.760xl019 
Brewster  2.017xl019 
Brewster  1.884xl019 
1.841xl019 
1.806xl019 
1.954xl019 
1.388X101' 
1.407xl018 
9.055xl017 
1.105xl018 
Brewster  2.150xl019 
Brewster  2.080xl019 
1.307xl019 
1.360xl019 
1.312xl019 
1.325xl019 


2.86x10  8 
2.86xl018 
3.66xl018 


2.88xl018 

2.76xl018 

8.49xl018 


3.14x10 

4.56x10* 

2.82x10* 

5.82x10* 


Concentration  Ann.time  Concentration 
calculated(Co)  forFe@  calculated(Fe) 


6xl018 


1000C 

17days 
17days 
17days 
66  days 
66  days 
66  days 
66  days 
16  days 


112x10* 

O.lxlO18 

38xl0*8 

l.lxlO18 

3-OxlO'8 

5.0xl018 

23xl018 
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J4110 

ZnS 

1.0 

0.54 

0.12 

6xl018 

16days 

0 

J4110 

ZnS 

1.0 

0.54 

0.12 

6xl018 

16  days 

0 

J4210 

ZnS 

1.0 

0.54 

0.12 

6xl018 

36  days 

0.5><  1 01S 

J4310 

ZnS 

1.0 

0.54 

0.12 

6xl018 

36days 

7.0xl0ls 

J4410 

ZnS 

1.0 

0.54 

0.13 

6xl0ls 

36days 

10xl018 

J2110 

ZnSe 

0.77 

0.75 

0.21 

6xl018 

16days 

0 

J2310 

ZnSe 

0.77 

0.75 

0.2 

6xl018 

36  days 

8.5xl0ls 

J2410 

ZnSe 

0.77 

0.75 

0.2 

6xl018 

36  days 

12.4xl018 

J2510 

ZnSe 

0.77 

0.75 

0.2 

6xl018 

36  days 

18.7xl018 

Figure  5.  Nickel  -  Cobalt  co-doped  ZnSe  and  ZnS  samples 


2.3  Fabrication  of  Fe2+  and  Co2+  doped  ZnSe/ZnS  samples 

Co-doped  Co:Fe:ZnSe  and  Co:Fe:ZnS  polycrystals  having  a  dimension  of  10x10x1.2  mm  were 
fabricated  using  a  two-stage  thermo-diffusion  method.  The  initial  thermal  diffusion  of  cobalt 
ions  was  perfonned  in  a  sealed  quartz  ampoule  at  10"5  Torr  at  an  annealing  temperature  of 
1000°C  for  16  days.  The  thermal  diffusion  of  iron  ions  was  performed  at  the  same  temperature 
for  20  days.  The  diffusion  procedure  was  optimized  to  provide  for  a  uniform  distribution  of  iron 
over  the  crystal  thickness  (~1  mm). 

A  controlled  thermal  diffusion  of  iron  ions  was  used  to  fabricate  co-doped  Co:Fe:ZnSe  and 
Co:Fe:ZnS  samples  with  an  Fe  concentration  range  of  8.5-19x10ls  cm'3  and  0.5-10x1018  cm"3, 
respectively,  while  the  cobalt  concentration  was  maintained  at  the  6.4x10  cm"  level  for  all 
crystals. 


10 

Approved  for  public  release:  distribution  unlimited. 


3.0  THEORETICAL  ANALYSIS  OF  ENERGY  LEVEL  DIAGRAMS,  RADIATIVE 
AND  NONRADIATIVE  TRANSITIONS  OF  3D  GROUP  IONS  (FE,  NI,  CO)  IN  II-VI 
CRYSTALS 


We  are  interested  in  the  energy  states  that  fall  within  the  forbidden  zone  of  II-VI  materials  since 
they  can  be  pumped  and  manifest  themselves.  The  purpose  of  the  modeling  in  this  section  is  to 
delineate  the  desired  energy  states  (that  fit  in  the  forbidden  zone)  and  the  other  energy  states  that 
cannot  be  pumped. 

3.1  Spectra  of  iron  group  ions  in  tetrahedral  semiconductors  of  the  AnB' 1  type 

Usually,  the  spectral  properties  of  the  dopant  ions  of  the  iron  group  are  investigated  in  dielectric 
crystals  with  a  high  crystal  field  potential  and  a  large  splitting  of  energy  levels.  As  a  result, 
transitions  between  the  primary  levels  of  absorption  and  luminescence  lay  in  the  visible  or  near 
infrared  wavelength  range. 

The  future  increase  in  the  data  traffic  of  telecommunication  optical  fiber  networks  calls  for  the 
creation  of  lasers  and  wideband  amplifiers  in  new  spectral  ranges.  In  particular,  for  the  second 
and  third  telecommunication  transmission  windows  (around  1.3  and  1.5  pm),  there  are  no 
effective  active  sources  that  are  compatible  with  common  silica-based  telecommunication  fibers. 
The  advancement  in  the  area  of  longer  wavelengths  (1.7-3. 5  pm  range)  is  limited  by  the 
absorption  properties  of  quartz  glass.  This  drives  the  utilization  of  a  different  type  of 
environment  such  as  telluride  (TeCfi)  or  AgBr  and  AgCl  materials  for  the  middle  infrared.  In 
addition,  powerful  light  sources  in  the  1.7-3. 5  pm  spectral  range  can  be  frequency  doubled  into 
the  0.85-1.75  pm  spectral  range  for  scientific,  technological  and  medical  purposes. 

The  basic  problem  of  modern  quantum  and  nano  electronics  is  the  creation  of  lasers  and 
amplifiers  in  semiconductor  materials.  Creation  of  such  devices  would  involve  a  symbiosis  of 
the  spectral  properties  of  the  active  centers  and  the  basic  semiconductor  material.  Taken  together 
with  the  ability  of  a  semiconductor  material  to  carry  a  current,  new  possibilities  in  both  quantum 
and  nano  electronics  can  be  opened  up. 

Among  the  large  set  of  possible  semiconductor  materials,  those  materials  having  an  ionic  bond 
and  a  wide  forbidden  zone  of  the  type  A11  BVI  are  of  specific  interest.  Thus,  it  is  necessary  to 
verify  that  the  absorption  levels  of  an  active  center  align  with  the  width  of  the  forbidden  zone  so 
that  sufficient  pumping  can  be  obtained.  The  intricate  problem  of  the  interaction  between  the 
electrons  of  the  active  center  with  the  electrons  of  the  valence  zone  and  the  conduction  band 
demands  special  research. 

3.1.1  Red  shift  of  iron  group  ion  spectra  in  tetrahedral  semiconductors  of  the  AnBVI  type. 

In  crystals,  ions  belonging  to  the  iron  group  may  occupy  the  coordinate  centers  of  polyhedrons 
having  orbits  with  coordinate  numbers  k  =  4,  6,  and  8.  The  crystal  field  parameter  Dq  (as  well  as 
its  sign)  depends  on  the  number  of  coordinates  associated  with  the  orbits: 

Dq{6)>  0,  (1) 

Dq(4)=-^Dq{6)  <2) 
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Semiconductor  crystals  having  tetrahedral  orbits  with  k  =  4  have  a  weaker  crystal  field  (CF) 
potential  than  dielectric  crystals  having  octahedral  orbits  with  k  =  6. 

In  crystals  of  the  AnBVI  type,  ions  of  the  iron  group  are  bivalent  and  have  a  smaller  CF  potential 
and  Dq  value  than  they  would  have  if  they  were  in  the  trivalent  state  in  an  octahedral  dielectric 
crystal.  Thus,  in  AnBVI  semiconductor  crystals  doped  with  ions  of  the  iron  group  having 
tetrahedral  orbits,  a  “red  shift”  should  be  exhibited  with  the  transition  energies  of  the  absorption 
and  fluorescence  spectra  of  the  ions  of  iron  group  characterized  by  lower  frequencies. 

The  parameter  Dq,  which  depends  on  the  charge  of  an  electron,  has  a  different  sign  for  an  iso 
electron-hole  configuration.  Thus,  the  Dq  parameter  for  the  hole  configuration  of  ions  of  the  iron 
group  in  tetrahedral  semiconductor  crystals  is  positive.  It  is  also  positive  for  the  case  of  an  iso¬ 
electron  configuration  of  ions  in  octahedral  crystals  of  a  dielectric: 

Dq(d'^",k  =  4)>0  Dq(d",k  =  6)>0  <3) 

The  parameters  of  the  electronic  electrostatic  interaction  depend  on  the  charge  squared  and  have 
the  same  sign  for  both  the  electron  and  hole  configurations.  The  spin-orbital  interaction 
parameter  is  linearly  proportional  to  the  sign  of  electronic  charge,  so  it  changes  sign  when  the 
hole  configuration  is  changed.  In  summary,  the  parameters  of  crystal  field  (Dq),  the  electrostatic 
interaction  (F2  and  F4)  and  the  spin-orbital  interaction  (cvi)  of  ions  of  the  iron  group  in 
semiconductor  crystals  of  the  AUBVI  type  should  have  the  following  signs: 

Dq  >  0,  F2  >  0,  F4  >  0,  <  0.  (4) 


3.1.2  Piezoelectric  effects  in  tetrahedral  semiconductors  of  the  AnBVI  type. 

The  tetrahedral  symmetry  of  the  group  Tti  (no  center  of  symmetry)  defines  the  existence  of  a 
piezoelectric  effect  in  semiconductor  materials.  This  allows  one  to  use  an  external  electric  field 
to  control  the  spectral  properties  of  the  laser  transitions.  It  is  difficult  to  expect  a  large  influence 
from  an  external  electric  field  on  the  frequencies  of  laser  transitions.  However,  mechanical 
deformations  caused  by  the  piezoeffect  can  generate  odd  harmonics  of  a  crystal  field  in  the 
environment  of  an  active  center  thereby  changing  the  intensity  of  the  transitions.  The  result  will 
be  a  redistribution  of  the  transition  intensities  on  the  Stark  levels  that  will  depend  on  the  size  and 
direction  of  the  applied  electric  field. 

3.2  The  peculiarities  of  spectra  of  iron  group  ions  in  AnBVI  semiconductor  crystals 

The  spectra  of  iron  group  ions  are  usually  analyzed  by  means  of  Tanabe-Sugano  (TS)  diagrams. 
However,  they  have  limited  application  for  tetrahedral  centers  in  semiconductor  crystals. 

TS  theory  is  based  on  the  strong  (intennediate)  field  case,  where  the  energy  of  interaction  with 
the  crystal  neighborhood  is  greater  than  the  electrostatic  interaction  between  electrons.  In  this 
case,  the  3d  energy  configuration  of  iron  group  ions  in  crystal  fields  having  Oh  group  symmetry 
is  split  into  e  and  t2  levels,  3dn  =>  emt"=m .  In  semiconductor  crystals  of  the  AnBVI  type,  the 
crystal  field  interaction  potential  is  smaller  than  the  electrostatic  interaction  and  there  is  therefore 
little  use  in  considering  configurations  on  the  basis  of  strong  CF  state  functions.  Using  this  basis, 
the  quantum  numbers  of  a  free  atom  and  the  representations  of  the  orthogonal  rotation  symmetry 
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group  O3  are  meaningless.  In  this  case,  quantum  numbers  are  replaced  by  the  indices  of  the 
irreducible  representations  of  a  point  symmetry  group  for  which  there  is  no  possibility  of 
defining  exact  quantum  numbers  and  perfonning  analytical  calculations  algebraically.  As  a 
result,  the  powerful  tool  of  angular  momentum  summation  theory  is  rendered  useless. 

TS  matrices  make  use  of  the  Racah  parameters  of  the  electrostatic  interaction  which  are  linear 
combinations  of  the  Slater  parameters.  The  Racah  parameters  arise  from  the  internal  symmetry 
of  the  free  atom  states.  They  facilitate  the  diagonalization  of  the  matrices  of  the  electrostatic 
interaction.  In  TS  theory,  combinations  of  the  Racah  parameters  are  used  in  the  basis  of  the 
strong  CF  wavefunctions.  When  studying  the  dependence  of  the  electrostatic  interaction 
parameters  on  covalent  effects  in  semiconductor  crystals,  it  is  better  to  use  the  Slater  parameters 
which  have  an  explicit  form  for  the  radial  wave  function  of  the  atom  rather  than  a  combination  of 
the  radial  wave  functions  of  various  orders. 

In  tetrahedral  semiconductors,  the  symmetry  group  is  of  the  type  Td  or  T.  While  even  parts  of  the 
crystal  potential  for  the  groups  Oh  and  (Td,  T)  do  not  change,  the  types  of  irreducible 
representations  do.  Non-diagonal  elements  (absent  in  crystal  fields  of  Oh  symmetry)  appear  as 
shown  in  [8],  In  addition,  odd  terms  (of  the  third  order)  of  the  crystal  potential  combine  with  an 
odd  electro-dipole  operator  to  form  an  even  operator  allowing  an  electro-dipole  transition.  That 
is  why  TS  matrices  and  diagrams  are  of  little  use  for  evaluating  oscillator  strengths. 

Another  point  is  that  relativistic  effects  must  be  taken  into  account  for  ions  of  the  3dn 
configuration  in  AUBVI  semiconductors.  The  effects  are  manifested  by  the  spin-orbit  (SO) 
interaction,  which  can  be  as  strong  as  the  CF  interaction  in  the  semiconductor  crystal  field.  The 
SO  interaction  splits  the  crystal  field  levels  and  is  seen  in  the  band  structure  of  transition  metals. 
When  considering  the  spin-orbit  interaction,  one  needs  to  employ  the  double  symmetry  group 
state  classification  which  cannot  be  accomplished  using  a  basis  comprised  of  the  strong  field 
wave  functions  of  the  TS  theory. 

In  addition,  the  actual  symmetry  of  a  crystal  field  can  be  described  by  lower  symmetry  groups  of 
the  trigonal  or  rhombic  type.  This  requires  additional  CF  terms  that  are  outside  the  framework  of 
the  TS  matrices. 

All  of  the  above  calls  for  the  development  of  a  new  theory  and  the  construction  of  new  matrices 
for  the  3d11  configuration  in  order  to  take  into  account  all  of  the  effects.  To  address  this  need,  a 
theory  and  concrete  fonnulas  to  enable  an  analytical  calculation  of  the  matrix  elements  of  the 
electrostatic  interaction,  spin-orbit  interaction,  and  crystal  field  interactions  in  the  basis  set  of 
wave-functions  with  the  quantum  numbers  /  =  2,  SLJnij  will  be  developed.  We  will  work  to 
enable  the  calculation  of  all  matrices  in  one  mathematical  program  without  referring  to  any 
external  data  for  the  3j,  6j,  and  9j  genealogical  coefficients.  Thus,  the  program  will  compute  all 
of  the  matrices  of  the  3d11  configuration  states. 

3.3.  State  classification  of  the  electronic  configuration  3dn 

Classification  of  the  states  of  the  3d11  electron  configuration  is  based  on  reduction  of  the 
irreducible  representations  of  the  following  groups: 
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^4/+2  13  3^2  X  ^2/+l  — 5  SU  2  X  R3 

M  3  [S]x[w]  3  [S]  X  [L]-  (5) 

where  U4/+2  is  a  unitary  group  having  a  dimension  of  4Z+2  that  represents  all  the  states  of  the 
anti-symmetric  representation  {/n}  which  describe  an  electron  in  the  /”  configuration,  W  is  a 
representation  of  the  orthogonal  group  R2/+1  having  a  dimension  of  2/+1,  and  S  and  L  are  the  spin 
and  orbital  quantum  numbers  representing  the  spin  and  orbital  groups  SU2  and  R3,  respectively. 
Ions  in  the  dn  configuration  with  1  =  2,  are  represented  by  the  unitary  group  U10  and  the 
orthogonal  group  R5  of  five-dimensional  space.  As  an  example,  one-electron  in  the  dl 
configuration  with  /  =  2  is  characterized  by  the  following  series: 

Ul0  id  SU2  x  R5  id  SU2  x  R3 

{1}  =>  [l]x[l]  =>  [l]x[2]=2  Z>  (6) 

All  states  of  any  dn  configuration  can  therefore  be  obtained  by  reducing  one  anti- symmetric  state 
{ 1''}  of  the  U10  group  as  in  (6). 

3.3.1  States  of  the  total  angular  momentum  J  of  the  spin-orbital  interaction. 

Due  to  the  spin  -  orbital  interaction,  the  spin  S  and  the  angular  momentum  L  are  mixed  with  a 
total  angular  momentum  J  being  obtained.  The  spin  S  and  the  orbital  momentum  L  of  a  given 
electronic  state  give  the  total  angular  momentum  J  according  to  the  following  rules  for  the 
multiplication  of  the  irreducible  representation  of  the  orthogonal  rotation  group  R3:  LxS  =  L+S, 
L+S-l . . . |L— S|.  For  example,  the  3F  state  has  the  following  multiplets:  S=l,  L=3  and  J=4,  3,  and 
2,  which  results  in  the  F4,  F3,  F2  spin  -  orbit  states.  Tables  of  the  states  of  the  three  electron 
configurations  dn  where  n  =  6,  7,  and  8  with  the  values  of  J  are  presented  below. 

3.3.2  States  of  Fe2+  (3d6  electronic  configuration). 

2_|_  /r 

The  Fe  ion  has  a  3d  electronic  configuration  and  all  its  states  are  detennined  by  the  following 
representations: 
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3  SI  I,  x  A',  =;  SU,  x  R, 

{l‘}  =>[2]x[l0]  => [2]x[2]=5  D 
[l]x[ll]  =[l]x[l]  ='/><'> 

=>[l]x[3]=’  Fm 

[l]x[2l]=[l]x[l]=»  Pm 

[l]x  [2]  — 3  D 

[l]x [3]='  Fm 
[l]x  [4]  =3  G 
[l]x[5]=  « 

[0]  x  [OO]  =>  [o]x  [0]  =‘  S(1) 

[0]  X  [22]  =3  [o]x  [0]  ='  S<2) 

[0]  x  [20]  Z)  [0]  x  [2]  =‘  Dm 
[ojx  [22]  =3  [0]  X  [2]  ='  D(2) 

[0]x[3]='  F 

[0] x  [20]  Z)  [ojx  [4]='  G(1) 

[0]  x  [22]  Z)  [0]  x  [4]  ='  G(2)  •  (7) 

[0]  x  [6] = 1  / 


The  resulting  states  are  designated  as  2S+1l(p).  The  upper-left  index  denotes  the  spin  (S) 
degeneration  and  is  represented  by  number  2S+1.  The  upper-right  index  (p)  =  (1)  or  (2) 
distinguishes  duplicate  states  having  the  same  quantum  numbers.  The  orbital  quantum  numbers 
are  denoted  by  the  letter  L  with  the  following  correspondence: 

S  P  D  F  G  H  I 

Quantum  number  LO  1  2  3  4  5  6 

OClt 

The  orbital  degeneration  is  equal  to  2L+1.  The  degeneration  numbers  of  the  L  states  are  equal 
to  (2S+1)(2L+1).  Finally,  the  states  and  their  degeneration  numbers  for  the  3d6  electron 
configuration  of  Fe2+  are  given  in  Table  3. 

Table  3.  States  and  their  degeneration  numbers  for  the  3d6  electronic  configuration  of  Fe2+ 


State  3d6 

5D 

3p(l) 

3p(2)’ 

3Fm 

3p(2> 

3D 

3G 

3H 

is<2) 

1D(1), 

id(2) 

’F 

1G(1), 

1G,2) 

’I 

L-degeneration 

5 

3 

7 

5 

9 

11 

1 

5 

7 

9 

13 

Total 

degeneration 

25 

9x2 

21x2 

15 

27 

33 

1x2 

5x2 

7 

9x2 

13 

Angular 
momentum  J 

4,3, 

2,0 

2,  1,0 

4,3,2 

3,2, 

1 

5,  4, 

3 

6,  5, 

4 

0 

2 

3 

4 

6 

15 
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The  states  with  quantum  numbers  S  and  L  are  obtained  by  taking  into  account  only  the 
electrostatic  interaction  between  the  electrons. 

Graphically,  a  diagram  of  the  full  matrix  of  the  electronic  configurations  d6  is  presented  in  Figure 
6.  The  white  spaces  correspond  to  elements  of  the  matrix  which  are  zero  while  the  black  spaces 
correspond  to  different  nontrivial  values  which  depend  in  a  complex  way  on  the  parameters  of  all 
the  interaction  types. 


21Q 

200 


150 


E 

-i  ioo 


50 


1 

1  50  100  150  200 

row 

Figure  6.  Diagram  of  the  matrix  describing  the  3d6  electronic  configuration  of  Fe2+. 

The  total  number  of  states  is  210. 


Base  Matrix 


2 1  j 

3.3.3  States  of  Co  (3d  electron  configuration). 

Co2+  ion  has  a  3d7  electronic  configuration  and  all  its  states  are  determined  by  the  following 
representations: 


16 

Approved  for  public  release:  distribution  unlimited. 


U10  3  SU2  xS5d  SU2  x  R3 
{ l7  }  3  [3/2]  X  [l  l]  3  [3/2]  X  [l]  =4  P 
[3/2]  X  [ll]  3  [3/2  ]x  [3]  =4  F 

[l/2]  X  [2l]  3  [l/2]  X  [l]  =2  P 

[l/2] x [lO]  3[l/2]x[2]=2  Dx  (g) 

[l/2]  X  [21]  3  [1/2]  X  [2]  =2  D2 
[1/2]  X  [2l]  3  [1/2]  X  [3]  =2  F 
[1/2] X  [21]  3[l/2]x[4]=2  G 
[1/2]  X  [21]  3  [1/2]  x  [5]  =2  H 

The  states  and  the  degeneration  numbers  associated  with  the  3d7  electronic  configuration  of  Co2+ 
are  given  in  Table  4.  The  diagram  of  the  matrix  is  shown  in  Figure  7. 


Table  4.  States  and  the  degeneration  numbers  for  the  3d7  electronic  configuration  of  Co2+ 


State  3d7 

r 

2P 

2d2 

2F 

2G 

ZH 

L-degeneration 

3 

7 

3 

5 

7 

9 

11 

Total  degeneration 

12 

28 

6 

10x2 

14 

18 

22 

Angular 
momentum  J 

5/2,  3/2, 

1/2 

9/2,  7/2,  5/2, 
3/2 

3/2, 

1/2 

5/2,  3/2 

7/2, 

5/2 

9/2, 

7/2 

11/2, 

9/2 
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Base  Matrix 
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Figure  7.  Diagram  of  the  matrix  for  the  3d7  electron  configuration  of  Co2+. 
The  total  number  of  states  is  120. 


3.3.4  States  ofNi2+  ( 3ds  electronic  configuration). 

.9+  o 

Ni“  has  a  3d  configuration  and  the  states  are  determined  by  the  following  representations: 

Ul0  z>  SU2  x  R5  DSt/2xi?3 

{l*MlMll]=,[l]x[l]=> 

=>[l]*[3]=3F 

[0]x[00]=»[0]x[0]='5  •  <9> 

[0]x[20]  z^Ojx^1  D 
=>[0]x[4]=‘G 

8  .  .  9_j_ 

States  and  the  degeneration  numbers  for  the  3d  electronic  configuration  of  Ni“  are  given 
Table  5.  The  diagram  of  the  matrix  is  shown  in  Figure  8. 
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Table  5.  States  and  the  degeneration  numbers  for  the  3d8  electron  configuration  of  Ni2+ 


State  3d8 

*S 

‘d 

*G 

L-degeneration 

3 

7 

1 

5 

9 

Total  degeneration 

9 

21 

1 

5 

9 

Angular  momentum  J 

2,1,0 

4,3,2 

0 

2 

4 

Base  Matrix 


□  □ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 
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□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□  □ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 
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□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 
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□ 
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Figure  8.  Diagram  of  the  matrix  for  the  3d8  electronic  configuration  of  Ni2+. 

The  total  number  of  states  is  45. 


3.4  Theory 

We  will  be  dealing  with  the  3d"1  electronic  configuration  of  M  atoms  in  an  MB p-comp\cx.  The 
goal  is  to  calculate  the  main  interactions  of  this  molecular  complex  to  include  the  spin-orbit  and 
crystal  field  interactions  of  M  with  the  ligand  electrons.  The  Hamiltonian  of  these  interactions  is 
described  by: 

H=Vqq+HCF+Hso  (10) 
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where  Vqq  is  the  matrix  representing  the  electrostatic  interaction  between  electrons,  H.so  is  the 
spin-orbit  interaction  and  HCf  is  the  crystal  field  interaction. 

The  phenomenological  method  of  molecular  orbital  (MO)  (PMMO)  theory  involves  calculation 
of  the  energies  by  using  the  symmetry  of  the  MO  states  and  the  interaction  operators.  The 
parameters  are  then  determined  by  comparing  theory  with  the  experiment. 

Because  of  this,  matrix  elements  of  the  operator  defined  by  equation  (10)  may  be  calculated  with 
a  formula  invoking  the  full  rotation  group  matrix  elements.  Therefore,  the  problem  consists  of 
the  operation  of  one  or  two  electron  interaction  operators  on  a  double  (l\,  h )  shell  electronic 
configuration. 

Matrix  elements  of  the  electrostatic  and  exchange  interactions  in  the  form  of  two  electron 
operators  are  represented  by  [9,  10] 

fl;  jl;(l]h-)F  +/l(lli2yfi(V24)|74. 

+  (-l)S*t('l'2i>2)Cgt(/1/22)G2  +(-l)S«t(/1/2L,l)Cst(/1/2l)G1  (n) 

where  F2,  F4  and  G2,  G\  are  the  parameters  of  the  electrostatic  and  exchange  interactions, 
respectively,  with  the  coefficients  fk(l\hL,k),  gk{l\hL,k),  Cfk{l\hk),  and  Cgk{l\hk )  being  expressed 
in  terms  of  the  3 j  and  6j  coefficients  [11]  as: 


c*  (V2O  =  (- 1)'1*1'*'  (2/,  + 1)(2/,  + 1) 


(h  k  IA 

72 

k  l2) 

"0 

0 

0 

lo 

0 

0 

Cjl,l,k)  =  (21,  +1X24+1) 

V 


h 

0 


k  /,y 

0  0y 


(12) 


Matrix  elements  for  the  exchange  interactions  with  the  parameters  Gi  and  G2  exist  only  for  the  p 
and  d  states,  respectively.  The  matrix  elements  for  the  exchange  operator  in  eq.  (10),  last  two 
tenns  in  equation  (11),  are: 

Hex  =  S2(dp2+  S\(pP\, 

«2  0)  =  (- 1  f  St  (2, 0,2, 2  )Cgt  (2,0,2)  S,  (/-)=(-  1)S  St  (lAU)Cgt  (l,0,l)  (]3) 

The  matrix  elements  for  the  spin-orbit  interaction  operator  of  eq.  (10)  are  [9]: 

H  so{l\l2SlS2LxL2J  )  =  %dVd(lxl2SxS2LxL2J  ), 

V,  {lll2SlS2LlL2j)=  (-  ir2+1+w  72(2/2  +  1X2^  +  ®  +  1X2 L2  +  1X2 52  +  1)  .  (14) 

\h  h  A  lfl/2  1/2  Sj  1  fLj  5,  J  | 

1 U  1  12\{S2  1  1/2J1S2  L2  lj 
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The  interaction  with  the  CF  potential  is  defined  by  matrix  elements  whose  values  may  be 
calculated  from  the  following  fonnula: 


vgf{ili2SLlL2jlj2qlq2)-  £  A 


k,q 


kq 


J\  k  J  2 
<7l  <7  <72 


Cgf  (l\l2SL\L2J\J  2  k). 


(15) 


where 


Cgf(l1l2SLlL2JlJ2k)  =  {-l)h+h+s^  a/(2Li+1)(2/1+1X2L2+1)(2/2+1) 

1 2  A  MjA  A  Sl 
l2  /2  kj|y2  fi2  kj 

and  the  square  brackets 

,/j  k  J  7 

q\  ci  <?2 


J ,  k  J- 


J2J,  +i  1 <?i  q  <?2 


denote  the  Klebsh-Gordon  coefficient.  Note  that  14,/  is  equivalent  to  Hcf  in  equation  10. 

Usually,  in  the  case  of  semiconductors  of  the  type  AUBV ,  the  CF  symmetry  of  the  active  center  is 
described  by  the  Td  group  with  distortion  from  the  C3V  group  determining  the  set  of  CF 
parameters:  k,  q  =  2,  0;  4,  0;  4,  3;  4,  -3.  It  is  more  convenient,  though,  to  use  well-known  CF 
parameters  corresponding  to  the  pure  tetrahedral  ( Dq  parameter)  and  trigonal  (Di  parameters) 
components  of  the  CF  potential  [9]: 


A2Q  -  A2E  > 
A40  =  2 1  Dq 


715 

~6~ 


D4E > 


A44 


A4_4  -  A44 


(16) 


Finally,  the  interaction  operators  in  (10)  are  parameterized  with  the  following:  £4  as  the  spin-orbit 
constant  of  the  d(M)  electrons;  A20,  A40,  and  A44  as  the  parameters  of  interaction  with  the  CF 
potential  of  tetrahedral  (Td)  (Dq)  symmetry;  and  D2E  ,  D4E  as  the  low  symmetry  parameters  of 
the  CF  potential. 


3.4.1.  Formula  for  calculation  of  the  matrix  elements  of  the  crystal  field  potential. 

q=k 

{SLxJxmx  \VCF\SL2J2m2)  =  ^  Y,AkqJn(Ji’k,J2,mx,q,m2Y:rd{nx,n2,S,Lx,L2,Jx,J2,k) 

£=2,4  q=-k  , 

The  calculation  of  this  is  described  in  Appendix  1 .  Ov  is  a  name  for  a  computer  procedure 
calculating  a  part  of  the  matrix  element  of  the  CF  potential. 

3.4.2.  Formula  for  calculation  of  the  matrix  elements  of  the  electrostatic  interaction 

(, nl,SLlJlml  \Vqq\n2,SL2J2m2)  =  qq(nl,n2 ) 
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(18) 


The  calculation  of  this  is  described  in  Appendix  1 . 


3.4.3.  Formula  for  calculation  of  the  matrix  elements  of  the  spin-orbit  interaction 

(nx,SxLxJ xmx  I^sqI n2 , S 2L2J 2m2 )  — 

d{j xmx ,  J 2^2  (^1  ’  ^2  ’  ^1  ’  ^2  ’  *^l  ’  *^2  (  C  5  ^2  ’  5  F2  ’  *^1  ’  ^2  ’  (19) 

W11  is  the  notation  of  the  reduced  double  (11)  tensor  operator  in  the  orbital  and  spin  coordinates 
of  the  electron.  The  calculation  of  this  is  described  in  Appendix  1 . 

3.5.  Intensities  of  transitions  between  states  having  the  3d11  electronic  configuration 

Electric  dipole  (ED)  transitions  between  states  of  the  dn  configuration  are  forbidden  by  parity 
selection  rules.  However,  in  crystals,  electro-dipole  transitions  may  exist  due  to  the  admixture  of 
the  states  having  opposite  parity  due  to  the  odd  tenns  of  crystal  field  potential  (CFP). 

The  classification  of  the  odd  terms  of  the  CF  potential  may  be  done  using  the  following 
symmetry  group  reduction  series: 


03  3  °h  =>  C/  =>  C 3v  3  C4v  3  C2v 

[3]  3  A2u  3  Ai 

['It3]  3  F\u  3F23Aj 

[3]  =>  F\u  =>  F2  3  =>  Aj 

[3]  =>  f2u  ^  F\  3  =>  A  (20) 

The  four  terms  of  the  series  (20)  has  the  following  symmetry:  Td,  C3V,  C4V  and  C2w,  respectively. 
This  is  determined  by  reduction  of  the  odd  representations  of  O3  on  the  unity  representation  of 
the  corresponding  groups. 

The  odd  part  of  CFP  may  arise  from  the  odd  static  and  dynamic  components  of  the  crystal  field. 
The  static  component  which  has  the  representation  [3]  of  the  O3  group  depends  on  the  positions 
of  the  surrounding  atoms  and  the  second,  dynamic  component  having  the  representation  [1]  of 
the  O3  group  may  be  due  to  the  vibrations  of  the  atoms. 

For  case  of  Td  symmetry,  the  resulting  even  ED  operator  may  be  obtained  as  the  product  of  the 
odd  tenn  A2U  operator  of  the  CFP  [3]  and  the  operator  [1]  of  the  electro  dipole  transition  with  the 
result:  [3] x [  1]=  [4]+[3]+[2].  For  the  p  and  d  states,  only  representations  [2]  and  [4]  have  non¬ 
zero  matrix  elements.  The  same  representation  in  the  Td  group  is  determined  by  the  next  product: 
[l]Fiux[3]A2u=[2]F2g+[4]F2g  and  only  the  [2]F2  and  [4]F2  operators  of  Td  group  may  be  used  as 
representations  of  the  transition  operators.  Thus,  for  all  of  the  symmetry  groups  (20),  the 
transition  operator  is: 

[2],  [4]  F2(0)  — >  a  polarization, 

[2],  [4]  F2(±l)  — ►  71  polarization.  (21) 

The  resulting  intensity  operator  is  an  even  with  the  matrix  coinciding  with  the  crystal  field 
matrix. 
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For  the  case  of  tetrahedral  symmetry  of  the  Td  group,  there  exists  only  one  independent 
parameter  (3A2),  with  the  other  parameters  being  a  linear  combination  of  a  spherical  function 
with  the  parameters  Bkq: 

BkVq  =  (l^l u  ’3A2 u  \kT\F\uVl ,  A2ue\^ci\  ^F\uMl  )|  ^A2ue)  ■>  (22) 

where  the  subgroup  and  Klebsh-Gordon  coefficients  are  presented  in  [12]  and  the  ratio  between 
the  intensities  of  the  spherical  parameters  is: 

B4q  =~^B2q-  (23) 


Expressions  for  the  B20,  B21,  and  B2-1  polarization  parameters  in  terms  of  the  spherical  functions 
are  the  following: 


B20  ~ 


—  3  z2~r2 
16k  r2 


’ B21  ~ 


[*  z(x  +  iy) 


IT5  z(~x  +  iy ) 

V  8k  r2 


(24) 


Intensities  of  the  MO  transitions  may  be  estimated  using  a  matrix  representing  the  transition 
operators  followed  by  normalization  of  this  matrix  by  calculating  the  interaction  between  the 
states  of  the  MO  configuration. 

The  energy  level  values  or  eigenvalues  are  obtained  by  diagonalization  of  the  matrix  Crp 
represented  by  the  operator  by:  E(SL)=MxCrpxM'1  where  M  is  the  basis  function 
transformation. 


The  nonnalized  expression  for  the  transition  intensities  will  be: 

IiSUSL^MxEopXM'1,  (25) 

and  the  intensities  of  the  electric  dipole  transitions  between  the  states  may  be  calculated  by  the 
following  formula: 


2q^q2  )  —  B 


2<1 


J]  2  J2 

<?i  q  q2 


CRf{l\l2SL\  L2J\J22)+  B4q 


J 1  4  -G 
q\  q  q2 


C.f(l\l2SL\L2J\J24),  (26) 


where 

Cgf{lll2SLiL2JlJ2k)={- 1  )'1+,2+Sl+72  V(2  Lx  + 1)(2/1  +  1)(2L2  +  lX2/2  +  l) 
|/2  L;  /JjL;  Jx  s)k  =  24 

\l2  i2  k}\j2  l2  ky 


Intensities  of  the  magnetic  dipole  transitions,  which  can  be  obtained  by  operating  L+2S  on  the 
MO  states  of  equivalent  electrons,  may  be  calculated  from  the  following  formula: 


l(lll2SLlL2JlJ2qlq2) -  s(li,l2 ]S(lu  L2)A]c/ 


h 

q\ 


1 

q 


J  2 

q2 


CL+2S 


(l\l2SLxL2J  1/2 1)  > 


(27) 


where 
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JL+2S 


{lxl2SLxL2JxJ2 1)  =  (- l)L>+5+^  JLl(Ll+  lX2L1+l)  r  1 


■  2(-  l)Li+s+7i  75(5  +  1X25  +  1) 


5  /,  Lx 

J2  5  1 


3.5./.  Formula  for  calculation  of  the  transition  intensity  matrix  elements  for  transitions  in 
iron  group  ions. 

q=k 

(5Lj7 xmx  |t/yi  j2 1 SL2J  2m2 )  =  kqJ iSf  \ ->k,J 2,mx,q, m2  )Crd (/7j ,n2,S ,LX,L2,  J J 2 ,  &) 

k=2,4  q=-k 


The  transition  intensity  matrix  elements  are  calculated  in  the  same  manner  as  the  matrix  elements 
of  the  crystal  field  with  the  exception  of  using  the  intensity  parameters  Bkq. 


3.6.  Simulation  of  the  absorption  spectra  of  dn  electronics  configurations 

Because  the  dn  states  have  many  levels  which  are  clumped  into  a  few  bands  of  levels,  the 
absorption  spectra  of  the  transitions  are  characterized  by  broad  bands.  This  is  why  it  is  preferable 
to  represent  the  spectral  distribution  of  the  transition  between  the  ground  and  exited  states  as  a 
Gaussian  which  depends  on  the  wavelength.  Thus,  theoretical  simulation  of  the  absorption 
spectra  may  be  performed  using  the  following  fonnula  for  the  effective  cross  section  of  the 
absorption  transitions: 


70  n 

QU)  =  Z  %exP 

v=5  A 


A107 

X 


+  E0~Ev 


(28) 


where  E{]  and  Ev  are  the  energies  of  ground  and  excited  states,  respectively,  Sv  is  the  dispersion 
width  parameter,  and  Qv  is  the  effective  cross-section  in  cm'1. 


3.7.  Oscillator  strength  and  decay  time  of  dn  electronics  configurations 

It  is  reasonable  to  analyze  the  absorption  spectra  of  ions  of  the  iron  group  using  parameters 
which  are  independent  of  the  form  of  absorption  band.  These  parameters  are  the  oscillator 
strength  of  the  absorption  transition  and  the  fluorescence  lifetime  (or  decay  time). 

The  oscillator  strength/and  decay  time  r  are  given  by  following  formulas: 

^  3^jr  3  4  co  -2 

f/r  =  — ^ — - —  — z —  =  4.52cm  5 

8;r  cr0Tjj,  8k  cr0 
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c2e 

f.,=  4.52— - =  1.13*10“ 


V~T 


JJ' 


fj  ^ 

AJJ' 

500 


=  1.13*10 


■  jj' 


20000 

v  EJJ'  J 


'  JJ' 


(29) 


fjj'  -4.52  <„■  =4.52-/— 

Ejj'tJJ'  jrJjr 

•> 

where  E jj-  is  the  energy  gap  between  the  states  J  and  J’,  the  speed  of  light,  c,  and  the  electron 
radius,  ro,  are  equal  to  3.0><1010  cm/sec  and  2.82*  10"13  cm,  respectively.  Also,  the  decay  time, 
wavelength,  and  energy  of  the  transitions  are  represented  by  x(sec),  X  (cm),  and  E  (cm1), 
respectively.  The  factor  representing  the  local  field  of  the  impurity  atom  arising  from  the  local 
polarization  effect  is  given  by 


where  n  is  the  index  of  refraction. 

The  transition  probability  A  is  connected  with  the  effective  cross  section  Q  according  to  the 
following  formula: 

=  W  A/i \E2Qjr(E}lE ,  (31) 

2J  + 1  J 

where  E  is  frequency  of  the  electron  transition  in  cm"1  and  Qjj-  is  the  effective  cross-section. 


a  = 


9  n 


2  ,  o 

n  +  2 


(30) 


The  oscillation  strength  for  electro-dipole  transitions  induced  by  the  crystal  field  is  equal  to 


fjj'  =  X jj’a2-‘~u'  “4.52  , - ,  (32) 

Ejj'tJJ' 

where  a2  is  an  intensity  parameter  (in  cm)  which  can  be  determined  from  eq.  (32)  and  U jj ■  is  the 
matrix  element  of  the  tensor  operator  between  the  state  J  and  J’.  For  the  theoretical  calculation  of 
the  probability  of  the  absorption  transitions,  matrix  elements  of  the  tensor  operator  are  used: 

a2[u22(jf)+U24(jf)]=a2Ujr 

where  U2(JJ')  and  U4(JJ')  are  the  matrix  elements  of  the  tensor  operator  for  electric  dipole 
transitions  induced  by  a  CF  potential  having  an  odd  order.  The  intensity  parameter  may  be 
calculated  using  the  following  formula: 


a2 


4.52- 


EJJ'rJJ'U JJ' 


(34) 


Using  a  known  value  for  a2  it  is  possible  to  calculate  the  oscillator  strength  and  decay  time: 
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s 


fjj>  -Ejj,a2U  JJ,,  Tjj'  -4.52 


JJ 


Ejrfjr 


(35) 


3.8.  Spectroscopic  diagrams  of  Co2+  in  tetrahedral  crystals  of  the  type  AnBVI 

3.8.1.  Diagrams  of  the  energy  states  of  Co2+  ions  in  tetrahedral  semiconductors  of  the  type 
AnBVI. 

To  qualitatively  describe  the  effect  of  the  CF,  electrostatic  and  spin-orbital  interactions  in  A11  BV1 
crystals,  we  have  presented  diagrams  of  the  energy  states  of  Co2+  in  terms  of  units  of  the 
interaction  parameters. 

Figure  9  shows  the  energy  levels  of  Co"  in  the  relative  units  E/Dq  and  E/F22  versus  the  relative 
parameters  F22/Dq  and  Dq/F22  when  ^  =  0,  and  the  Slater  parameters  F44  =  F22/I2.5.  The  left 
hand  side  of  the  diagram  corresponds  to  the  case  of  a  strong  CF  with  the  electrostatic  interaction 
parameter,  F22/Dq,  changing  from  0  to  1 .  At  F22/Dq  =  0,  all  energy  levels  combine  into  the  four 
CF  configurations:  t  ,  t  e  ,  t  e  ,  and  e  having  the  relative  energies  E/Dq  =  0,  10,  20,  and  30. 


Figure  9.  Energy  level  diagram  Co2+,  showing  E/Dq  and  EIF22  versus  the  relative 
parameters  F22/Dq  and  Dq/F22when  ^3d  =  0. 
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The  right  hand  side  of  the  diagram  corresponds  to  the  case  of  a  weak  CF,  where  the  electrostatic 
interaction  parameter,  Dq/  F22,  changes  from  0  to  1 .  When  Dq/  F22  =  0,  the  energy  levels  are 

2d- 

equivalent  to  those  of  a  free  Co”  ion. 

In  AnBVI  crystals,  the  parameter  Dq/  F22  changes  from  0.25  to  0.5.  The  energy  states  of  Co2+  ions 
in  terms  of  the  relative  parameter  Dq!F2i  are  shown  in  Figure  10  and  1 1  over  the  range  relevant 
for  the  three  crystal  types:  ZnS,  ZnSe,  and  ZnTe.  We  will  analyze  these  in  the  sections  which 
follow  starting  from  the  baseline  case  of  Co'  in  CdS  for  illustrative  purposes. 


Figure  10.  The  energy  levels  of  Co2+  in  the  relative  units  E/F22  versus  the  relative  parameter 

Dq/F22  =  0-0.5  at  ^3d  =  0. 


0,30  0,32  0,34  0,36  0,38  0,40  0,42  0,44  0,46  0,48  0,50 

ZnSe  ZnTe  Dq/F22 

Figure  11.  The  energy  levels  (with  labels)  of  Co2+  in  the  relative  units  E/F22  versus  the 
relative  parameter  DqtFn  =  0.3-0.5  at  £3 d  =  0. 
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3.9.  Spectra  of  Co2+  in  CdS 

In  order  to  obtain  a  simulation  of  the  spectra  of  Co  ,  it  is  necessary  to  calculate  the  intensity 
parameter,  ch,  for  the  4A2-  4Fi  luminescence  transition  having  a  degeneration  factor  2J+ 1  =  12, 
a  decay  time  of  1 1 15  usee,  an  energy  E  =  3000  cm"1,  and  an  internal  polarization  factor  of  s  = 
1/8.  We  derived  the  following  value  for  the  intensity  parameter:  ci2  =  2054682  cm”.  Using  this 
value  of  «2,  it  is  possible  to  calculate  the  cross  section  for  any  transition. 

3.9.1.  Absorption  spectra  ofCo2+  in  CdS 

The  Absorption  spectra  of  both  the  a-  and  n—  transitions  of  Co  in  CdS  in  the  2000  -  14000  cm" 
1  and  2000  -  6000  cm"1  energy  region  are  shown  in  Figures  12  and  13,  respectively.  The  major 
transitions  are  labeled  in  both  diagrams. 


Figure  12.  The  Absorption  spectra  of  both  the  a-  and  n-  transitions  of  Co2+  in  CdS  in  the 

2000  -  14000  cm'1  energy  region. 
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Figure  13.  The  absorption  spectra  of  both  the  o-  and  n-  transitions  of  Co2+  in  CdS  in  the 

2000  -  6000  cm'1  energy  region. 


3.9.2.  Luminescence  spectra  of  Co2+  in  CdS 

The  Luminescence  spectra  upon  exciting  both  n-  and  a-  transitions  of  Co  in  CdS  are  presented 
in  Figures  14-19  over  various  spectral  regions.  Figures  14  and  15  show  luminescence  spectra  of 
excitations  on  n-  and  o-  4A2-4Fi  transitions  of  over  the  0-6000  cm"1  spectral  range.  Figures  16 
and  17  show  luminescence  spectra  of  excitations  on  n-  and  a-  4A2-(4P)4Fi  transitions  of  over  the 
8000-14000  cm'1  spectral  range.  Figures  18  and  19  show  luminescence  spectra  of  excitations  on 
7i-  and  a-  4A2-(2D)2F2  transitions  of  over  the  0-18000  cm'1  spectral  range. 
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6000 


6000 


Figure  15.  The  luminescence  spectra  in  the  0-6000  cm'1  region  of  Co2+  in  CdS  upon 
excitation  of  the  o-  transition  4A2-(4F)4Fi. 
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13000 


Figure  16.  Luminescence  spectra  in  the  8000-14000  cm'1  region  of  Co2+  in  CdS  upon 

excitation  of  the  ^-transition  4A2-(4P)4Fi. 


13000 


Figure  17.  Luminescence  spectra  in  the  8000-14000  cm'1  region  of  Co2+  in  CdS  upon 

excitation  of  the  o-transition  4A2-(4P)4Fj. 
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1BOOO 


Figure  18.  Luminescence  spectra  in  the  0-18000  cm'1  region  of  Co2+  in  CdS  upon  excitation 

of  the  n-  transition  4A2-(2D)2F2. 

18000 


Figure  19.  Luminescence  spectra  in  the  0-18000  cm'1  region  of  Co2+  in  CdS  upon  excitation 

of  the  a-  transition  4A2-(2D)2F2. 
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3.10.  Absorption  spectra  from  excited  states  of  the  (4F)4F2  multiplet 

2_j_ 

The  absorption  spectra  of  Co  in  CdS  for  both  the  a-  and  n-  transitions  from  the  excited  states 
of  the  (4F)4F2  multiplet  are  presented  in  Fig.  20. 


3000 


Figure  20.  Absorption  spectra  of  Co2+  in  CdS  for  both  the  a-  and  n-  transitions  from  the 

excited  states  of  the  (4F)4F2  multiplet. 


3.11  Selection  rule  for  the  transitions  (4F)4A2-  (4F)4F2  and  (4F)4A2-  (4F)4Fi 

The  representation  of  the  electric  dipole  operator  in  crystals  with  Td  CF  symmetry  is  F2. 
Transitions  from  the  A2  initial  state  are  only  possible  to  the  states  having  the  representation 
A2xF2  =  Fi.  This  is  why  the  transition  4A2-  (4F)4Fi  is  allowed  and  the  transition  4A2-  (4F)4F2  is 
forbidden.  Yet,  this  rule  holds  only  for  the  case  of  a  strong  CF.  For  the  kind  of  CF  commonly 
found  in  semiconductors,  this  rule  is  violated  because  of  the  spin  -  orbit  interaction.  A  very 
interesting  task  is  to  investigate  the  dependence  of  the  ratio  of  cross  sections  for  the  two 
transitions  G  =  <2(4A2-  (4F)4F2)/<2(  4A2-  (4F)4Fi)  on  the  ratio  of  the  two  main  intensity  parameters 
B4  IBi.  With  an  increase  in  the  fourth  order  intensity  parameter  relative  to  the  second  order 
intensity  parameter,  B41/B11,  from  0  to  0.8,  the  ratio  G  increases  from  0.535  to  1.426  (see  Table 
6). 
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Table  6.  Dependence  of  the  ratio  of  the  cross  section  for  the  two  transitions 
G=Q(4 A2-  (4F)4F2)/0(  4A2-  (4F)4Fi)  on  the  ratio  of  the  two  main  intensity  parameters  B4/B2. 


B 21-  B21 

■u 

II 

bs 

4- 

B41I B21 

G 

1 

0.795 

0.795 

1.426 

1 

0.696 

0.696 

1.29 

1 

0.596 

0.596 

1.14 

1 

0.496 

0.496 

0.998 

1 

0.396 

0.396 

0.86 

1 

0.298 

0.298 

0.727 

1 

0.198 

0.198 

0.625 

1 

0 

0 

0.535 

3.12  Configuration  coordinate  model  of  transition  dynamics 

To  understand  the  dynamics  of  the  electronic  transitions  and  the  influence  of  fluctuations  of  the 
atoms  in  the  crystal  lattice  on  these  dynamics,  a  system  of  configuration  coordinates  is  used  to 
define  the  ion  energy.  The  origin  of  the  configuration  coordinate,  R,  is  defined  to  be  where  the 
minimum  in  the  potential  energy  curve  for  the  ground  state  of  the  active  ion  occurs.  For  excited 
states,  there  is  a  shift,  Ro,  associated  with  the  minimum  in  the  potential  energy  curve  relative  to 
the  ground  state. 

Ug(r)  =  \kg(r/Ry  (36) 

Ue{r)  =  \ke((--^)2  +  U0.  (37) 

The  potential  functions  for  the  ground  ( Ug(r ))  and  excited  state  ( Ue(r ))  in  equations  (36)  and  (37) 
have  spherical  symmetry  and  form  one  dimensional  configuration  curves  for  the  ground  and 
exited  states  of  the  active  ions  in  semiconductors.  These  curves  determine  the  dynamics  of  the 
transition. 


In  addition,  crystal  fluctuations  will  modulate  the  potential  of  a  crystal  field  as  well  as  the  Dq 
parameter  whose  value  can  be  represented  as 

Dq{R ) 


Dq(R,r)  — 

where  k  is  the  coefficient  of  electron  phonon  coupling. 


(38) 


3.13  Spectra  of  Co2+  in  ZnSe 
3.13.1  Absorption  spectra  ofCo2+  in  ZnSe 

We  calculated  the  absorption  spectra  of  Co"  ions  using  the  following  parameters: 

Electrostatic  and  spin-orbital  interactions:  F22  =  965  cm'1,  F 44 =  76  cm'1,  B  =  F22-5F44  =  5  8  5  cm'1, 
C  =  35F44=  2660  cm'1,  ^d=  -535  cm'1; 
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Crystal  field  parameters:  Dq  =  361  cm'1,  A40  =  3^70Dq  =  9061  cm'1,  A44=  15Dq  =  5415  cm'1; 
A20  =  A21  =  A2.1  =  A22  =  A2.2  =  0,  A41  =  A44  =  A42  =  A4.2  =  A43  =  A4.3  =  0. 

Intensity  parameters  for  a  polarization:  B20  =  0,  B22  =  1 ,  B40  =  0,  B41  =  0,  B42  =  -2/V5  B22*0. 1 ,  B43 
=  0,  B44  =  0. 

Intensity  parameters  for  n  polarization:  B2o=  0,  B21  =  1,  B22=  0,  B41  =  1/V 1 0  B2i*0.1,  B42=  0;  B43 
=  -V7/Vl0*0.1  B2i,  B44=  0. 

Parameter  s  of  the  local  polarization  at  n  =  2.44:  s  =  0.347. 

Shown  in  Table  7  are  the  energy  levels  and  oscillator  strengths  of  the  absorption  spectra 
transitions  of  Co  in  ZnSe  (with  polarizations  k,  o,  and  total  polarization). 

Table  7.  Energy  levels  and  oscillator  strengths  of  the  absorption  spectra  transitions  of  Co2+ 
in  ZnSe  (with  polarizations  n ,  a,  and  total  polarization). 


JV2 

J' 

r 

E 

fg  (tt)*107 

fg(G)*  107 

fg(a)+2fg( n)  *107 

5 

(4f)4t2 

r6 

3518 

15 

0,7 

31 

7 

r8 

3593 

161 

1,4 

323 

11 

r8 

3670 

343 

1,5 

688 

15 

r7 

3732 

40 

1,8 

82 

17 

(4F)4T  1 

r6 

6030 

1214 

56,8 

2486 

19 

r8 

6088 

3436 

26,1 

6898 

23 

r8 

6418 

630 

44,9 

1304 

27 

r6 

6609 

277 

8,8 

562 

29 

(2G)2E 

r8 

11655 

24 

3,4 

52 

33 

(2G)2T  \ 

r8 

12192 

542 

27,1 

1111 

37 

r6 

12243 

309 

27,8 

645 

39 

(4P)4Ti 

r6 

13208 

20497 

31,2 

41025 

41 

r8 

13306 

33632 

174,6 

67438 

45 

r8 

13598 

22393 

73,8 

44860 

49 

r6 

13601 

11912 

182,6 

24007 

51 

(2G)2A, 

r6 

14074 

122 

1,3 

245 

53 

(2G)2T2 

r6 

14432 

68 

1,0 

137 

55 

r8 

14596 

47 

0,3 

94 

59 

(2H)2T2 

r8 

15967 

14 

0,2 

28 

63 

r6 

16009 

135 

0,1 

270 
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65 

(WTi 

r6 

16172 

19 

0,0 

38 

67 

r8 

16458 

467 

1,7 

936 

71 

(2H)2E 

r8 

17660 

63 

0,0 

125 

75 

(2P)2T2 

r6 

18443 

62 

14,2 

139 

77 

r8 

18631 

223 

10,7 

457 

81 

(WTi 

r6 

19712 

1 

1,7 

4 

83 

r8 

19797 

19 

1,0 

39 

87 

(2D)2T2 

r8 

20359 

10 

0,2 

20 

91 

r6 

20585 

6 

0,3 

12 

93 

(2D)2E 

r8 

22719 

10 

0,0 

19 

97 

(2f)2a2 

r6 

25683 

0 

0,0 

0 

99 

("F)2T2 

r6 

25963 

2 

2,5 

7 

101 

r8 

26129 

3 

1,5 

7 

105 

TfTtT 

r6 

27262 

6 

30,6 

42 

107 

r8 

27369 

36 

34,3 

107 

111 

(2D)2E 

r8 

39836 

6 

0,0 

12 

115 

(W  t2 

r8 

40477 

111 

0,4 

222 

119 

r6 

40517 

90 

0,0 

180 

— • - - - - . - . - - t - ; -  9_j_ 

Shown  in  Figure  21-23  are  the  absorption  spectra  in  the  3000-42000  cm"  region  of  Co“  in 
ZnSe.  The  major  absorptions  are  labeled  in  each  figure 


E  (sm  ') 

Figure  21.  Absorption  spectra  in  the  3000-42000  cm'1  range  of  Co2+  in  ZnSe. 
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E  (sm ') 

Figure  22.  Absorption  spectra  in  the  3000-24000  cm'1  range  of  Co2+  in  ZnSe. 


Figure  23.  Absorption  spectra  in  the  11000-15000  cm'1  range  of  Co2+  in  ZnSe. 
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3.14  Configuration  coordinate  model  of  the  transition  dynamics  of  Co2+  in  ZnSe 

2+ 

Shown  in  Figure  24  are  the  energy  levels  of  Co  in  ZnSe  versus  the  relative  configuration 
coordinate  r/R.  The  major  states  are  labeled  in  the  diagram. 


Figure  24.  Energy  level  diagram  of  Co2+  in  ZnSe  versus  the  relative  configuration 

coordinate  r/R . 

Parameters  associated  with  the  minima  in  the  configuration  curves  (CC)  for  the  ground  and 
excited  states  as  well  as  the  absorption  energies  from  the  (4F)4T2  ground  state  with  k  =  0.2  are 
given  in  Table  8.  These  parameters  were  calculated  using  the  formulas  given  in  equations  (36- 
38). 

Table  8.  Parameters  associated  with  the  minimum  in  the  configuration  curves  (CC)  as  well 
as  the  absorption  energies  from  the  (4F)4T2  ground  state  with  k  =  0.2. 


Parameter  of  CC  minimum 

Absorption  energies 

States 

r/R 

E(cm'1) 

(4F)4T2 

(2G)zT!rE 

(4P)4Tj 

(4F)4A2 

-0.341 

-25061 

-20145 

-16355 

-12989 

(4F)4T2 

-0.03 

-21034 

(2G)2T,;2E 

0.11 

-17875 

(4P)4Ti 

-0.338 

-13367 

The  dependence  of  configuration  curve  energies  vs  r/R  for  the  ground  level  4A2(r8)  and  the 
excited  level  4r2(r8)  are  given  by  the  following  equations: 

4A2(r8):  £a=1462x3+6821  x2+4143  x-24383;  (39) 

Xmin  =  -0. 341 ;  Ea>min  =  -25061;  Ea(x  =  - 0.03  )  =  -24501. 

4T2(r8):  £>=394x3+9336  x2+557.7  v -2 1026;  (40) 
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xm,n  =  -0.03;  Eejnin  =  - 21034;  Ee(x  =  -0.341  )  =  -20145. 


A  diagram  of  the  energy  levels  of  the  4A2(r8)  state  (the  lower  curve)  and  the  ^(Ts)  state  (the 
upper  curve)  of  Co2+  in  ZnSe  versus  the  relative  configuration  coordinate  x  =  r/R  are  shown  in 
Figure  25. 


Figure  25.  Diagram  of  the  energy  levels  of  the  4A2(r8)  state  (the  lower  curve)  and  the 
4T2(r8)  state  (the  upper  curve)  of  Co2+  in  ZnSe  versus  the  relative  configuration  coordinate 

x=r/R. 

Shown  in  Figure  26  is  a  diagram  of  the  energy  levels  4A2(r8)  (the  lower  curve)  and  4T2(r8)  (the 

2”F 

upper  curve)  of  Co  in  ZnSe  with  the  energies  and  configuration  coordinate  values  for  the 
absorption  and  emission  transitions. 


Figure  26.  Diagram  of  the  energy  levels  4A2(r8)  (the  lower  curve)  and  4T2(r8)  (the  upper 
curve)  of  Co2+  in  ZnSe  with  the  energies  and  configuration  coordinate  values  for  the 

absorption  and  emission  transitions. 
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Generally,  the  configuration  curve  may  be  represented  as  a  polynomial  line  as  shown  in  Fig  27 
for  ground  and  upper  states. 


(cm'1) 


Figure  27.  Diagram  of  the  energy  levels  4A2(r8)  (the  lower  curve)  and  4T2(r8)  (the  upper 
curve)  of  Co2+  in  ZnSe  along  with  the  photon  and  phonon  parameters  of  the  absorption  and 

emission  transitions. 


The  energies  of  the  ground  and  excited  states  in  the  parabolic  case  may  be  represented  by 
following  formula: 

Ea{x)  =  ka  (x-x0 a)2  and  Ee(x)  =  ke(x-x0ef,  (41) 

where  ka  and  ke  are  the  coefficients  of  the  harmonic  vibrations  (CHV)  associated  with  the 
absorption  and  emission  transitions  in  cm'1. 

The  eigenfunctions  of  the  absorption  and  emission  states  are  given  by  the  following 
eigenfunctions  for  harmonic  oscillations 

KO)  =  Nmexp[-[^]2  /2]Hm(x/Aa)  (42) 


where  Hm  and  Nm  are  Hermite  functions  and  N  is  the  normalization  factor.  The  vibration 
amplitude  Aa  at  the  point  of  the  energy  minimum  for  the  ground  state  is  equal  to: 

h  h 


A 2  - 
a  R/Mk^ 


(43) 


R2Ma)a  ’ 

the  vibration  amplitude  Ae  at  the  point  of  the  energy  minimum  for  the  excited  state  is  equal  to: 

*2  — - —  — - —  (44) 


A 2  — _ - _ 

e  RjMke 


R2McoP 


where  R  is  the  M  -  ligand  distance. 

The  phonon  frequencies  for  the  absorption  and  emission  transitions  are  equal  to 
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Vn  - 


2nR  ^  M 


—  V  = 
a  /t  >  ve 


2nR 


(45) 


The  energies  of  an  optical  phonon  are  defined  as  0.03  leV  =  250cm'1  for  longitudinal  phonons 
and  0.0263eV  =  212cm'1  for  transverse  phonons  in  [13]. 

The  Gaussian  form  of  the  absorption  band  around  the  central  frequency  vfi  is  described  by 


/(v)  = 


V27rcra 


exp[-(v-vj})/ 2ol], 


(46) 


where  the  band  widths  are 


2  _  GoR)2  ke  _  Go)2  fcf  _2  _  Gofi)2  fca  _  Go)2  ka 

a  27 ih  R3JkaM  27 ih  R^JkaM’  e  27 ih  R3JkeM  27 ih  RjkeM’  ^  ' 

2 

and  M  =  4 (Mse+orMzn)  ,  a  is  the  coefficient  of  a  chemical  bond  Zn-Se  (0.2-0. 4),  and  ka  and  ke  are 
the  parabolic  coefficients  of  the  vibration-electron  coupling  of  the  absorption  and  emission 
bands,  respectively. 

The  energies  Sa  and  Se  are  defined  as: 

Sa=Nahve ,  Se=Nehva,  (48) 

where  Na  is  the  average  number  of  quanta  (of  heat  energy  or  phonons)  with  energy  hve  emitted 
after  photon  absorption;  and  Ne  is  the  average  number  of  phonons  with  energy  hva  which  are 
absorbed  after  photon  emission. 


The  bandwidths  of  the  emission  and  absorption  transitions  are  equal  to 


^  =  Nef, 

v  P. 


<=Na, 


(49) 


2_|_ 

According  to  equations  (39)  and  (40)  for  Co'  ,  the  CHV  coefficients  of  the  emission  and 
absorption  transitions  are  equal  to:  ka  =  6824  cm'1  and  ke  =  9336  cm'1.  The  coordinates  of  the 
energy  minimums  of  the  absorption  and  emission  transitions  are  x(“  =  —0.341,  and  Xq  = 
—0.03.  Therefore,  starting  with  our  parameters,  the  energy  of  the  phonons  for  the  absorption  and 
emission  transitions  is  calculated  to  be  va=  43  cm'1  and  ve=  51cm'1. 


3.15  Spectral  properties  of  Ni2+  in  ZnSe 

3.15.1  Ni2+  energy  state  diagrams  in  tetrahedral  semiconductors  of  the  type  AnBVI. 

To  qualitatively  find  the  parameter  range  of  the  crystal  field  as  well  as  the  electrostatic  and  spin- 
orbital  interactions  in  AnBVI  crystals,  we  calculated  energy  state  diagrams  of  Ni2+  in  terms  of  the 
relative  units  of  the  interaction.  An  energy  level  diagram  for  Ni2+  in  ZnSe  is  shown  in  Figure  28. 

The  energy  level  diagram  of  Ni~  in  ZnSe  in  terms  of  the  relative  units  E/Dq  and  E/F22  versus 
the  relative  units  of  the  Eii/Dq  and  Dq/Fn  parameters  with  zero  spin-orbital  interaction  £3 d  =  0, 
and  F44  =  F22/I2.5  is  shown  in  Figure  29.  The  left  side  of  the  diagram  corresponds  to  a  strong 
crystal  field  where  electrostatic  interaction  parameters  F2i/Dq  changes  from  0  to  1 .  When  Fii/Efq 
=  0  all  of  the  energy  levels  converge  to  three  crystal  field  configurations  having  relative  energies 
E/Dq  =  0,  10,  and  20. 
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Figure  28.  Diagram  of  Ni2+:ZnSe  energy  levels  (in  cm'1). 
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Figure  29.  Ni2+:ZnSe  energy  level  diagram  of  E/Dq  and  E/F22  versus  the  relative  parameter 
values  F22/Dq  and  Dq/F22  with  no  spin-orbital  interaction  ^d=  0. 
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The  right  hand  side  of  the  diagram  shown  in  Figure  29  corresponds  to  a  weak  crystal  field  where 
the  electrostatic  interaction  parameters  Dq/F22  changes  from  0  to  1 .  At  DC//F22  =  0,  the  energy 
levels  are  equal  to  the  energy  levels  of  a  free  Ni2+  ion. 

Because  the  parameter  Dq/  F22  changes  from  0.25  to  0.5  in  crystals  of  the  type  AnBVI,  we  plotted 
this  region  separately,  see  Figure  30. 

E/F22 


Figure  30.  Energy  levels  of  Ni2+:ZnSe  in  terms  of  E/F 22  versus  the  relative  parameter 

Dq/F 22  =  0.25-0.5  at  i;3d  =  0. 

.9+ 

Shown  in  Table  9  are  the  strengths  of  the  absorption  transitions  of  Ni  in  ZnSe  for  both  the  n 
and  a  polarizations  as  well  as  the  total  polarization.  State  nominations  are  given  in  the  following 
form:  (2S+1L)(  2S+1Tji)  2S+1r.  A  quantum  number  of  the  fonn  2S+1L  for  the  state  of  a  free  Ni2+  ion 
is  given  in  first  parentheses  in  Table  9.  Quantum  numbers  of  the  form  “  Tji  for  a  state  of  Ni~ 
ion  in  a  crystal  field  with  consideration  of  the  spin  -  orbit  interaction  is  given  in  second 
parentheses  of  Table  9  and  the  resulting  states  2S+1r  are  given  without  parentheses. 


Table  9.  Strengths  (f*107)  of  the  absorption  transitions  ( J-J ')  of  Ni2+  in  ZnSe  for  both  n  and 
o  polarizations  as  well  as  the  total  polarization  (transitions  occur  from  the  base  state  J  = 

(3F)(3T1)3A1). 


N 

J' 

g 

£(cm‘) 

f-n  f  Bn  ) 

/-o  (Bn) 

/-(<r+2;r) 

1 

(JF)(jT1)iA1 

1 

0 

0 

0 

0 

2 

(JF)  (3Ti)jT! 

3 

309 

16387 

19131 

51905 

5 

(JF)  (3Ti)jE 

2 

870 

22288 

34966 

79542 

7 

(JF)  (jTi)jT2 

3 

1020 

81867 

60862 

224596 

10 

(JF)  (jT2)3 Tj 

3 

4555 

99376 

53359 

252112 

13 

C*F)  (3T2)jE 

2 

4593 

178213 

36583 

393009 

15 

(JF)  (3T2)3T2 

3 

4658 

222948 

24919 

470815 

18 

(JF)  (jT2)3A2 

1 

4882 

88633 

9449 

186715 
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19 

(JF)jA2 

3 

8593 

182096 

23538 

387730 

22 

(‘GyE 

2 

9482 

32908 

7358 

73173 

24 

~7g)% 

3 

9640 

298994 

36078 

634066 

27 

(JP)JT! 

2 

11806 

356359 

438772 

1151490 

29 

(JP)3Ti 

3 

11858 

688707 

601059 

1978473 

32 

(3P)3T! 

3 

12175 

968688 

492849 

2430224 

35 

(3P)JTi 

1 

12308 

316133 

134627 

766892 

36 

(Wt, 

3 

14603 

19210 

11303 

49724 

39 

('G)'A, 

1 

15730 

5914 

2790 

14618 

40 

(‘D)1^ 

3 

15930 

5700 

2071 

13471 

43 

TWe 

2 

18741 

2715 

3825 

9255 

45 

(*S)1A1 

1 

34011 

3152 

7647 

13951 

•  2+ 

The  energies  E  and  the  oscillator  strength/ of  the  absorption  spectra  of  Ni  in  ZnSe  for  the  total 
polarization  are  shown  in  Figure  31.  Also  shown  in  the  figure  are  the  transitions  involved. 
Shown  in  Figure  32  are  the  energies  E  and  oscillator  strength /of  the  absorption  spectra  of  Ni~ 
in  ZnSe  for  the  tc  and  a  polarizations. 

The  oscillator  strength /is  given  by: 

/=  hvfm\04exp[-(hv-Em2)/s2]  (50) 

where  5  =  350  cm fm  =  QmEm,  Qm  =  a2'Zn=l[U% (n, rri)  +  [/f  (n,m)]  ,  and  a2  =  3.56><10'5. 
Shown  in  Figure  33  is  the  oscillator  strength/ of  the  theoretical  spectra  for  the  n  polarization  (the 
upper  red  line)  and  the  experimental  spectra  (the  lower  black  line)  in  arbitrary  units  [14]  for  the 
absorption  of  Ni"  in  ZnSe  in  the  wavelength  region  0  -  2500  mn.  Shown  in  Figure  34  is  the 
oscillator  strength  /  of  the  theoretical  spectra  for  the  k  polarization  (the  upper  red  line)  and  the 
experimental  spectra  (the  lower  black  line)  in  arbitrary  units  [14]  for  the  absorption  of  Ni"  in 
ZnSe  in  the  wavelength  X  region  of  1000  -  2600  mn. 


■  a+2TC-polarisation 


E  (sm'1) 

Figure  31.  The  energies  E  and  oscillator  strength/for  the  absorption  spectra  of  Ni2+  in 

ZnSe  for  the  total  polarization. 
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E  (srrf1) 


Figure  32.  The  energies  E  and  oscillator  strength /of  the  absorption  spectra  of  Ni2+  in  ZnSe 

for  the  n  and  o  polarizations. 


Figure  33.  Oscillator  strength /of  the  theoretical  spectra  for  the  n  polarization  (the  upper 
red  line)  and  the  experimental  spectra  (the  lower  black  line)  in  arbitrary  units  [14]  for  the 
absorption  of  Ni2+  in  ZnSe  in  the  wavelength  region  0  -  2500  nm. 
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Figure  34.  Oscillator  strength /of  the  theoretical  spectra  for  the  n  polarization  (the  upper 
red  line)  and  the  experimental  spectra  (the  lower  black  line)  in  arbitrary  units  [14]  of 
absorption  of  Ni2+  in  ZnSe  in  the  wavelength  k  region  of  1000  -  2600  nm. 


3.16  Configuration  coordinate  model  for  the  transition  dynamics  of  Ni2+  in  ZnSe 

The  energy  levels  of  Ni  in  ZnSe  versus  the  relative  configuration  coordinate  r/R  is  presented  in 
Figures  35-36.  The  major  states  are  labeled  in  the  diagram. 


E  (cm1) 


r/R 

Figure  35.  Energy  level  diagram  of  Ni2+  in  ZnSe  versus  the  relative  configuration 

coordinate  r/R. 
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r/R 

Figure  36.  Energy  level  diagram  of  Ni2+  in  ZnSe  versus  the  relative  configuration 

coordinate  r/R . 


To  get  a  true  picture  of  the  influence  of  fluctuations  on  the  energy  characteristics  of  active  ions, 
it  is  necessary  to  include  in  the  Hamiltonian  an  operator  representing  the  spherically  symmetric 
interaction  energy  of  the  fluctuations  according  to  formulas  (36)  and  (37).  In  addition,  it  is 
necessary  to  include  a  modulation  multiplier  (38)  in  the  crystal  field  parameter  Dq. 

Minimum  values  for  the  energies  of  the  states  and  the  associated  configuration  coordinates  r/R, 
Table  10,  were  calculated  using  the  electron  phonon  parameters  kg=  2000  cm'1  and  k  =  0.2. 

Table  10.  Minimum  values  of  the  energies  of  the  states  and  the  associated  configuration 
coordinates  r/R  calculated  with  the  following  electron  phonon  parameters:  kg  =  2000  cm'1 

and  k  =  0.2 


N 

States 

g 

r/R 

E^cm"1) 

1 

(3F)(iT1)JA1 

1 

-0.198 

-11907 

2 

(JF) 

3 

-0.198 

-11907 

5 

(3F)  (jTi)jE 

2 

-0.198 

-11907 

7 

(JF)  (3T03T2 

3 

-0.198 

-11907 

10 

(3F)  (jT2)jT! 

3 

0.0725 

-6883 

13 

(jF)  (3T2)jE 

2 

0.0725 

-6883 

15 

(3F)  (3T2)jT2 

3 

0.0725 

-6883 

18 

(3F)  (jT2)3A2 

1 

0.0725 

-6883 

19 

(3F)3A2 

3 

0.198 

-2914 

22 

(1G)‘E 

2 

0.0437 

-1591 

24 

(1G)1T2 

3 

0.0437 

-1591 

27 

TpTt, 

2 

0.0338 

536 

29 

TpTt, 

3 

0.0338 

536 

32 

(3P)jT, 

3 

0.0338 

536 

35 

(3P)JT, 

1 

0.0338 

536 
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LUO 


•  •  •  2-t-  •  •  • 

Shown  in  Figure  37  are  the  energy  levels  of  Ni  in  ZnSe  for  various  kg  (quoted  in  the  figure  as 
“diag”)  and  k  values. 


Figure  37.  Energy  levels  of  Ni2+  in  ZnSe  for  various  kg  (quoted  in  the  Fig.  as  “diag”)  and  k 

values. 
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•  •  •  2H“  •  •  • 

Finally,  shown  in  Table  1 1  are  the  values  for  the  energy  levels  of  Ni“  in  ZnSe  for  the  minimas 
of  the  relative  configuration  coordinate  xo  of  two  states. 

Table  11.  Values  for  the  energy  levels  of  Ni2+  in  ZnSe  for  the  minimas  of  the  relative 
configuration  coordinate  jco  of  the  two  states  (3F)(3Ti)Ai,  and  (3F)  (3T2)Ti. 


N 

SO  levels 

g 

jco  —0.218 

jco  =0.186 

1 

(jF)(JT1)A1 

1 

-11703 

-10523 

10 

(:iF)CT2)T, 

3 

-5411 

-7129 

3.17  Spectral  properties  of  Fe2+  in  ZnSe 

2+. 

The  full  interaction  matrix  of  Fe  in  ZnSe  has  a  dimension  of  210x210.  Eigenvalues  of  this 
matrix  and  the  intensity  matrix  are  solved  with  the  following  parameters  values: 


Electrostatic  and  spin-orbital  interactions:  F2  =  965;  F4=  76;^3d  =  -480;  B  =  F2-5F4;  C  =  35F4  = 
2660;  C/B  =  4.547; 

Crystal  field  parameters:  Dq  =  300;  A4o=  3V70Dq  =  7530;  A44  =  15Dq  =  4500;  A2o  =  A21  =  A22 
A41  =  A42  =  A43  =  0; 

The  intensity  parameters:  B2o  =  0;  B4o=  0;  644=  0;  S2=  10;  S3=  100;  s  =  9n/(n2+2)2;  g  =  12;  x  = 
108/ 1 06;  n  =  2.44; 

Intensity  parameters  for  71  polarization:  B2i=  1;  B4i  =  -2/a/40B2i  =  -0.316;  B43=  V  7/VlOB2i  = 
0.837; 

Intensity  parameters  for  a  polarization:  B22  =  1 ;  B42  =  2N5B22  =  0.894; 

Shown  in  Table  12  are  the  absorption  energies  and  strengths  of  Fe  in  ZnSe  for  both  the  total 
and  (n  and  a)  polarizations.  The  upper  state  excited  along  with  its  degeneracy  is  shown  when 
known.  Also,  the  transitions  are  assumed  to  occur  from  the  base  state  (5D4)5Ai. 

Table  12.  Absorption  spectra  of  Fe2+  ZnSe  for  both  the  total  and  (n  and  o)  polarizations 
(transitions  occur  from  the  base  state  (5D4)5Ai) 


N 

States 

g 

£'(cm‘1) 

/(7T)*107 

f(a)  *107 

f(2n+a)  *107 

1 

5d4 

1 

Ai 

0 

0 

0 

0 

2 

3 

Ti 

29 

35 

82 

153 

5 

2 

E 

57 

46 

107 

198 

7 

3 

t2 

97 

122 

275 

519 

10 

sd3 

1 

A2 

139 

59 

133 

250 

11 

3 

Ti 

2770 

3131 

52 

6315 

14 

3 

t2 

3029 

2585 

150 

5319 

17 

5d2 

2 

E 

3051 

790 

159 

1740 

19 

3 

t2 

3372 

1964 

30 

3958 
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Shown  in  Figure  38  are  the  energy  E  and  the  oscillator  strength /for  the  absorption  spectra  of 
Fc2  in  ZnSe  for  the  tc  and  a  polarizations  as  well  as  the  total  polarization.  The  transitions 
involves  are  labeled  when  known.  In  Figure  39,  the  energy  E  and  oscillator  strength  /  of  the 
absorption  spectra  of  Fe2+  in  ZnSe  are  shown  for  only  the  total  polarization.  Known  transitions 
are  labeled  in  the  figure.  Shown  in  Figure  40  are  the  energy  E  and  oscillator  strength /of  the 
absorption  spectra  of  Fe2+  in  ZnSe  for  the  a  polarization.  In  Figure  41  are  the  energies  E  and 
oscillator  strength  /  of  the  absorption  spectra  of  Fe  in  ZnSe  for  the  tc  polarization.  Finally, 
Figure  42  shows  the  oscillator  strength  /  of  the  theoretical  spectra  for  the  a  polarization  (the 
upper,  red  line)  and  the  experimental  data  in  arbitrary  units  (the  lower,  black  line)  [15]  for  the 
absorption  spectra  in  the  2500  -  3500  nm  spectral  region  of  Fe2+  in  ZnSe. 


E  (cm  ) 

Figure  38.  The  energy  E  and  oscillator  strength/ of  the  absorption  spectra  of  Fe2+  in  ZnSe 
for  the  n  and  o  polarizations  as  well  as  the  total  polarization. 


A  27T+CT  (Pol) 


Figure  39.  Energy  E  and  oscillator  strength/of  the  absorption  spectra  of  Fe2+  in  ZnSe  for 

the  total  polarization. 
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Figure  40.  Energy  E  and  oscillator  strength/of  the  absorption  spectra  of  Fe2+  in  ZnSe  for 

the  a  polarization. 
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Figure  41.  Energy  E  and  oscillator  strength/of  the  absorption  spectra  of  Fe2+  in  ZnSe  for 

the  n  polarization. 
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Figure  42.  The  oscillator  strength /of  the  theoretical  spectra  for  the  a  polarization  (the 
upper,  red  line)  and  the  experimental  data  in  arbitrary  units  (the  lower,  black  line)  [15]  for 
the  absorption  spectra  in  the  2500  -  3500  nm  region  of  Fe2+  in  ZnSe. 


3.18  Configuration  coordinate  model  of  transition  dynamics  of  Fe2+  in  ZnSe 

2+ 

Shown  in  Figure  43  is  a  diagram  of  the  energy  levels  of  Fe“  in  ZnSe  versus  the  relative 
configuration  coordinate  r/R.  In  Figure  44  is  a  diagram  of  the  energy  levels:  5D4(r8),5D3(Fg)  and 
3H6(r8)  of  Fe2+  in  ZnSe  versus  the  relative  configuration  coordinate  x  =  r/R. 
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Figure  43.  Diagram  of  the  energy  levels  of  Fe2+  in  ZnSe  versus  the  relative  configuration 

coordinate  r/R. 


Figure  44.  Diagram  of  the  energy  levels:  5D4(r8)  (the  lower  curve)  ,5D3(r8)  and  3H6(r8)  (the 
upper  curve)  of  Fe2+  in  ZnSe  versus  the  relative  configuration  coordinate  x  =  r/R. 
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Shown  in  Table  13  are  equations  describing  the  energy  levels  of  Fe2+  in  ZnSe  as  a  function  of 
the  relative  configuration  coordinate  x  =  r/R  when  k  =  0.2.  In  addition,  the  coordinates  xq  and 
values  of  the  energy  minima  E(x o)  for  these  levels  are  also  shown  in  Table  13.  Finally,  shown  in 
Table  14  are  values  of  the  energy  levels  of  Fe"  in  ZnSe  for  the  minima  of  the  relative 

5  5  3 

configuration  coordinate  xo  of  three  states  D4,  D4,  and  H6. 


Table  13.  Equations  of  the  energy  levels  of  Fe2+  in  ZnSe  versus  the  relative  configuration 
coordinate  x  =  r/R  when  k  =  0.2.  Also  shown  are  the  coordinates  xo  and  values  of  the 

energy  minimums  E(x o)  for  these  levels. 


SO  levels 

Energies  (cm1) 

Xo 

F(x0) 

2d4 

735.2xJ+8430.r+2069x-36924 

-0.125 

-37052 

2d3 

-5204x4+3485x3+13155.r-2622x-33695 

-0.096 

-33322 

3h6 

-4666x5+4648x4+2927x3+4448x2+5 1 1 3x-25752 

-0.461 

-27143 

Table  14.  Values  of  the  energy  levels  of  Fe2+  in  ZnSe  for  the  minima  of  the  relative 
configuration  coordinate  xo  of  the  three  states  5D4, 5D4,  and  3H6. 


SO  levels 

Xo  =  -0.125 

xo  =  -0.096 

Xo  =  -0.461 

sd4 

-37052 

-37046 

-36158 

sd3 

-33171 

-33322 

-30267 

3h6 

-26326 

-26207 

-27143 

3.19  Overview  of  experimental  data  in  literature 

The  goal  of  this  section  is  to  compare  spectra  from  different  sources  in  order  to  obtain  values  of 
the  energy  levels  and  linewidths.  These  values  were  then  used  in  calculations  throughout  section 
3  and  the  resulting  spectra  are  given  in  the  main  body  of  the  report. 

We  have  considered  the  spectral  properties  of  Co2+,  Ni2+  and  Fe2+  in  ZnSe  from  a  number  of 
experimental  works  and  have  decided  to  digitize  and  plot  them  together  on  the  same  graph  to 
enable  a  detailed  analysis. 


3.19.1  Experimental  spectra  ofCo2+  in  ZnSe. 

2_|_ 

For  Co-  ions  in  ZnSe,  the  experimental  data  as  well  as  fits  of  this  data  are  shown  across  various 
wavelength  ranges  in  Figure  45. 

For  the  lowest  energy  range,  see  Figure  45,  we  can  see  one  peak  at  3420-3450  cm'1  (orange  and 
green  lines  -  data  from  ref.  [1]).  The  black  line  (ref.  [16])  has  a  main  peak  at  3300  cm  1  and  some 
smaller  ones  at  3400,  3500  and  3660  cm'1.  Similarly,  the  violet  line  (ref.  [17])  has  a  main  peak  at 
3280  cm'1  and  three  smaller  ones  at  3350,  3580  and  3700  cm'1. 


55 

Approved  for  public  release:  distribution  unlimited. 


Wavelength,  cm 

Figure  45.  Low  energy  experimental  absorption  spectra  of  Co2+  in  ZnSe:  orange  line,  (a)  - 
ref.  [1];  green  line,  (b)  -  ref.  [1];  black  line,  (c)  -  ref.  [16];  and  violet  line,  (d)  -  ref.  [17], 


In  the  next  region,  5000  -7500  cm'1,  see  Figure  46,  one  more  experimental  curve  appears,  i.e., 
the  cyan  line  which  is  data  from  ref.  [18].  However,  this  data  displays  a  flat  behavior  because  of 
the  lack  of  cooling.  This  quite  a  bit  different  from  the  other  work  where  there  was  cooling,  i.e., 
the  temperature  was  77°K  in  the  work  of  ref.  [1],  and  the  temperature  was  10°K  in  the  work  of 
ref.  [16]  and  [17]. 


Wave  number,  cm’1 

Figure  46.  Infrared  absorption  spectra  in  the  5000-7500  cm'1  spectral  region  of  Co2+  in 
ZnSe:  orange  line,  (a)  -  ref.  [1];  green  line,  (b)  -  ref.  [1];  black  line,  (c)  -  ref.  [16];  violet 

line,  (d)  -  ref.  [17];  cyan  line,  (e)  -  ref.  [18], 
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We  start  with  fitting  the  experimental  curves  (a)  and  (b)  in  Figure  46.  In  Figure  47,  the 
experimental  curves  (a)  and  (b)  are  shown  together  with  curves  obtained  by  fitting  the  data  with 
the  sum  of  two  Gaussian  functions.  The  expressions  used  in  the  fits  are: 

Fit(a,  [1])  =  0.00538+0.5 123*Exp(-0.5*[(x-5935)/255]A2)+  (51) 

0.9941*Exp(-0.5*[(x-6724)/312]A2) 

Fit(b,  [1])  =  1 .2278+0.9433 *Exp(-0. 5 * [(x-5854)/267]A2)+  (52) 

2.1 1 1  l*Exp(-0.5*[(x-6640)/498]A2) 
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Figure  47.  Infrared  absorption  spectra  in  the  5000-7500  cm'1  spectral  region  of  Co2+  in 
ZnSe:  orange  line,  (a)  -  ref.  [1];  green  line,  (b)  -  ref.  [1];  together  with  curve  fits  (black 

dashes). 


Now  we  separately  present  fit  results  for  the  data  in  curves  (c)  and  (d),  see  Figure  48  and  49, 
respectively.  In  Figures  48  and  49,  we  show  the  fit  of  the  experimental  data  in  reference  [16] 
(Figure  46,  curve  (c)).  We  also  demonstrate  the  similarity  of  the  spectra  from  reference  [16], 
curve  (c)  and  reference  [17],  curve  (d)  in  Figure  49.  The  expressions  used  in  the  fits  are: 

Fit(c,  [16])  =  0.0532+1. 8245  *Exp(-0. 5  *[(x-5574)/42]A2)+0.8648*Exp(-0. 5  *[(x- 

5744)/46]  A2)+ 1 .6035  *Exp(-0.5  *  [(x-5963)/ 1 93]  A2)  +  (53) 

.5057*Exp(-0.5*[(x-6347)/25]A2)  +2.0427*Exp(-0.5*[(x-6642)/330]A2) 

Fit(d,  [17])  =  0.409+2.077*Exp(-0.5*[(x-5588)/13]2)+0.54*Exp(-0.5*[(x-5635)/68]2)+ 

0.99*Exp(-0.5*[(x-  5792)/14]2)  +0.9124*Exp(-0.5*[(x-5966)/174]A2)+  (54) 

0.50*Exp(-0.5*[(x-6364)/20]A2)  +1.46*Exp(-0.5*[(x-6654)/323]A2) 
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Summarizing,  in  this  infrared  interval,  we  observe  two  peaks  (at  about  5580  and  6650  cm'1)  with 
greater  intensity  plus  3-4  peaks  of  lesser  intensity. 
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Figure  48.  Infrared  absorption  spectra  in  the  5000-7500  cm'1  spectral  region  of  Co2+  in 
ZnSe.  Experimental  data  is  taken  from  ref.  [16]  and  is  represented  by  the  black  line  and 

the  curve  fit  is  represented  by  the  red  line. 


Wave  number,  cm1 


Figure  49.  Infrared  absorption  spectra  in  the  5000-7500  cm'1  spectral  range  of  Co2+  in 
ZnSe.  The  experimental  data  from  ref.  [17]  is  represented  by  the  blue  line  with  its  fit 
represented  by  the  green  line.  The  experimental  data  from  ref.  [16]  (curve  c)  is  also  fit  and 

is  represented  by  the  red  line. 


There  is  also  some  evidence  of  a  peak  found  only  in  paper  [17],  see  Figure  50,  where  the 
Gaussian  fit  is  defined  to  be: 
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Fit([  1 7])  =  0.005+0.01  l*Exp(-0.5*[(x-11550)/40]2). 


(55) 


11200  11300  11400  11500  11600  11700  11800 

Wave  number,  cm'1 

Figure  50.  Spectra  of  Co2+  in  ZnSe  from  11200  to  11800  cm'1  from  ref.  [17]  (blue  line)  and 

the  Gaussian  fit  to  the  peak  (green  line). 


Further,  there  are  two  more  peaks  that  are  only  found  in  refs.  [16]  and  [17],  see  Figure  51.  The 
fits  to  these  peaks  are  as  follows: 


Fit(c,  [16])  =  0.3246+0. 0584*Exp(-0. 5 *[(x-12151)/27]2)  (56) 

+0.0908*Exp(-0.5*[(x-12361)/134]2) 

Fit(d,  [17])  =  0.04+0.0337*Exp(-0.5*[(x-12145)/31]2)  (57) 

+0.0908*Exp(-0.5*[(x-12386)/207]2) 
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Figure  51.  Spectra  of  Co2+  in  ZnSe  from  1180  to  12800  cm'1:  (c),  black  line  -  experimental 
data  from  ref.  [16]  and  its  fit  (red  line);  (d)  -  experimental  data  from  ref.  [17]  (blue  line) 

and  its  fit  (green  line). 
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Shown  in  Figure  52  are  various  experimental  absorption  spectra  in  the  visible  region  from  the 
scientific  literature.  In  order  to  visualize  all  of  the  spectra  on  one  graph,  part  of  the  spectrum 
from  ref.  [  1 7]  was  divided  by  a  factor  of  ten.  The  spectrum  in  ref.  [  1 7]  demonstrates  the  richest 
behavior,  and  thus  in  the  next  two  figures,  we  present  a  mathematical  fit  to  it. 


Wave  number,  cm' 

Figure  52.  Absorption  spectra  of  Co2+  in  ZnSe  in  visible  region:  orange  line,  (a)  -  ref.  [1]; 
green  line,  (b)  -  ref.  [1];  black  line,  (c)  -  ref.  [16];  violet  line,  (d)  -  ref.  [17],  blue  line  (d’)  - 
data  from  ref.  [17]  multiplied  by  a  factor  of  0.1;  cyan  line,  (e)  -  ref.  [18], 


For  the  spectrum  shown  in  Figure  53  between  13000  and  14600  cm’1,  five  main  peaks  are 
identified  and  the  fit  to  these  peaks  is  described  as: 


Fit(  [  1 7])  =  0.0275+0. 287  l*Exp(-0. 5  *[(x-13137)/12]2)+ 

0. 1663Exp(-0.5*[(x-  13227)/38]2)  +0.0413*Exp(-0.5*[(x-13371)/23]A2)  +  (58) 

0.03 16*Exp(-0.5*[(x-13778)/5 1]A2)  +0.0640*Exp(-0.5*[(x-14306)/14]A2), 


For  the  spectrum  shown  in  Figure  54  between  14400  and  15400  cm’1,  two  peaks  have  been 
identified  and  their  fit  is  described  as: 


Fit(  [  1 7])  =  0.02+0.4722*Exp(-0.5*[(x-14382)/44]2)+  0.4411* 
Exp(-0.5*[(x-  14538)/67]2). 


(59) 
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Figure  53.  Absorption  spectra  of  Co2+  in  ZnSe  in  the  visible  region  from  13000  to  14600  cm 
experimental  data  from  ref.  [17]  (blue  line)  and  its  fit  by  the  sum  of  five  Gaussian 

functions  (green  line). 
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Figure  54.  Absorption  spectra  of  Co2+  in  ZnSe  in  visible  region  from  14400  to  15400  cm'1: 
experimental  data  from  ref.  [17]  (violet  line)  and  its  fit  by  the  sum  of  two  Gaussian  function 

(green  line). 


Absorption  spectra  are  given  in  references  [16]  and  [17]  for  the  next  frequency  range  from  15000 
-  18000  cm'1.  For  the  experimental  data  from  reference  [16]  shown  in  Figure  55,  the 
mathematical  fit  is  described  by: 
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Fit([16])  =  0.2702+0.0742Exp(-0.5*[(x-16076)/l  1 1]2)+  0.031 1*  (60) 

Exp(-0.5*[(x-  16580)/22 1  ]2)  +  0.021  lExp(-0.5*[(x- 17123)/1 1 1]2); 

and  the  mathematical  fit  to  similar  data  from  reference  [17],  shown  in  Figure  56,  is  given  by: 

Fit([  1 7])  =  0. 0083+0. 0535Exp(-0. 5  *  [(x-16067)/l  18]2)+  0.0247*  (61) 

Exp(-0.5*[(x-  16602)/284]2)  +  0.0122Exp(-0.5*[(x-  17096)/103]2). 
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Figure  55.  Absorption  spectra  of  Co2+  in  ZnSe  from  15500  to  18000  cm'1:  experimental  data 
from  ref.  [16]  (black  line)  and  its  fit  by  the  sum  of  three  Gaussian  functions  (red  line). 
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Figure  56.  Absorption  spectra  of  Co2+  in  ZnSe  from  15500  to  18000  cm'1:  experimental  data 
from  ref.  [17]  (violet  line)  and  its  fit  by  the  sum  of  three  Gaussian  functions  (green  line). 
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The  experimental  spectra  of  the  high  energy  region  from  18750  to  20750  cm'1  is  given  in 
references  [1,  16,  17]  and  is  shown  in  Figure  57. 


Figure  57.  High-energy  absorption  spectra  of  Co2+  in  ZnSe  from  18750  to  20750  cm'1: 
orange  line,  (a)  -  ref.  [1];  black  line,  (c)  -  ref.  [16];  violet  line,  (d)  -  ref.  [17]. 


We  have  separately  fit  the  three  experimental  lines  presented  in  Figure  57  and  have  obtained  the 
following  expressions  for  the  fits: 

For  line  (a),  ref.  [1]  (see  Figure  58): 


Fit(a,  [1])  =  0.8414+0.065Exp(-0.5*[(x-19099)/104]2)+  (62) 

0.1418Exp(-0.5*[(x-  19660)/70]2)  +  0.1514*Exp(-0.5*[(x- 
19986)/362]2)  +  0.5602Exp(-0.5*[(x-  20629)/!  1 1]2); 


For  line  (c)  -  ref.  [16]  (see  Figure  59): 


Fit(c,  [16])  =  0.2843+0.089Exp(-0.5*[(x-19124)/85]2)+  0.4457*  (63) 

Exp(-0.5*[(x-  19634)/38]2)  +  0.235  lExp(-0.5* 

[(x-  19933)/286]2)  +  0.8643Exp(-0.5*[(x-  20572)/82]2); 


For  line  (d)  -  ref.  [17]  (see  Figure  60): 


Fit(d,  [17])=0. 0678+0. 3308Exp(-0. 5 *[(x-19041)/2]2)+  0.9657*  (64) 

Exp(-0.5*[(x- 1961 1)/7]2)  +  0.1465*Exp(-0.5* 

[(x-  19700)/! 05]2)  +1.0048Exp(-0.5*[(x-  20528)/9]2). 
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Due  to  the  sharp  peaks  in  the  last  spectrum  (line  (d)),  it  is  probably  more  accurate  to  use  the 
peaks  as  specified  in  Figure  57:  19040,  19058, 19614,  19682,  19928,  20527,  and  20635  rather 
than  the  maximums  of  the  fit. 
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Figure  58.  Absorption  spectra  of  Co2+  in  ZnSe  from  18500  to  21000  cm'1:  orange  line  - 
experimental  data  taken  from  ref.  [1];  green  line  -  the  mathematical  fit. 


Wave  number,  cm 

Figure  59.  Absorption  spectra  of  Co2+  in  ZnSe  from  18500  to  21000  cm'1:  black  line  - 
experimental  data  taken  from  ref.  [16];  red  line  -  the  mathematical  fit. 
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Figure  60.  Absorption  spectra  of  Co2+  in  ZnSe  from  18500  to  21000  cm"1:  violet  line  - 
experimental  data  taken  from  ref.  [17];  green  line  -  the  mathematical  fit. 


And  finally,  we  digitized  and  converted  the  experimental  data  on  luminescence  taken  from 
reference  [16]  to  cm"1,  see  Figure  61,  where  the  maximum  values  for  each  peak  are  shown.  (The 
red  plus  signs  signify  the  location  of  the  maximum  of  each  peak.). 


Figure  61.  Luminescence  spectra  of  Co2+  in  ZnSe  over  the  2500  to  2800  cm"1  spectral  range 
(experimental  data  were  taken  from  paper  [16]). 
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3.19.2  Experimental  spectra  ofNi2+  in  ZnSe. 

Unfortunately,  the  literature  does  not  contain  much  experimental  absorption  and  luminescence 
spectral  data  for  Ni“  in  ZnSe.  Of  the  three  ions  in  this  study,  the  smallest  amount  of 
experimental  work  has  been  perfonned  on  Ni2+:ZnSe. 

Let  us  consider  an  experimental  spectrum  from  the  classical  work  in  reference  [14]  on 
Ni2+:ZnSe,  see  Figure  62.  As  before,  we  perfonned  a  fit,  in  this  case,  by  the  summation  of  three 
Gaussian  curves,  see  Figure  63. 

Fit  =  7.97*Exp(-0.5*[(x-4871)/639]2)+16.98Exp(-0.5*[(x-  (65) 

8827)/523]2)+86.54Exp(-0.5*[(x-  12083)/767]2) 


Wavelength,  nm 

Figure  62.  Absorption  spectra  of  Ni2+  in  ZnSe  over  the  500  to  2500  cm'1  spectral  range, 

taken  from  reference  [14]. 


Figure  63.  Absorption  spectra  of  Ni2+  in  ZnSe  from  4000  to  14000  cm'1  taken  from 
reference  [14]  in  cm'1  (black  line)  and  its  mathematical  fit  (red  line). 
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3.19.3  Experimental  spectra  ofFe2+  in  ZnSe. 

2+ 

In  [19],  a  source  material  involved  in  the  creation  of  Fe  :ZnSe  was  99.9%  FeSe.  In  this 
reference,  data  are  presented  for  two  samples  of  Fe2+:ZnSe  having  Fe  concentrations  of  3 x  1 0 1 7 
and  4x10  cm'  .  The  concentrations  of  iron  were  calculated  by  using  absorption  data  at  3.15  pm 
and  14°K  along  with  the  corresponding  absorption  cross  section  of  a  =  1.5x10"  cm-  [2],  For  our 
present  analysis,  we  chose  to  digitize  experimental  data  for  the  sample  having  the  greater 
concentration  of  iron  since  it  demonstrated  a  somewhat  richer  behavior  relative  to  the  other 
sample.  Reference  [19]  also  contained  spectra  of  the  shorter  wavelength  region  which  indicated 
that  the  incorporation  of  Fe  in  ZnSe  does  not  produce  any  significant  shift  of  the  absorption  edge 
in  the  visible  region  relative  to  the  absorption  edge  of  pure  bulk  ZnSe  which  occurs  around 
440nm.  The  absorption  spectrum  of  the  sample  with  the  smaller  concentration  of  iron  showed  a 
small  “bump”  near  540  nm  at  low  temperature  of  an  allegedly  unknown  nature.  And,  finally,  in 
both  samples,  there  were  absorption  band  peaks  around  1.5  pm  in  the  temperature  region  from  5 
to  300°K.  The  authors  attributed  this  to  either  Fe  or  Ni,  deducing  that  the  samples  were 
contaminated  with  Ni  from  the  PL  data  which  showed  a  2.6  pm  (3850  cm'1)  band  that  the 
authors  assigned  to  the  internal  transition  T2(F)— >  Ti(F)  and  a  3  pm  (3340  cm'  )  band  that  was 
assigned  to  the  transition  3E2(D)— »2T2(D)  of  Ni2+  in  ZnSe. 

We  also  have  digitized  the  absorption  and  luminescence  spectra  presented  in  reference  [2]  (see 
Figure  64).  The  luminescence  spectral  region  in  this  work  is  much  shorter  than  that  presented  in 
[19]  with  the  peak  being  located  at  about  3.75  pm  (2660  cm'1).  One  of  the  advantages  of  this 


Wavelength  (pm) 

Figure  64.  Absorption  and  luminescence  spectra  of  Fe2+  in  ZnSe  taken  from  [2]. 

work  is  that  there  are  absorption  spectra  at  room  temperature  as  well  as  at  14°K  presented  in 
absolute  values  rather  than  arbitrary  ones,  which  allows  us  to  compare  the  area  under  the  spectral 
curves. 
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2-t-  •  1 

Shown  in  Figure  65  are  the  absorption  spectra  of  Fe  in  ZnSe  over  the  2600  to  4000  cm" 
spectral  range  taken  from  the  following  references  [1,2,  and  19].  The  locations  of  the  maxima  of 
the  peaks  for  the  data  from  reference  [1]  are  shown. 


Figure  65.  Absorption  spectra  of  Fe2+  in  ZnSe  over  the  2600  to  4000  cm'1  spectral  range:  (a) 

ref.  [2],  (b)  -  ref.  [1],  (c)  -  ref.  [19]. 


We  have  fit  the  absorption  spectra  from  reference  [2]  by  the  sum  of  Gaussian  functions,  and 
obtained  the  following  expression  for  the  mathematically  generated  spectral  curve: 

Fit(a,[2])  =  1 .3Exp(-0. 5  *  [(x-2728)/ 1 7]2)+0. 5  8Exp(-0.5  *  [(x-2796)/ 1 8]2)  +  (66) 

1.26Exp(-0.5*[(x-2900)/101]2)+0.70*Exp(-0.5*[(x- 
3145)/77]2)+1.24Exp(-0.5*[(x-3350)/202]2). 

The  fit  along  with  the  experimental  data  are  shown  in  Figure  66. 
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Figure  66.  Absorption  spectra  of  Fe2+  in  ZnSe  over  the  2600  to  4100  cm'1  spectral  range: 
experiment  from  ref.  [2]  -  black  line  and  its  fit  -  red  line. 


•  •  •  •  "2, 1 
Finally,  the  location  of  the  peaks  of  the  emission  spectra  of  Fe“  :ZnSe  are  shown  for  the  data 

from  references  [2,  19]  in  Figure  67. 


Figure  67.  Luminescence  spectra  of  Fe2+  in  ZnSe,  blue  line  -  ref.  [19],  black  line  -  ref.  [2], 
Spectroscopic  characterization  of  TM  doped  samples 
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4.0  SPECTROSCOPIC  CHARACTERIZATION  OF  TM  DOPED  SAMPLES 


4.1  Concentration  and  temperature  quenching  of  mid-IR  photoluminescence  in 
Fe:ZnSe  and  Fe:ZnS  crystals 

In  this  project,  we  also  focused  on  concentration  and  temperature  quenching  of  mid-IR 
luminescence  in  Fe:ZnS  and  Fe.'ZnSe  crystals.  Parameters  associated  with  these  processes  are 
very  important  for  the  application  of  these  materials  in  lasers. 

Iron  doped  ZnSe  and  ZnS  crystals  were  prepared  using  the  thermo-diffusion  method  [1,  20].  The 
thickness  of  the  metallic  layer  and  the  diffusion  procedure  were  optimized  to  provide  a  unifonn 
distribution  of  iron  in  the  samples  studied.  Iron  doped  ZnS  and  ZnS  samples  having  an  Fe 
concentration  of  0.1-112x10  cm'  were  prepared  for  spectroscopic  experiments.  The 
concentration  of  Fe  ions  in  the  samples  was  estimated  using  a  value  for  the  absorption  cross- 
section  of  Gab=  0.65xl018  cm'2  at  2.698  pm. 

An  actively  Q-switched  Er:Cr:YSGG  laser  operating  at  2.78  pm  with  a  80  ns  pulse  duration  was 
utilized  for  optical  excitation  of  the  studied  samples.  The  Fe:ZnSe  and  Fe:ZnS  samples  were 
placed  in  a  closed-cycle  refrigerator  system  (Janis  Research  Co.,  Inc.,  Model  CCS-450)  for 
temperature  dependent  spectroscopic  measurements  over  the  14-300°K  temperature  range.  The 
photoluminescence  spectra  were  measured  using  an  Acton  Research  ARC-300i  spectrometer. 
The  kinetics  of  luminescence  for  the  given  samples  were  measured  with  an  InSb  detector 
(Teledyne  Judson  J10D-M204-R250)  equipped  with  a  Perry  amplifier  which  provides  a  10  ns 
(10/90%)  rise  time  and  DC  to  50  MHz  response. 

Room  temperature  absorption  and  luminescence  spectra  on  the  5E  <->572  transition  of  Fe2+  ions  in 
ZnSe  and  ZnS  crystals  are  shown  in  Figure  68.  A  dip  in  the  photo-luminescence  curves  at  4300 
nm  is  caused  by  absorption  of  atmospheric  CCF  in  the  experimental  setup. 


Wavelength  (nm) 

Figure  68.  Room  temperature  absorption  (A)  and  emission  (B)  spectra  of  Fe:ZnSe  (solid 

line)  and  Fe:ZnS  (dashed  line)  crystals. 
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4.1.1  Fe:ZnSe  crystal. 

Luminescence  lifetime  studies  over  the  15-300°K  temperature  range  were  reported  in  [21]  and 

IQ  2 

[3].  In  [21],  Fe:ZnSe  single-crystal  samples  with  an  iron  concentration  from  3.4  to  6.5x10  cm' 
were  reported.  In  the  temperature  range  from  12  to  120°K,  the  lifetime  unexpectedly  increased 
initially  from  33  to  105  ps  and  then  decreased  as  the  temperature  was  further  increased  due  to 
temperature  induced  non-radiative  processes.  Similarly,  an  unusual  growth  of  the  luminescence 
lifetime  with  an  increase  in  temperature  from  14  to  120°K  was  also  reported  in  [3]  where  it  was 
attributed  to  the  thermal  population  of  iron  sub-levels  having  a  smaller  oscillator  strength.  In  the 
current  study,  we  also  observed  an  increase  in  the  luminescence  lifetime  with  an  increase  in 
temperature  from  14  to  120°K  in  an  Fe:ZnSe  sample  having  an  iron  concentration  of  CFe  = 

1  o  o 

5  x  1 0  cm'  .  The  luminescence  kinetics  of  this  sample  at  20  and  77°K  are  shown  in  the  Figure  69 
(curve  iii  and  iv).  As  one  can  see  from  the  Figure  69,  the  luminescence  decay  reveals  a  single 
exponential  behavior  with  luminescence  lifetimes  of  x  =  54  ps  at  T  =  20°K  and  x  =  68  ps  at  T  = 
77°K.  However,  in  the  crystal  having  the  smallest  iron  concentration  CFe  =  0.1x10  cm'  ,  the 
kinetics  of  the  luminescence  at  low  temperature  were  similar  in  the  14-1 10°K  temperature  range 
with  the  luminescence  lifetime  being  ~x  =  57±4  ps  (see  Figure  69,  curves  i  and  ii).  Therefore,  we 
consider  the  unusual  increase  in  the  fluorescence  lifetime  between  14  and  1 10°K  reported  in  the 
papers  [3,  21]  to  be  attributed  to  iron-iron  interactions  as  well  as  luminescence  re-absorption  at 
low  temperatures.  The  highly  doped  sample  (CFe  =  86x10  cm'  )  also  revealed  an  asymptotic 
decay  at  T  =  110°K  (curve  v)  which  is  longer  than  that  at  T  =  20°K  (curve  vi),  but  the  low 
temperature  kinetics  of  this  sample  were  essentially  non-exponentials.  One  of  the  reasons  behind 
the  non-exponential  behavior  could  be  related  to  multi-center  formation  in  highly  doped  crystals. 

Figure  70  shows  the  low  temperature  (T  =  14°K)  luminescence  spectra  of  samples  having  iron 
concentrations  of  CFe  =  5X1018  cm'3  and  Cpe=  86x1018  cm'3.  The  spectrum  of  the  sample  with  an 
iron  concentration  of  Cpe  =  5 x  1 0  cm'  is  similar  to  what  has  been  published  before.  The 
luminescence  spectrum  of  the  highly  doped  sample  (CFe  =  86x10  cm'  )  is  slightly  red-shifted 
and  demonstrates  profound  peaks  on  the  long  wavelength  slope  of  the  luminescence  spectra.  The 
nature  of  multi-center  formation  in  a  highly  doped  Fe:ZnSe  crystal  is  a  subject  for  separate 
studies. 
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Figure  69.  Luminescence  kinetics  of  Fe:ZnSe  samples  having  iron  concentrations  of:  C  = 
0.1  xlO18  cm'3  (i,ii);  C  =  5><1018  cm 3  (iii,iv);  and  C  =  86><1018  cm 3  (v,vi);  at  T  =  20°K 
(ii;iv;vi),  T  =  110°K  (iii,v);  and  T  =  120°K  (i). 


3.5  4.0  4.5 

Wavelength,  gm 

Figure  70.  Luminescence  spectra  of  Fe:ZnSe  samples  at  T  =  14°K  having  the  following  iron 

10-1  i  o  -1 

concentrations:  (i)  CFe=  112x10  cm"  and  (ii)  CFe=  5x10  cm'  ;  and  an  Fe:ZnS  sample  at 
T  =  14°K  having  an  iron  concentration  of:  (iii)  CFe=  5xl018  cm'3. 

1  o  o 

The  temperature  dependence  of  the  luminescence  lifetime  of  a  low  doped  (CFe  =  0.1x10  cm' ) 
Fe:ZnSe  sample  is  shown  in  the  Figure  71.  As  it  was  mentioned  earlier,  in  the  14-120°K 
temperature  range,  the  luminescence  lifetime  was  unchanged  initially  at  t  =  57±4  ps  and  then 
decreased  with  temperature  due  to  thermally  activated  non-radiative  decay.  We  utilized  two 
models  of  the  non-radiative  process  to  fit  the  experimental  data.  Initially,  we  used  a  model 
describing  thermally  activated  radiationless  processes  based  on  the  intersection  of  the 
configuration  coordinates  of  the  upper  and  low  levels.  These  processes  are  typical  for  non- 
radiative  relaxation  in  optical  centers  having  a  strong  electron-phonon  coupling. 
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Figure  71.  Luminescence  lifetime  versus  temperature  for  Fe:ZnSe  (+  -  CFe=  O.lxlO18  cm'3) 
and  for  Fe:ZnS  crystals  (x-CFe  =  3X1018  cm'3,  □-  CFe=  73xl018  cm'3);  curves  describing  the 

fits  are  described  in  the  text. 

This  approach  leads  to  a  temperature  dependence  of  the  upper  level  lifetime  given  by  the 
following  equations  (see  [22]): 


=  x~ld  +  Wnr  (67) 

Wnr  =  W0exp(-AEa/kT),  (68) 

where  rmd  is  the  radiative  lifetime,  Wnr  is  the  nonradiative  relaxation  rate,  AEa  is  the  energy  gap 
(activation  energy)  between  the  intersection  of  the  adiabatic  potentials  of  the  ground  and  excited 
states  and  the  minimum  of  the  adiabatic  potential  of  the  excited  state,  Wo  is  the  relaxation 
parameter;  and  kB  is  the  Boltzmann  constant,  Figure  72.  The  best  fit  was  obtained  using  the 
parameters  AEa  =  2500  cm'1  and  1/Wa  =  0. 1  ns  (see  Figure  69,  curve  i). 


Figure  72.  Configuration  coordinate  diagram  for  linear  coupling  [23] 
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Temperature  measurements  of  the  low  doped  sample  allow  us  to  estimate  the  radiation  lifetime 
of  the  5E<rJ>5T2  transition  to  be  Tracj  =  57  ±  4  ps.  Using  this  value  for  the  radiative  lifetime,  the 
emission  cross  sections  at  RT  were  calculated  from  the  FUchtbauer-Ladenburg  equation  [1]  with 
the  emission  cross  section  equaling  aem  =  1.2x10'  cm'  at  4.3  pm. 

Figure  73  shows  the  room  temperature  luminescence  kinetics  of  a  Fe:ZnSe  crystal  with  an  iron 
concentration  varying  from  0.1x10  to  112x10  ‘  cm'  .  Luminescence  kinetics  of  the  low  doped 
Fe:ZnSe  crystal  (curve  1)  was  single  exponential  at  RT  with  a  luminescence  lifetime  Trt  =  382 
ns.  This  decay  time  is  significantly  smaller  than  the  radiative  lifetime  (w)  of  the 
transition  and  is  detennined  by  thermally  induced  non-radiative  processes.  The 
photoluminescence  kinetics  of  the  highly  doped  crystal  reveals  a  non-exponential  decay  (curve 
4).  However,  asymptotic  decay  curves  were  exponential  with  a  small  decrease  of  the  lifetime  to 
309  ns  for  a  sample  having  an  iron  concentration  of  CFe  =  112x10  cm"  in  comparison  with  382 
ns  for  a  sample  having  an  iron  concentration  of  Cpe  =0.1  x  10  cm" .  Decreases  in  the  asymptotic 
decay  rate  of  highly  doped  crystals  relative  to  lower  doped  crystals  as  well  as  the  non¬ 
exponential  decay  feature  are  characteristics  of  the  energy  migration  process.  The  decay  rate 
(Wq)  induced  by  concentration  quenching  can  be  calculated  from  the  difference  of  the  asymptotic 
decay  rates  of  the  low  and  highly  doped  samples.  But,  even  for  the  highly  doped  sample,  Wq  = 
0.6  ps’1  is  still  smaller  than  that  associated  with  temperature  quenching  of  the  luminescence. 
Simulation  of  the  photoluminescence  decay  of  highly  doped  samples  requires  a  numerical 
simulation  of  the  energy  migration  processes  to  include  those  associated  with  the  multi-center 
complexes  of  the  iron  in  the  highly  doped  crystal  (a  good  approximation  can  be  obtained  through 
the  double  exponential  decay  fit  of  the  photoluminescence  decay  curve  with  time  constants  of  r/ 
=  309  ns  and  Z2  =  138  ns).  However,  it  is  more  convenient  for  laser  applications  such  as  passive 
Q-switching,  mode-locking  and  gain-switching  to  estimate  the  average  relaxation  time  <xav> 
asj  l (t) / Imaxdt . 


Figure  73.  A)  Luminescence  kinetics  of  Fe:ZnSe  samples  with  different  iron 
concentrations:  (1)-  CFe  =  0.1  xlO18  cm'3;  (2)-CFe=  14xl018  cm'3;  (3)-CFe=  45xl018  cm'3;  and 

i  o  1 

(4)-  CFe=  112  xlO  cm'  ,  B)  same  luminescence  kinetics  with  an  expanded  time  scale. 


Figure  74  shows  the  dependence  of  the  average  relaxation  rate  on  the  iron  concentration.  As 
shown  in  Figure  74,  the  decay  rate  decreases  from  382  ns  for  a  lower  concentration  Fe:ZnSe 
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18  3 

sample  to  186  ns  for  a  Fe:ZnSe  sample  having  an  iron  concentration  of  CFe  =  1 12x10  cm'  .  The 
concentration  dependence  of  the  average  decay  time  can  be  fit  with  the  empirical  equation; 


Tr 

l+wFe/w0 


(69) 


2_j_ 

where,  NFe  is  the  concentration  of  Fe  ions,  and  No  is  the  critical  concentration  parameter  where 
the  decay  time  is  reduced  by  a  factor  of  2.  A  least  squares  fit  of  the  experimental  lifetime  data  is 

1  o 

also  shown  in  Figure  72  (dashed  curve).  The  best-fit  value  for  N0  was  determined  to  be  106x10 
cm'  .  Since  these  decay  times  are  longer  than  the  typical  pulse  duration  for  Q-switched  lasers 
(-100  ns),  lasers  based  on  Fe:ZnSe  crystals  are  able  to  lase  at  room  temperature  under  short- 
pulse  gain-switched  excitation. 


We  observed  a  decrease  in  the  average  relaxation  time  for  highly  doped  samples  due  to 
concentration  quenching  at  low  temperature.  The  average  relaxation  time  of  the  lower  doped 
sample  (Cpe  =  O.lxlO18  cm'3)  was  equal  to  the  decay  time  obtained  from  a  single  exponential  fit 
(57  ps)  at  T  =  20°K.  The  sample  having  an  iron  concentration  of  CFe  =  5X 10  cm'  had  a  slightly 
faster  relaxation  time  of  50  ps.  A  further  increase  in  the  concentration  results  in  a  shortening  of 

1  o  o 

the  relaxation  time  to  7  ps  when  the  iron  concentration  reaches  86x10  cm"  . 


18  3 

Fe  concentration  (xlO  ,  cm'  ) 

Figure  74.  Average  relaxation  time  of  the  5T2  — »SE  transition  versus  iron  concentration  in 
Fe:ZnSe  samples  measured  at  room  temperature  (dashed  line  shows  least  square  fit  using 

eq.  70). 


4.1.2  Fe:ZnS  crystal 

2_|_ 

The  Zinc  Selenide  crystal  is  one  of  the  most  frequently  used  active  element  hosts  for  Cr  and 
Fe  based  mid-IR  lasers.  However,  the  iron  doped  Zinc  Sulfide  crystal  is  also  very  attractive  for 
mid-IR  laser  applications.  First  of  all,  due  to  the  larger  crystal  field  splitting  of  the  5D  ground 
state  in  the  Fe:ZnS  crystal,  the  absorption  and  emission  spectra  are  shifted  to  a  shorter 
wavelength  in  comparison  with  the  ZnSe  host  (see  Figures  68  and  70).  As  a  result,  the  room 
temperature  luminescence  spectra  of  Fe:ZnS  cover  the  3-5  pm  spectral  range.  Due  to  the  strong 
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absorption  bands  of  hydrocarbons  and  other  organic  molecules  in  this  spectral  range,  Fe:ZnS 
lasers  are  very  promising  for  sensing  applications.  In  this  work,  we  report  the  temperature 
dependence  of  the  mid-IR  photoluminescence  of  Fe:ZnS  samples  having  iron  concentrations 
which  vary  over  the  (3-73)x10  cm"  range.  Figure  75A  shows  the  luminescence  kinetics  of 
Fe:ZnS  crystals  measured  at  T  =  220°K  for  three  different  iron  concentrations.  The  luminescence 
kinetics  were  found  to  be  single  exponential  with  a  luminescence  lifetime  of  r=  182  ns  for  all  of 
the  studied  samples.  Due  to  limitations  in  the  duration  of  the  pump  source  (80  ns),  the  lifetime  of 
the  Fe:ZnS  samples  could  not  be  obtained  above  250°K.  Decreasing  the  sample  temperature  to 
14°K  was  found  to  result  in  a  monotonous  increasing  of  the  luminescence  lifetime  to  t=  6.2  and 
5.5  ps  for  iron  concentrations  of  CFe  =  3x10  cm"  and  CFe  =  73  x  10  cm"  ,  respectively  (see 
Figure  71).  These  results  are  close  to  what  was  reported  earlier  in  [23]  where  a  lifetime  of  5.5  ps 
and  a  quantum  yield  of  10%  were  measured  for  this  transitions  at  low  temperature.  Using  these 
results,  we  can  estimate  a  radiative  lifetime  of  Trad  =  62-55  ps  for  the  Fe:ZnS  crystal. 

The  temperature  dependence  of  the  luminescence  cannot  be  fit  by  a  thermally  activated 
radiationless  relaxation  model  based  on  a  single  configuration  coordinate.  A  more  accurate 
model  should  consider  the  additional  physical  processes  along  with  a  more  accurate  energy 

structure.  Since  it  is  known  that  the  Jahn-Teller  interaction  results  in  additional  splitting  of  the 
upper  and  ground  state  D  levels,  a  more  accurate  model  should  therefore  include  several 
radiationless  processes  having  different  relaxation  parameters.  The  solid  curve  in  the  Figure  71 
shows  a  fit  using  the  two  relaxation  processes  described  by  equation  (68)  with  the  following 
parameters:  AEa/  =818  cm'1;  1/Woi  =  1  ns;  and  AEca  =  126  cm"1  ;  and  I/W02  =  420  ns.  At  high 
temperature,  a  relaxation  model  should  also  include  multi -phonon  processes  [22].  A  more 
accurate  model  will  require  additional  experimental  study  of  the  luminescence  lifetime  at  high 
temperatures  under  short  pulse  excitation.  Extrapolation  of  the  luminescence  lifetime  to  room 
temperature  allow  us  to  estimate  the  lifetime  as  r=  50  ns. 


Figure  75.  A)  Luminescence  kinetics  of  Fe:ZnS  samples  having  different  iron 
concentrations:  (1)-  CFe=  3  xlO18  cm'3;  (2)-CFe=  15x1018  cm'3;  (3)-CFe=  73x1018  cm'3  at  T  = 

T  8 

220°K;  B)  Luminescence  kinetics  of  Fe:ZnS  having  an  iron  concentration  of  CFe=  3x10 

cm'3  at  different  temperatures  (4-6) 
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4.2  Cobalt  doped  samples 


The  cobalt  concentration  gradient  was  measured  using  an  Er-fiber  laser  operating  at  1.56  pm.  To 
measure  the  concentration  gradient  with  high  spatial  resolution,  a  1  mm  output  beam  from  an 
Er-fiber  laser  was  focused  by  lens  with  a  focal  length  of  F  =  5  cm  to  a  beam  with  a  waist 
diameter  of  d1/e2  —  1.23  ■  A  ■  F/D1/e2  —  97  pm  and  a  confocal  depth  of 
ZR  —  n  d\/e2ll  —  1.9  cm.  To  calculate  the  concentration  gradient,  we  measured  the  variation  of 
the  transmission  of  the  laser  beam  through  a  1  cm  long  crystal  as  shown  in  Figure  76.  The 
cobalt  concentration  gradient  of  a  CoiZnSe  sample  annealed  for  7  days  is  also  shown  in  the 
figure  and  is  found  to  be  negligible. 
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Figure  76.  Experimental  setup  and  results  for  the  measurement  of  the  cobalt  concentration 

gradient  in  a  Co:ZnSe  crystal. 

Most  of  previous  studies  of  cobalt  doped  ZnSe  and  ZnS  crystals  were  focused  on  the 
4Ao(F)— >4Ti(F)  transition  near  1.6  pm.  This  particular  transition  of  Co2+  has  attracted  significant 

T_i_ 

interest  due  to  possible  utilization  in  saturable  absorber  crystals  in  Er  lasers  operating  at  1.5 
pm.  Several  authors  have  reported  the  absorption  cross-section  of  this  transition  (see  Table  6  for 
details).  Based  on  reported  values  for  the  absorption  cross-sections  of  Co2+,  we  used  absorption 
cross-sections  of  7.5x10'  cm'-  at  1625  nm  and  7.8x10'  cm'"  at  1500  mn  for  the  ground  state 
absorption  4A2(F)— >4Ti(F)  transition  in  ZnSe  and  ZnS  crystals,  respectively. 

The  4A2(F)^4T2(P)  near  infrared  transition  has  an  absorption  cross-section  approximately  one 
order  of  magnitude  bigger  (e.g.  64x10'  cm  at  768  nm  for  Co:ZnSe).  It  should  be  noted  that 
because  the  absorption  cross-sections  for  the  entire  700-800  nm  band  are  higher  than  the 
absorption  cross-section  at  1625  nm,  this  band  could  be  used  for  optical  excitation.  However, 
because  of  the  high  absorption  coefficient  associated  with  this  band,  only  a  side-pumping 
configuration  is  feasible. 
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Another  absorption  band  of  practical  importance  is  the  4A2(F)— >4T2(F)  mid-IR  transition  around 
3000  nm.  First,  this  transition  could  be  used  to  enable  passive  Q-switching  of  Er-lasers  operating 
at  3  pm.  Second,  this  transition  is  promising  for  the  development  of  tunable  lasers  operating  over 
the  3-4  pm  range.  Third,  the  4T2(F)  upper  level  is  the  donor  level  in  Co-Fe  co-doped  gain 
elements.  However,  this  transition  is  symmetry  forbidden  for  the  electric-dipole  interaction  and 
the  maximum  value  of  the  absorption  cross  section  was  estimated  to  be  5.2x10'“  cm'“  at  2956 
nm  in  Co:ZnSe  samples.  This  is  more  than  order  of  magnitude  less  than  the  value  of  the 
absorption  cross-section  values  for  the  4A2(F)— >4Ti(F)  transition  near  1.6  pm.  The 
experimentally  measured  absorption  cross-sections  for  Co  ions  in  ZnSe  and  ZnS  crystals  on  the 
A)  4A2(F)— »4T2(P);  B)  4A2(F)— >4Ti(F);  and  C)  4A2(F)— >4T2(F)  transitions  are  depicted  in  Figures 
77  and  78,  respectively  and  are  summarized  in  Table  15. 


Table  15,  Summary  of  the  spectroscopic  characteristics  of  Co:ZnSe  and  Co:ZnS. 


Crystal 

Transition 

Wavelength, 

[nm] 

Ground  state 
absorption  cross 
section 

maximum 

cross-section 

Ref 

Co:ZnSe 

4A2(F)^4T!(P) 

768 

64x1  O'19  cm'2 

[cs] 

4A2(F)^4T!(F) 

1535 

6.5xl0'19  cm'2 

[24] 

4A2(F)^4T!(F) 

1598 

5.3xl0'19  cm'2 

[25] 

4A2(F)^4Ti(F) 

1534 

1 1.5-9. 7  xlO'19  cm'2 

[26] 

4A2(F)^4T!(F) 

1540 

7.6xl0'19  cm'2 

[27] 

4A2(F)^4T!(F) 

1573 

5. 7x1  O'19  cm'2 

[28] 

4A2(F)^4Ti(F) 

1625 

7.5xl0'19  cm'2 

[1] 

4A2(F)— >4T2(F) 

2956 

7.8xlO'20  cm'2 

[1] 

4a2(F)^4t2(F) 

2790 

llxlO'20  cm'2 

[19] 

4A2(F)^4T2(F) 

2956 

5. 2x1  O'20  cm'2 

[cs] 

Co:ZnS 

4A2(F)^4T!(P) 

726 

56xl0'19  cm'2 

[cs] 

4A2(F)^4T!(F) 

1534 

10-8  xlO'19  cm'2 

[26] 

4A2(F)^4Ti(F) 

1540 

7.7xl0'19  cm'2 

[27] 

4A2(F)^4Ti(F) 

7.5xl0'19  cm'2 

[1] 

4A2(F)— >4T2(F) 

7.8xlO'20  cm'2 

[1] 

4a2(F)^4t2(F) 

2790 

5. 6x1  O'20  cm'2 

[19] 

4A2(F)^4T2(F) 

2760 

4.8xlO'20  cm'2 

[cs] 

cs  -  current  study 
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Figure  77.  The  absorption  cross-sections  of  Co2+  ions  in  ZnSe  crystals  on  the  4A2(F)-»4T2(P) 
{Upper  part);  4A2(F)— »4T  i(F)  {Middle  part  part);  and  4A2(F)->4T2(F)  transitions  {Bottom 
part).  The  black  curve  refers  to  a  cobalt  concentration  of  Nc0=  lxlO18  cm'3  and  the  blue 
curves  refers  to  a  cobalt  concentration  of  Nc0=  2.2xl019  cm'3. 
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Figure  78.  The  absorption  cross-sections  of  Co2+  ions  in  ZnS  crystals  on  the  4A2(F)-»4T2(P) 
(Upper part)'  4A2(F)— >4T|(F)  (Middle part);  and  4A2(F)->4T2(F)  transitions  (Bottom  part). 
The  black  curve  refers  to  a  cobalt  concentration  of  Nc0  =  0.5x10  cm'  and  the  blue  curve 
refers  to  a  cobalt  concentration  of  Nc0=  1.9xl019  cm'3. 
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4.3 


Spectroscopic  characterization  of  Co2+  and  Fe2+  co-doped  samples 


The  lack  of  convenient  pump  sources  is  a  major  reason  hampering  the  broad  application  of  mid- 
IR  lasers  based  on  Fe:ZnSe  crystals.  One  possible  way  of  exciting  iron  is  via  energy  transfer 
from  other  conveniently  pumped  TM  ions.  The  donor-acceptor  energy  transfer  rate,  according  to 
Forster-Dexter  theory,  is  proportional  to  the  overlap  integral  of  the  donor  emission  and  acceptor 
absorption  bands.  Therefore,  among  TM  ions,  chromium  and  cobalt  are  very  attractive 
candidates  for  this  application.  Both  ions  have  strong  absorption  bands  in  the  1. 4-2.0  pm  spectral 
range  where  a  number  of  convenient  pump  sources  are  available.  Figure  79  shows  possible 
energy  transfer  channels  among  Fe2+,  Co2+  and  Cr2+  ions  in  II-VI  crystals.  As  one  can  see,  the 
spectral  overlap  between  cobalt  and  iron  is  practically  ideal.  To  study  energy  transfer 
characteristics,  several  cobalt  and  Co-Fe  co-doped  samples  were  fabricated. 


cr*+ 

Fe* 

CoJ' 

Figure  79.  Diagram  of  the  energy  transfer  mechanisms  in  the  iron,  chromium  and  cobalt 

co-doped  II-VI  materials. 


Initially,  a  one  week  thermal  diffusion  from  a  cobalt  thin-film  was  carried  out  to  fabricate  cobalt 
doped  ZnSe  and  ZnS.  After  that  in  order  to  fonn  Co:Fe:ZnSe  and  Co:Fe:ZnS,  iron  thin  films 
were  deposited  on  both  crystals.  To  provide  a  homogenous  distribution  of  ions  over  the  crystal 
thickness,  the  crystals  were  annealed  for  an  additional  2  weeks.  The  absorption  spectra  of  the 
Co:Fe:ZnSe  samples  are  depicted  in  Figure  80.  As  one  can  see  from  Figure  80,  the  Co  ion 
concentration  is  the  same  (Nc0  =  4.5x10  cm'  )  in  all  samples,  while  the  concentration  of  Fe“ 
changes  from  0toNFe=  18x1018  cm'3. 
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Figure  80.  Transmission  and  absorption  spectra  of  Co:Fe:ZnSe  crystals 


Similar  results  were  obtained  for  the  Co:Fe:ZnS  sample  shown  in  the  Figure  81.  As  one  can  see 
from  Figure  81,  the  absorption  bands  of  the  Fe"  and  Co  ions  are  shifted  slightly  to  shorter 
wavelengths  in  comparison  with  the  ZnSe  host  due  to  a  stronger  crystal  field  splitting.  From 

1  o 

absorption  measurements,  the  cobalt  concentration  was  estimated  to  be  Nc0~  8><10  cm'  while 

1  O  T 

the  iron  concentration  varied  from  0  to  10><10  cm'  . 

The  luminescence  and  PL  kinetics  were  measured  using  an  injection  seeded,  single  frequency  Q- 
switched  Nd:YAG  laser  (GCR-230-10,  Spectra  Physics)  with  a  maximum  output  energy  of  1.5  J 
at  1.064  pm,  a  linewidth  of  0.003  cm'1,  a  10  ns  pulse  duration  and  a  repetition  rate  of  10  Hz. 
After  passing  through  an  optical  isolator,  the  1.064  pm  Nd:YAG  radiation  was  focused  by  a  25 
cm  focal  length  lens  into  a  50  cm  long  Raman  cell  filled  with  D2  at  a  pressure  of  600  psi.  In  a 
backscattering  geometry,  the  output  energy  of  the  first  Stokes  line  at  1.56  pm  exceeded  300  mJ 
with  a  pulse  duration  of  7  ns  at  the  FWHM. 


Figure  81.  Absorption  spectra  of  Co:Fe:ZnS  crystals 
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Room  temperature  time-resolved  photoluminescence  spectra  of  the  Co:Fe:ZnSe  crystal  (with 
cobalt  and  iron  concentrations  equal  to  Cc0=  6x10  cm'  and  CFe  =  8.5x10  cm' ,  respectively) 
are  shown  in  Figure  82.  Pulsed  radiation  at  1.56  pm  with  a  15  ns  pulse  duration  was  used  for 
direct  excitation  of  Co2+  on  the  4A^— >4Ti  transition.  It  is  interesting  to  note  that  when  the  time 
delay  is  a  minimum,  the  luminescence  spectrum  (curve  i)  consists  of  both  Co  and  Fe  bands 
with  approximately  equal  intensities.  A  dip  in  the  luminescence  spectrum  at  4.3  pm  results  from 
atmospheric  CCF  absorption  in  the  experimental  setup.  After  a  40  ps  time-delay,  the 
luminescence  spectrum  reveals  only  a  Co"  photoluminescence  signal  (ii).  This  indicates  that  the 
most  efficient  energy  transfer  occurs  only  directly  after  excitation  and  that  the  energy  transfer 
rate  from  the  4T2  cobalt  level  is  not  significant  for  the  TM  concentrations  studied. 


wavelength  (nm) 

Figure  82.  Time  resolved  spectra  of  Co:Fe:ZnSe  under  1.56  pm  excitation  with  a  4  ps  gate- 

width  and  0  ps  (i)  and  40  ps  time  delays  (ii). 


4.4  Energy  transfer  in  Co:Fe:ZnS/ZnSe  crystals  under  excitation  at  720  nm 

Energy  transfer  studies  were  associated  with  measurements  of  the  mid-IR  spectra  and  the 
kinetics  of  photoluminescence  in  several  cobalt  iron  co-doped  samples.  The  absorption  spectra  of 
Co:Fe:ZnS  (J4510)  are  shown  in  the  Figure  83. 
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Figure  83.  Absorption  spectrum  of  the  Co:Fe:ZnS  sample  measured  at  RT  (red  curve)  and 

T  =  16°K  (blue  curve). 

18  3  19  3* 

Average  concentrations  of  cobalt  and  iron  ( Nco  =  1.8x10  cm";  NFe  =  3.5x10  cm'  )  in  the 
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samples  were  estimated  using  the  absorption  cross-sections  for  cobalt  (aab=  7.8x10'  cm'  )  and 
iron  (aab  =  6.7x10'  cm'')  at  1.5  and  2.7  pm,  respectively.  The  right  hand  side  of  Figure  83 

shows  room  and  low  temperature  absorption  on  the  4A2(F)— >4T2(P)  transition.  As  one  can  see 
from  Figure  83,  the  long-wavelength  edge  of  the  absorption  band  is  shifted  to  726  nm  at  low 
temperature.  Therefore  the  excitation  wavelength  should  be  between  670  and  725  nm  for 
effective  excitation  of  mid-IR  emission. 

To  demonstrate  the  feasibility  of  exciting  the  mid-IR  photoluminescence  on  the  4A2(F)— >4T2(P) 
transition,  we  studied  mid-IR  photoluminescence  kinetics  and  spectra  using  radiation  from  an 
Alexandrite  laser  for  optical  pumping  of  the  Co  ions.  We  used  a  slightly  focused  Alexandrite 
laser  (PAL  101,  Light  Age,  Inc.)  beam  having  a  maximum  pulse  energy  of  80  mJ  at  720-780  nm, 
a  beam  diameter  of  2  mm  (FWHM),  a  pulse  duration  of  70  ns,  and  a  repetition  rate  of  13  Hz. 

The  kinetics  of  the  mid-IR  photoluminescence  signal  of  a  Co:Fe:ZnS  crystal  under  excitation  by 
an  Alexandrite  laser  were  measured  at  T  =  20  and  77°K  and  are  shown  in  Figure  84.  Curve  i  in 
the  upper  part  of  the  figure  shows  the  low  temperature  kinetics  of  non-selective  mid-IR 
luminescence  of  the  Co:Fe:ZnS  (J4510)  sample.  As  one  can  see  from  the  Figure  84,  the  kinetics 
shown  in  curve  i  features  a  non-exponential  behavior.  On  a  long  time  scale  (t  >  300  ps),  the 
kinetics  asymptotically  approach  the  decay  curve  of  Co2+  kinetics  in  a  Co:ZnS  (J4110,  NCo  = 
1.8x10  cm'  )  crystal  with  a  decay  time  of  x  =  1.2  ms  (curve  ii).  Strong  non-exponential  decay 
results  from  the  photoluminescence  and  excitation  migration  over  Fe2+  ions  in  the  initial  stage  of 
the  kinetics.  For  comparison,  the  insert  in  Figure  84  shows  the  kinetics  (curve  iii)  of  a  Fe:ZnS 
(J3310,  NFe=  5x10  cm'  )  crystal  under  direct  optical  excitation  at  2.8  pm  with  a  decay  time  of  x 
=  6  ps  at  T  =  20°K.  The  rise  time  of  the  kinetics  is  shorter  than  1  ps  and  is  limited  by  the 
temporal  resolution  of  the  registration  system.  Similar  mid-IR  photoluminescence  dynamics 
were  measured  at  T  =  77°K  and  are  shown  in  the  bottom  part  of  Figure  84. 
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Figure  84.  Kinetics  of  the  mid-IR  photoluminescence  of  Co:Fe:ZnS  (curve  i);  Co:ZnS 
(curve  ii);  and  Fe:ZnS  (curve  iii)  at  T  =  20°K  and  77°K  under  Co2+  excitation  (curve  i  and 
ii)  at  720  nm;  and  under  Fe2+  excitation  (curve  iii)  at  2.8  pm. 

The  kinetics  of  the  mid-infrared  photoluminescence  signal  of  Co:Fe:ZnSe  under  excitation  by  an 
Alexandrite  laser  were  measured  at  T  =  14°K  and  are  shown  in  Figure  85.  For  comparison 
purposes,  the  kinetics  of  the  Fe~  mid-infrared  photoluminescence  signal  in  ZnSe  having  a  decay 
time  of  ~  50  ps  is  also  presented  in  the  Figure  85.  As  one  can  see  from  Figure  85,  the  decay 
curve  of  the  Co:Fe:ZnSe  (Nc0  =  6.4x10  cm";  Npe  =  1.9x10  cm'  )  sample  kinetics  reveals  a 
strong  non-exponential  behavior  which  asymptotically  approaches  the  decay  curve  of  the 
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Fe:ZnSe  (J1209,  NCo=  9x10  cm'  )  crystal  kinetics.  It  is  interesting  to  note  that  in  contrast  to  the 
kinetics  of  the  Co:Fe:ZnS  crystal,  the  kinetics  of  the  Co:Fe:ZnSe  crystal  does  not  have  the  long 

1  o 

decay  time  typical  of  a  single-impurity  center.  The  kinetics  of  the  CoiZnSe  (J21 10,  Nco  =6x10 
cm'3)  crystal  under  the  same  excitation  conditions  is  also  shown  in  Figure  85.  The  much  longer 
decay  time  of  the  Co:ZnSe  curve  means  that  the  photo  luminescence  of  “practically  all  of  the 
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Co  ions”  is  quenched  by  the  energy  transfer  to  Fe  ions  in  spite  of  the  fact  that  the  iron 
concentration  in  the  Co:Fe:ZnS  crystal  was  higher.  One  possible  way  of  explaining  this  is  that 
the  formation  of  the  Fe-Co  pair  center  could  be  more  efficient  in  the  ZnSe  crystal  than  in  the  ZnS 
crystal.  One  example  of  a  system  with  strong  impurity  pairing  is  the  Nd:CaF2  system  where  Nd- 
Nd  pair  formation  is  the  predominant  process  even  at  small  impurity  concentrations. 


Figure  85.  Kinetics  of  mid  infrared  photoluminescence  of  Co:Fe:ZnSe  (curve  i);  Co:ZnSe 
(curve  ii);  and  Fe:ZnS  (curve  iii)  under  760  nm  excitation  of  Co2+  (curves  i  and  ii);  and  2.8 
pm  excitation  of  Fe2+  excitation  (curve  iii)  at  T  =  14°K. 


Shown  in  Figure  86  are  the  photoluminescence  spectra  of  Co:Fe:ZnSe  and  CoiZnSe  samples 
under  direct  excitation  of  the  4A2(F)— >4Ti(P)  (-750  nm)  transition  in  cobalt  at  14°K.  In  these 
experiments,  we  used  a  crystal  with  Fe  and  Co  concentrations  of  19x10  cm'  and  6.4x10  cm' 
3,  respectively.  In  the  Co.'ZnSc  crystal,  one  can  see  a  well-structured  photoluminescence 
spectrum  corresponding  to  the  4T2(F)  — »4A2(F)  symmetry-forbidden  transition  of  cobalt  with  a 
luminescence  maximum  at  3040  nm.  In  the  Co:Fe:ZnSe  co-doped  samples,  the 
photoluminescence  shown  in  Figure  86B  was  significantly  shifted  to  longer  wavelengths  and 
coincides  with  earlier  reported  results  of  the  low  temperature  photoluminescence  spectrum  of 
Fe2+  ions.  A  maximum  in  the  photoluminescence  signal  of  the  iron-cobalt  co-doped  samples  was 
found  to  occur  at  3750  nm.  Also,  as  one  can  see  from  Figure  86,  the  direct  photoluminescence 
signal  from  cobalt  was  significantly  quenched  and  was  hardly  detectable  in  comparison  with  the 
iron  photoluminescence  signal. 
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Figure  86.  Mid-infrared  luminescence  spectra  of  iron  and  cobalt  co-doped  ZnSe  crystals  at 

T  =  14°K  under  760nm  excitation. 


4.5  Spectroscopic  characterization  of  Co2+  and  Ni2+  co-doped  samples 

The  absorption  spectrum  of  an  unpolished  Co:Ni:ZnSe  sample  (811-65)  is  shown  in  Figure  87 
below.  For  comparison,  the  absorption  spectrum  of  a  Co.'ZnSe  crystal  is  also  presented  in  Figure 
87.  Two  additional  bands  at  870  and  1130  nm  were  observed  due  to  the  ’Ti(F)  — »JTi(P)  and 
Ti(F)— »  A2(F)  transitions,  respectively,  of  nickel  in  the  nickel-cobalt  co-doped  crystal. 
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Figure  87.  Absorption  spectra  of  the  Co:Ni:ZnSe  (black  curve),  and  Co:ZnSe  (red  curve) 

crystals. 
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The  photoluminescence  spectra  of  Ni:Co:ZnS  and  Ni:Co:ZnSe  were  measured  at  room 
temperature  under  excitation  of  a  1064  nm  free  running  Nd:YAG  laser.  While  the  Ni:Co:ZnSe 
sample  did  not  feature  a  mid-infrared  luminescence,  a  characteristic  Co  luminescence  was 
observed  in  the  Ni:Co:ZnS  sample  (see  Figure  88). 


Figure  88.  RT  photo-luminescence  spectra  of  Ni:Co:ZnS  (black  curve)  and  Co:ZnS  (red 

curve)  samples. 

However,  a  similar  PL  band  was  also  observed  in  a  Co:ZnS  sample  under  1064  nm  excitation. 
This  luminescence  is  not  believed  to  be  due  to  energy  transfer  from  Ni  to  Co,  but  rather  to  be  due 
to  a  direct  excitation  into  the  short  wavelength  tail  of  the  absorption  band  of  Co2+. 
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5.0  LASER  CHARACTERIZATION  OF  TM  DOPED  SAMPLES 


5.1.  Lasing  of  Co:Fe:ZnSe  and  Co:Fe:ZnS  crystals  via  effective  Co— »Fe  energy  transfer 
under  1.56  pm  wavelength  excitation 

The  high  gain  associated  with  highly  doped  Co:Fe:ZnS  crystals  allowed  us  to  test  the  lasing 
properties  of  these  crystals  by  pumping  at  1.56  pm  even  in  the  absence  of  an  external  cavity.  For 
this  experiment,  we  used  Co:Fe:ZnS/ZnSe  crystals  (10x10x1.2  mm3)  having  cobalt  and  iron 
concentrations  of  0.5  x  10  cm'  and  10x10  cm'  ,  respectively. 

Pump  radiation  from  a  Q-switched  Raman  shifted  YAG  laser  at  1.56  pm  was  slightly  focused 
into  a  spot  having  a  -3  mm  diameter.  The  crystal  was  placed  in  a  Close-Cycle  Refrigerator 
operating  in  the  14-300°K  temperature  range. 

Figure  89  shows  the  luminescence  spectra  of  several  crystals  measured  at  14°K.  The  first  plot 

2_|_ 

involving  Co:ZnS  shows  the  typical  Co”  luminescence  which  has  characteristic  peaks  around 
2800  and  3100  mn  upon  excitation  at  1.56  pm.  The  third  plot  shows  the  luminescence  spectrum 
of  a  Co:Fe:ZnS  crystal  having  a  low  iron  concentration  of  0.5X10  cm"  under  the  same 
excitation  with  1.56  pm.  As  one  can  see  from  this  graph,  in  addition  to  the  Co2+  luminescence, 
there  is  also  a  strong  Fe  luminescence  having  characteristic  peaks  at  3440  mn,  3570  nm,  and 
3820  nm.  For  comparison  purposes,  the  luminescence  spectrum  of  an  Fe:ZnS  sample  under 
direct  excitation  of  the  iron  ions  at  2.9  pm  was  measured  at  14°K  and  is  depicted  in  the  second 
graph.  The  fourth  graph  shows  the  luminescence  spectrum  of  a  Co:Fe:ZnS  crystal  under  1.56  pm 
excitation  having  an  iron  concentration  x50  times  greater.  As  one  can  see,  there  is  no  Co2+ 
luminescence  in  this  crystal.  The  measured  luminescence  spectrum  shows  only  Fe“  bands  that 
are  slightly  red-shifted  in  comparison  with  the  spectrum  of  the  low  doped  sample.  This  result 
proves  that  there  is  a  fast  energy  transfer  from  Co  to  Fe  ions  in  this  crystal.  We  used  this 
crystal  for  the  following  laser  experiments. 
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Figure  89.  Mid-infrared  luminescence  spectra  of  iron  and  cobalt  doped  crystals  at  T  = 

14°K. 

Two  basic  experiments  on  the  pump  energy  dependence  of  the  luminescence  kinetics  and 
luminescence  spectra  were  perfonned  to  demonstrate  lasing  of  Fe  via  effective  energy  transfer 
from  the  Co2+  ion.  Positive  feedback  for  the  oscillation  of  Fe  was  enabled  by  the  Fresnel 
reflections  from  the  crystal  facets  which  provided  a  reflectivity  of  18%  from  each  surface.  The 
large  Fresnel  Number  NF  —  d2 /LA  ~104  of  the  resonator  provided  for  minimum  diffraction 
losses  in  the  cavity. 

The  results  of  the  experiments  are  shown  in  the  Figure  90.  Figure  90A  shows  the 
photoluminescence  spectra  of  a  Co:Fe:ZnSe  crystal  for  different  pump  energies.  As  one  can  see 
in  Figure  90A,  at  low  pump  energy,  the  measured  PL  spectra  were  typical  for  the  luminescence 
of  iron.  The  profile  of  the  photoluminescence  spectrum  also  demonstrated  a  threshold  behavior 
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that  was  dependent  on  the  pump  energy.  A  stimulated  emission  band  around  3850  nm  as  shown 
in  Figure  90  appeared  when  the  threshold  was  exceeded. 

The  change  of  the  luminescence  profile  above  the  laser  threshold  was  accompanied  by  a 
shortening  of  the  decay  time  of  the  luminescence  kinetics  as  shown  in  the  Figure  90B.  When 
pumped  at  low  energy,  one  can  see  a  low  temperature  Fe"  kinetic  with  a  luminescence  rise  time 
of  t  ~  1  ps  followed  by  an  exponential  decay  with  a  lifetime  of  t  ~  2.5  ps.  Above  the  laser 
threshold,  the  kinetics  demonstrate  a  typical  laser  spike  with  a  ~200  ns  pulse  duration  followed 
by  luminescence  decay  kinetics.  The  threshold  dependence  of  the  emission,  a  significant  line 
narrowing,  and  the  shortening  of  the  lifetime  are  clear  evidence  of  Co:Fe:ZnSe  lasing  via  the 
Co— >Fe  energy  transfer  mechanism. 
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Figure  90.  A)  Mid-infrared  luminescence  spectra  of  a  Co:Fe:ZnS  crystal  under  different 
pump  energies  at  T  =  14°K,  B)  photoluminescence  kinetics  of  a  Co:Fe:ZnS  crystal  for 
different  pump  energies  below  (red)  and  above  (pink)  the  laser  threshold. 
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5.2.  Lasing  of  Co:Fe:ZnSe  and  Co:Fe:ZnS  crystals  via  effective  Co— >Fe  energy  transfer 
under  0.76  pm  wavelength  excitation 

It  was  shown  in  the  previous  section  that  the  strong  absorption  band  between  650  and  800  mn  in 
cobalt  doped  ZnS  and  ZnSe  crystals  could  be  utilized  for  the  optical  excitation  of  cobalt  by 
commercially  available  high  power  laser  diodes  operating  in  this  spectral  range.  The  high 
concentration  of  co-dopants  in  Co:Fe:ZnSe  and  Co:Fe:ZnS  crystals  as  well  as  the  fast  rate  of 
energy  transfer  will  provide  the  necessary  high  gain  to  enable  lasing  under  pumping  at  0.76  pm 
without  the  need  for  an  external  cavity.  The  Co:Fe:ZnSe  crystal  used  in  these  experiments  had  a 
cobalt  concentration  of  6.4x10  cm'  and  an  iron  concentration  of  19x10  cm'  ,  while  the 
Co:Fe:ZnS  crystal  that  was  used  had  a  cobalt  concentration  of  8x10  cm'  and  an  iron 
concentration  of  10xl018  cm'3. 

In  the  laser  experiments,  the  pump  radiation  of  an  alexandrite  laser  was  slightly  focused  to  a  spot 
diameter  of  ~2  mm  and  the  Co:Fe:ZnSe  crystal  was  pumped  at  760  mn.  The  crystal  was 
contained  in  a  closed-cycle  refrigerator  to  enable  adjustment  of  the  temperature  from  14  to 
300°K.  The  luminescence  spectra  and  kinetics  were  measured  over  a  range  of  pump  energies  to 
show  threshold  behavior,  Figure  89.  Feedback  in  the  crystal  was  provided  by  Fresnel  reflections 
of  -18%  from  each  of  the  polished  facets. 

Figure  91  shows  the  photoluminescence  spectra  of  the  Co:Fe:ZnSe  crystal  for  different  pump 
energies.  As  one  can  see,  at  low  pump  energy,  the  measured  PL  spectra  were  typical  for  the 
luminescence  of  iron  ions.  The  dependence  of  the  PL  spectrum  profile  on  the  pump  energy 
demonstrated  a  threshold  behavior  accompanied  by  the  appearance  of  a  stimulated  emission  band 
around  3.85  pm.  The  laser  threshold  was  measured  to  be  19  mJ  of  pump  power. 


Wavelength  (|xm) 

Figure  91.  Mid-infrared  luminescence  spectra  of  the  Co:Fe:ZnSe  crystal  for  different  pump 
energies  under  760  nm  excitation  at  T  =  14°K  measured  below  and  above  the  laser 

threshold. 
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The  change  of  the  luminescence  profile  above  the  laser  threshold  was  accompanied  by  a 
shortening  of  PL  kinetics  as  shown  in  Figure  92.  For  low  pump  energies,  one  can  see  the  low 
temperature  Fe“  kinetics  with  a  luminescence  rise  time  ~0.5  ps  followed  by  an  exponential 
decay.  Above  the  laser  threshold,  the  kinetics  exhibit  the  typical  laser  spike  with  a  -100  ns  pulse 
duration  followed  by  below  threshold  luminescence  kinetics.  The  temporal  profile  of  the  laser 
pulses  were  obtained  by  taking  the  difference  between  the  measured  kinetics  and  the  normalized 
luminescence  signals  measured  below  the  laser  threshold.  The  results  are  shown  in  Figure  92B. 
As  one  can  see  from  Figure  92,  an  increase  in  the  pump  energy  results  in  a  shortening  of  the 
onset  of  lasing  as  well  as  the  pulse  duration.  Relaxation  spikes  with  -100  ns  between  laser  spikes 
were  observed  at  high  pump  energy.  The  threshold  dependence  of  emission,  a  significant  line 
narrowing,  as  well  as  the  shortening  of  the  lifetime  are  clear  evidence  of  Co:Fe:ZnSe  lasing  via 
the  Co— ►Fe  energy  transfer  mechanism. 


Figure  92.  A)  Co:Fe:ZnSe  crystal  photoluminescence  kinetics  for  different  pump  energies 
at  an  excitation  wavelength  of  760  nm  measured  below  and  above  the  laser  threshold  at 
14°K.  B)  Difference  between  the  measured  kinetics  and  the  normalized  luminescence 
signals  measured  below  the  laser  threshold. 
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5.3.  Preliminary  laser  experiments  with  Co:ZnSe/ZnS  crystals 

Preliminary  laser  experiments  were  perfonned  with  Co:ZnSe  and  Co:ZnS  Brewster  cut  samples 
under  transverse  pumping  by  a  free  running  Alexandrite  laser  at  739.7  mn.  The  chosen  samples 
are  of  the  dimensions  and  concentrations  shown  in  Table  16. 

Table  16.  Properties  of  Co:ZnSe  and  Co:ZnS  samples  studied  by  pumping  with  a  739.7  nm 

laser. 


Length 

(mm) 

Width 

(mm) 

Height 

(mm) 

Co  Concentration 
[cm3] 

Co:ZnS  (#811-14-06) 

37 

6 

2.5 

1.76xlOiy 

Co:ZnSe  (#811-29-11) 

40 

6.8 

2.4 

2.15xlOiy 

RT  mid-infrared  PL  measurements  at  739.7  nm  under  free  running  excitation  of  each  sample  are 
shown  in  Figure  93. 


wavelength,  nm 

Figure  93.  Middle-Infrared  PL  spectra  of  Co:ZnS  and  Co:ZnSe. 


Figure  94  shows  the  RT  photo  luminescence  kinetics  of  Co:ZnS  and  Co:ZnSe  under  excitation  of 
the  4A2(F)^4T2(P)  transition.  As  one  can  see  from  Figure  94,  the  luminescence  lifetime 
constants  are  equal  to  208  and  167  ps  at  RT  for  the  Co:ZnSe  and  Co:ZnS  crystals,  respectively. 
The  kinetic  rise-time  constants  were  measured  to  be  less  than  the  response  time  of  the  detection 
system  (<0.5  ps).  These  results  indicate  that  the  upper  laser  level  should  accumulate  a  population 
and  thereby  create  an  inversion  even  under  the  relatively  long  free  running  pulses  of  the 
Alexandrite  laser.  Let’s  estimate  the  possible  optical  gain  of  Co:ZnSe/ZnSe  active  elements  on 
the  4 A2(F) — >4T2(F)  transition  under  excitation  with  100  mJ  (Np h  =  3.7xl017  photon  @  0.74  pm) 
pulses  from  an  Alexandrite  laser.  We  will  assume  that  quantum  efficiency  of  the  pump  is  equal 
to  1 .  By  using  a  cylindrical  lens  in  addition  to  taking  advantage  of  the  strong  absorption  of  the 
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0.1 


Figure  94.  PL  kinetics  of  Co:ZnS  and  Co:ZnSe  crystals 


pump  radiation  from  an  Alexandrite  laser  (k  >  30  cm'1),  one  can  design  a  pump  transfer 
geometry  with  an  active  channel  area  of  A  =  300x300  pm.  In  this  geometry,  the  single  pass  gain 
could  be  estimated  as: 


G  =  exp (Nphaem/A)  ~  exp(7)  (70) 

20  2 

with  cTem  ~  gi/g2'Oah  ~  1.7x10'“  cm'“.  The  laser  performance  of  the  Co:ZnSe/ZnS  gain  elements 
was  tested  in  the  experimental  setup  shown  in  Figure  95.  In  this  setup,  the  output  coupler  had  a 
reflectivity  R  >  98%  for  wavelengths  longer  than  3.5  pm.  In  addition,  the  Alexandrite  pump  laser 
operated  in  the  free  running  regime  with  an  average  power  of  1.4  W  with  a  repetition  rate  of  13 
Hz  at  739.7  nm. 


Alexandrite 

Laser 


Wave 

Meter 


Figure  95.  Laser  Experimental  Setup. 
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Unfortunately,  under  these  experimental  conditions  Co:ZnS/ZnSe  lasing  was  not  detected.  There 
are  two  major  factors  that  could  have  prevented  lasing  at  RT.  First  of  all,  it  can  be  due  to 
possible  excited  state  absorption  from  the  upper  laser  level  (4T2~^4Ti).  The  second  reason  could 
have  to  do  with  the  branching  ratio  of  relaxation  processes  from  the  pumped  "^(P)  level.  In 
addition  to  relaxation  to  the  upper  laser  level  4T2(P)  — >  4  73  (F),  relaxation  from  the  pumped  level 
could  also  include  relaxation  to  the  ground  level  4T2(P)  — »  4Aj  (F).  The  total  pump  efficiency 
depends  on  this  ratio. 


5.4.  Co:ZnSe  saturable  absorber  for  passively  Q-switched  Er:YAG  lasers 

Er:YAG  lasers  operating  in  the  1.6pm  eyesafe  spectral  range  (4Ii3/2  — >41 1 5/2  transition)  are  widely 
used  in  a  variety  of  medical,  scientific,  and  military  applications.  Direct  resonant  pumping  into 
the  41 13/2  manifold  by  commercially  available  1.532pm  Er-fiber  lasers  results  in  a  small  quantum 
defect  of  -7%  as  well  as  a  high,  more  than  80%,  efficient  laser.  Cobalt  and  chromium  doped  II- 
VI  crystals  have  been  successfully  utilized  as  passive  saturable  absorbers  due  to  the  absence  of 
excited  state  absorption  as  well  as  having  a  high  cross-section  of  saturation  of  -  10'  cm”.  The 
simplicity  and  compactness  of  TMTI-VI  passive  Q-switches  is  ideal  for  many  applications. 
However,  a  1.6  pm  Er-laser  passively  Q-switched  by  a  Co  doped  II-VI  material  under 
continuous  wave  (CW)  excitation  was  not  reported  yet.  The  purpose  of  this  work  was  to  verify 
whether  or  not  a  Co:ZnS  crystal  having  a  relatively  long  upper  state  relaxation  time  (-200  ps) 
could  be  effective  as  a  passive  Q-switch  in  erbium  lasers. 

The  laser  cavity  configuration  used  in  this  experiment  is  depicted  in  Figure  96.  A  1.1  mm  pump 
beam  from  an  ELR-20-1532.6-LP  fiber  laser  (IPG  Photonics  Corp.)  was  focused  into  an  Er:YAG 
laser  rod  (0.25  %  at.  Er:YAG,  5  mmx50  mm)  through  a  dichroic  flat  45-degree  folding  mirror 
(HR  @  1600-1700  nm,  AR  @  1532  nm)  with  a  300  mm  AR-coated  focusing  lens.  The  laser 
scheme  is  based  on  a  semiconcentric  cavity  design  with  a  250  mm  end  mirror  and  an  80%  output 
coupler.  ZnS  doped  with  cobalt  having  an  initial  transmission  of  -92%  at  1617  nm  served  as  a 
saturable  absorber  and  was  placed  at  Brewster’s  angle  50  mm  away  from  the  dichroic  mirror. 
The  transmission  spectrum  of  the  Co:ZnS  crystal  utilized  is  shown  in  Figure  97  at  normal 
incidence  . 


FM 


Er-Fiber 

Laser 

1532.6  nm 


0.25  at%  Er:YAG  05x50  mm 


Co:ZnS 


OC 

R=80% 


Back  Mirror 
R=250  mm 


Figure  96.  Experimental  setup. 
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wavelength,  nm 

Figure  97.  Transmission  spectrum  of  the  Co:ZnS  crystal. 


In  Figure  98  are  shown  the  input-output  characteristics  of  the  passively  Q-switched  Er:YAG 
laser  with  an  initial  transmission  of  the  Co:ZnS  passive  crystal  of  92%.  The  output  power 
reached  a  level  of  2.1  W  when  pumped  with  16  W  of  incident  power.  The  slope  efficiency  was 
24%  slope  with  a  real  optical  efficiency  of  14%.  For  comparison,  CW  oscillations  of  the  Er:YAG 
laser  were  measured  to  have  a  38%  slope  efficiency  in  the  same  cavity  without  the  Co:ZnS 
crystal.  The  slope  efficiency  for  operation  in  the  passively  Q-switched  regime  relative  to  CW 
operation  was  as  high  as  63%. 


Figure  98.  Output  power  of  a  passively  Q-switched  and  CW  Er:YAG  laser. 


Figure  99  represents  a  characteristic  train  of  Q-Switched  Er:YAG  pulses  obtained  from  the  laser 
with  an  input  power  of  7.4  W.  As  one  can  see  for  this  pump  power,  the  passively  Q-switched 
Er:YAG  laser  operates  at  repetition  rate  of  1.25  kHz  and  features  a  stable  output  pulse  energy  of 
0.1  mJ. 
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time,  ms 

Figure  99.  Train  of  output  pulses  from  a  passively  Q-switched  Er:YAG  laser  at  a  pump 

power  of  7.4  W. 


According  to  a  simple  model  of  the  passively  Q-switched  regime  of  operation,  the  pulse  energy 
is  defined  by  the  initial  transmission  of  the  absorber  and  is  not  sensitive  to  a  variation  in  the 
pump  power.  As  one  can  see  from  Figure  100A,  the  output  energy  of  a  passively  Q-switched 
Ho:YAG  laser  is  approximately  stable  with  increases  in  the  pump  power  to  more  than  30%  of  the 
laser  threshold.  Figure  100B  shows  the  repetition  rate  of  a  passively  Q-switched  Er:YAG  laser 
versus  the  incident  pump  power.  The  maximum  estimated  repetition  rate  was  38  kHz. 
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Figure  100.  Output  energy  (A)  and  repetition  rate  (B)  of  the  passively  Q-switched  Er:YAG 

laser. 


Finally,  oscillation  pulses  of  duration  0.5- 1.2  ps  from  the  passively  Q-switched  cavity  showed  no 
evidence  of  the  mode  beating  associated  with  multimode  operation  (see  Figure  101.) 
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Figure  101.  Temporal  profile  of  the  passively  Q-switched  output  pulses. 
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6.0  CONCLUSIONS 


The  technology  of  high  optical  quality  gain  elements  based  on  Co,  Ni,  and  Fe  doped  ZnSe  and 
ZnS  crystals  which  are  promising  for  room-temperature  lasing  was  developed.  Multiple  co¬ 
doped  ZnSe  and  ZnS  crystals  have  been  fabricated  and  studied  via  absorption  and 
photoluminescence  spectroscopy.  For  the  first  time,  lasing  was  achieved  on  Fe  ions  in  ZnS  and 
ZnSe  crystals  co-doped  with  Co  via  effective  Co-Fe  energy  transfer.  For  the  first  time,  Co:ZnS 
crystals  were  successfully  used  for  passive  Q-switching  of  CW  pumped  Er:YAG  1645  nm  lasers. 

We  developed  a  theory  to  describe  the  properties  of  TMs  in  II-VI  materials  taking  into  account 
all  of  the  operators  necessary  to  describe  the  ion  interactions  associated  with  3dn  electronic 
configurations  (n  =  2,  3,  4,  6,  7,  and  8)  to  include  the  electrostatic  and  spin-orbit  interactions  as 
well  as  the  influence  of  the  crystal  field  and  external  magnetic  field  on  the  electronic  energy.  For 
the  ions  Co2+,  Fe2+,  and  Ni2+,  we  derived  explicit  formulas  for  the  energy  interaction  matrices. 

We  constructed  a  system  for  calculating  the  eigenvalues  of  the  energy  matrices  based  on  the 
quantum  numbers  of  an  orthogonal  group  of  three-dimensional  rotation  states  of  the  ions  as  well 
as  two  electrostatic  parameters  (F22,  F44,  (B,  C)),  the  spin-orbit  parameter  ^3d,  plus  the  13 
parameters  Akq  (k  =  2,  4)  which  describe  a  crystal  field  of  any  symmetry. 

The  matrices  which  were  calculated  are  a  modem  alternative  of  the  Tanabe-Sugano  matrices, 
which  happen  to  be  a  particular  case  of  the  matrices  that  we  developed,  i.e.,  the  case  where  -Ad  =0 
and  the  crystal  field  has  cubic  symmetry  with  a  single  parameter  Dq  that  depends  on  the  strong 
crystal  field  wave-function  basis. 

We  created  matrices  to  calculate  the  intensities  of  electro-dipole  transitions  arising  from  odd 
potentials  of  a  crystal  field  having  tetrahedral  symmetry.  The  matrices  enabled  calculation  of  the 
oscillator  strengths  and  lifetimes  for  all  transitions  of  the  3dn  (n  =  2,  3,  4,  6,  7,  and  8)  electronic 
configuration.  In  addition,  they  allow  the  computation  of  the  oscillator  strengths  and  lifetimes  of 
all  transitions  in  the  absorption  and  fluorescence  spectra,  as  well  as  fluorescence  and  absorption 
spectra  for  transitions  to  and  from  excited  states.  We  further  developed  a  method  to  find  the 
intensity  parameters  from  the  experimental  values  of  the  fluorescence  transition  lifetimes.  These 
values  are  then  used  to  calculate  the  oscillator  strengths  and  lifetimes  of  all  the  other  transitions 
in  spectra  of  various  types. 

We  developed  the  principles  for  modeling  spectra.  To  detennine  the  bandwidths  of  the 
transitions  that  identify  the  spectra  of  3d11  (n  =  2,  3,  4,  6,  7,  and  8)  ions  in  semiconductor  crystals, 
we  considered  the  oscillator  strengths  versus  the  transition  energies.  We  coordinated  the 
experimental  and  theoretical  results  by  fitting  the  modeled  spectra  to  the  experimental  data. 

We  created  a  method  to  calculate  the  configuration  curves  of  the  energy  levels  taking  into 
account  the  vibration  amplitudes  of  the  surrounding  atoms.  We  developed  energy  and  intensity 
matrices  and  approximated  the  crystal  field  parameters  of  the  zeroth  and  fourth  order  by  the 
parabolic  vibration  amplitude  dependencies.  Configuration  diagrams  are  calculated  for  the  Co2+, 
Fe"  ,  and  Ni"  ions  and  as  well  as  various  anharmonic  cases. 

The  theory  that  was  created  was  used  to  analyze  experimental  data  from  the  literature  for  Co  , 

2_|_  •2-b 

Fe  ,  an  Ni"  ions  in  ZnSe  and  good  agreement  was  obtained  between  the  experimental  and 
modeled  spectra. 

To  conclude,  all  of  the  theory  and  application  methods  were  originally  developed. 
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APPENDIX 


A.l  Program  for  calculation  of  the  spectra  of  TM2+:ZnSe 

The  program  for  calculating  the  spectra  of  Co2+,  Fe2+,  and  Ni2+  in  ZnSe  is  realized  in  the  Maple 
package  with  the  following  sequence  of  operations.  First,  an  array  of  states  (STS)  states  is 
formed.  An  example  is  the  array  of  Co2+:  ZnSe  states: 


STS=([(4F)4A2(G’),  (4F)4A2(G’),  (4F)4A2(G’),  (4F)4A2(G’),  (4F)4T2,  (4F)4T2, 
(4F)4T2,  (4F)4T2,  (4F)4T2,  (4F)4T2,  (4F)4T2,  (4F)4T2,  (4F)4T2,  (4F)4T2,  (4F)4T2, 
(4F)4T2,  (4F)4Th  (4F)4T1?  (4F)4T1?  (4F)4T!,  (4F)4Th  (4F)4Tj,  (4F)4Th  (4F)4T,, 
(4F)4T1?  (4F)4T1?  (4F)4T1?  (4F)4T1?  (2G)2E,  (2G)2E,  (2G)2E,  (2G)2E,  (2G)2Tj, 
(2G)2T1?  (2G)2Ti,  (2G)2Th  (2G)2Ti,(2G)2Th  (4P)4Th  (4P)4Th  (4P)4T,,  (4P)4T!, 
(4P)4T1?  (4P)4T1?  (4P)4T1?  (4P)4T1?  (4P)4T,,  (4P)4T1?  (4P)4Tj,  (4P)4Tj,  (2G)2A1? 
(2G)2Ai,  (2G)2T2,  (2G)2T2,  (2G)2T2,  (2G)2T2,  (2G)2T2,  (2G)2T2,  (2H)2T2,  (2H)2T2, 
(2H)2T2,  (2H)2T2,  (2H)2T2,  (2H)2T2,  (2H)2Ti,  (2H)2TI;  (2H)2T1?  (2H)2T1;  (2H)2T1? 
(2H)2T1?  (2H)2E,  (2H)2E,  (2H)2E,  (2H)2E,  (2P)2T2,  (2P)2T2,  (2P)2T2,  (2P)2T2, 
(2P)2T2,  (2P)2T2,  (2H)2Th  (2H)2T1?  (2H)2T1?  (2H)2Th  (2H)2Th  (2H)2Th  (2D)2T2, 
(2D)2T2,  (2D)2T2,  (2D)2T2,  (2D)2T2,  (2D)2T2,  (2D)2E,  (2D)2E,  (2D)2E,  (2D)2E, 
(2F)2A2,  (2F)2A2,  (2F)2T2,  (2F)2T2,  (2F)2T2,  (2F)2T2,  (2F)2T2,  (2F)2T2,  (2F)2T1? 
(2F)2T!,  (2F)2T!,  (2F)2T1?  (2F)2T1?  (2F)2Th  (2D)2E,  (2D)2E,  (2D)2E,  (2D)2E, 
(2D)2T2,  (2D)2T2,  (2D)2T2,  (2D)2T2,  (2D)2T2,  (2D)2T2]). 


(71) 


We  also  formed  an  array  of  experimental  energy  level  values.  Cited  below  is  an  example  of  the 

2d-  •  •  1 

experimental  energy  level  values  for  Co“  :ZnSe  in  eV  translated  into  cm'  : 


ZnSelE 

0.43,0.43,0.75,0.75,0.75,0.75,0.75,0.75,0.83,0.83,0.83,0.83,0.83,0.83,0,0, 

0,0,0,10,0,0,0,0,1.64,1.64,1.64,1.64,1.64,1.64,1.71,1.71,1.71,1.71,1.78,1.78, 

0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2.36,2.36,2.43,2.43, 

2.43,2.43,2.55,2.55,2.55,2.55,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 

0,0,0,0,0]). 


(72) 


A.2  Matrix  of  energy  interactions  for  iron  group  ions 

The  program  outlined  below  calculates  the  matrix  elements  of  interaction. 

1 .  The  factorials  are  calculated  as: 

>  F(j,k);  if  j<k  then  0  else  (j -k) ! 

>  FS(ji,j2,j3,li,l2,l3,k)=  k!*F(ji+j2-j3,k)*F(li+l2-j3,k)*F(ji+l2-l3,k)*F(l1+j2- 

l3,k)*F(j  3+13+kj  i+li)*F(j  3+l3+k,j2+l2);  (73) 

>  Ci(jl  J2  j3M(F(j !+j2  j3)*F(j  1  +j 3  j 2 ) * F (j 2  +j 3 , j  l)/(j  l+j2+j3+l)!)]0'5; 

>  FSi(ji,j2j3,li,l2,l3,k)=  if  FSOi  J2,j3,li,l2,l3,k)=0  then  0 
else(-lf*F(ji+j2+li+l2+l,k)/FS(ji,j2,j3,li,l2,l3,k). 


2.  The  6J  coefficients  are  calculated  by: 

>6j(j'i,j2,j3,ii,i2,i3)=(-iy1+j2+11+12* 

Cl0lj2,j3)*Cl(l1,l2,j3)*Cl(ll,j2,l3)*Ci0‘i,l2,l3)*(FSi(ji,j2,j3,ll,l2,l3,O)+ 
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FSi(jl,j2,j3,ll,l2,l3,4)+FSi(jl,j2,j3,ll,l2,l3,5)+FSi(jl,j2,j3,ll,l2,l3,6)+ 
FSi(ji,j2,j3,ll,l2,l3,’7)+FSi(ji,j2,j3,ll,l2,l3,8)+FSi(ji,j2,j3,ll,l2,l3,9)+ 
FSi(jl,j2,j3,ll,l2,l3,10)+FSi(ji,j2,j3,ll,l2,l3,ll)+FSi(ji,j2,j3,ll,l2, 13,12)+ 

FSi(ji,j2,j3,ll,l2,l3,13)+FSi(ji,j2,j3,ll,  12,13,14)). 


3.  The  six  momenta  vector  summation  coefficient  is  given  by: 

6JJ0l,j2j3,ll,  l2,l3)=  6J(j  l,j2,j  3,1 1,12,13)* 
(-1yi+J2+n+i2*[(2j3+1)*(2l3+1)]0'5. 

4.  The  Clebsch-Gordon  coefficients  are: 


GiG  i,j2,j3,mi,m2,m3)=[(F(j  i+j2,j3)*F0  1  +j  3  ,j  2)  *  F  (j  2+j  3  ,j  1)* 

O'  1  ■ +m  1 ) !  (j  1  ■ -mi) !  0  2+m2) !  O'2-nfi) !  03+m3) !  0  3-m3) !  /(j  1  ■ +j2+j3+ 1 ) !  ] ' °'5; 

FG  (j  1  j2  j3 ,m  1  ,m2,m3 ,z)=z ! *F(j  l+j2-j3,z)*F(j  1  -m  1  ,z)* 

F  0  2+m2  ,z)  *F0  3 +z,j  2-m  1  )*F0  3 +z  ,j  1  +m2) ; 

FG  1 0 1  ,j 2 ,j 3  ,m  1  ,m2 ,m3  ,z)=if  FGO 1  ,j 2 ,j 3  ,m  1  ,m2 ,m3  ,z)=0  then  0 
else  (- 1  )z/FG0 1  ,j 2 ,j 3 ,mi  ,m2,m3 ,z). 

5.  The  3J  coefficients  are  given  by: 

3  JO  1  ,J2  ,j  3  ,iti  1  ,m2  ,m3)=(  - 1  y 1  "j2"m3  *  G 1 0 1  ,j  2  ,j  3  ,m  1  ,m2  ,m3)  * 

(FG  1 0 1  ,j2,j  3  ,mi  ,m2,m3 ,0)+F  G 1  0ij2j3,mi,m2,m3,l)+FGi0ij2j3,mi,m2,m3,2) 
+FGi0i,j2,j3,mi,m2,m3,3)+FGi0i,j2,j3,mi,m2,m3,4)+FGi(ji,j2,j3,m1,m2,m3,5) 
+FG 1  (j  1  ,j  2  ,j  3  ,m  1  ,m2  ,m3 ,6)+FG  1  (j  1  ,j  2  ,j  3  ,m  1  ,m2,m3 , 7)+F  G 1 0 1  ,j  2  ,j  3  ,m  1  ,m2  ,m3 , 8) 
+FGiOi,j2,j3,mi,m2,m3,9)+FGiOi,j2,j3,mi,m2,m3,10)+FGiO'i,j2,j3,mi,m2,m3,ll) 
+FGi0i,j2,j3,mi,m2,m3,12)+FGi0i,j2,j3,mi,m2,m3,13)+FGi0i,j2,j3,mi,m2,m3,14)). 


6.  The  zero  projection  3J  coefficients  are: 

3  JoO  1 J  2  ,j  3  )=(- 1  )9 1 +F+J  3)/2  * 

[FO 2+j  3 ,j  1 ) *F0 1  +j  3 ,j 2) * F 0 1  +j 2 ,j  3)/0 1  +j 2+j 3 + 1 ) !  ] 0  5 * (0 1  +j  2+j 3>/2) !  / 

F(0  l+j2+j3)/2  j  l)/F(0 1  +j  2+j  3>/2  ,j  2)/F  (0 1  "Fj  2+j  3)/2  ,j  3)  • 

7.  The  matrix  elements  of  electrostatic  interaction  are: 

fd(l1,l2,S,L)=fk(l1,l2,L,0)*Cfk(li,l2,0)*F0+fk(li,l2,L,2)*Cfk(l1,l2,2) 

*F2+fk(l1,l2,L,4)*Cfk(li,l2,4)*F4+(-l)s*(gk(li,l2,L,l)*Cgk(l1,l2,l)* 

G1+gk(li,l2,L,3)*Cgk(l1,l2,3)*G3+gk(l1,l2,L,5)*Cgk(li,l2,5)*G5). 

8.  More  specifically,  the  Clebsch-Gordon  coefficients  are  calculated  as: 

KGz0i,j2,j3,nii,m2,m3,z)=if  ji+j2-j3<0  then  0  elif  mi+m2-m3<>0  then 
0  elif  j  1  +j  3-j2<0  then  0  elif  j  2+j  3-j  1  <0  then  0  elif  j  1  +m  1  <0 
then  0  elif  ji-mi<0  then  0  elif  j2+m2<0  then  0  elif  j2-m2<0  then 
0  elif  j3+m3<0  then  0  elif  j3-m3<0  then  0  elif  ji+j2-j3-z<0  then 
0  elif  ji-mi-z<0  then  0  elif  j2+m2-z<0  then  0  elif  j3+z-j2+mi<0 
then  0  elif  j3+z-ji-m2<0  then  0  else  [(2 *j 3+ 1  )* 0 1  +j 2-j  3) !  * 0 1  +j 3-j 2) ! 
*02+j3-j  1)!  *0 1  +m  i ) !  *  0  i-m,)!5Tj2+m2)!5Tj2-m2)!5Tj3+m3)!* 

0  3-m3) !  /(j  i+j2+j3+ 1 ) !  ]°'5*((- 1  )z/(z !  *0 1  +j2-j  3-z) !  *0 1  -mi  -z) !  * 

0  2+m2-z) !  *  0  3+z-j  2+m,) !  *  0  3+z-j  1  -m2) ! )); 
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KGCji,j2,j3,mi,m2,m3)=KGz(ji,j2,j3,mi,m2,m3,0)+KGz(ji,j2,j3,mi,m2,m3,l) 
+KGz(ji,j2,j3,mi,m2,m3,2)+KGz(ji,j2,j3,mi,m2,m3,3)+KGz(j1,j2,j3,mi,m2,m3,4) 
+KGz(ji,j2,j3,mi,m2,m3,5)+KGz(ji,j2,j3,mi,m2,m3,6)+KGz(ji,j2,j3,mi,m2,m3,7)  (81) 

+KGz(ji,j2,j3,mi,m2,m3,8)+KGz(ji,j2,j3,mi,m2,m3,9)+KGz(ji,j2,j3,mi,m2,m3,10) 

+KGz(ji  j2  j3,mi,m2,m3, 1  l)+KGz(ji,j2,j3,mi,m2,m3, 12)+KGz(j  i,j2,j3,mi,m2,m3,13) 
+KGz(j'i,j2,j3,mi,m2,m3, 1 4)+KGz(j ,  j2  j3,m  i  ,m2,m3, 1 5)+KGz(j  i  j2  j3,mi,m2^n3,16) 
+KGzO‘ij2j3, mi, m2,m3, 17)+KGz(ji  j2  j3,mi,m2,m3, 18)+KGz(j  i  j2  j3,mi,m2,m3,19) 

+KGz(ji  j2  j3,mi,m2,m3,20). 

9.  The  Clebsch-Gordon  coefficients  from  the  Wigner  theorem  of  one-electron  operator  matrix 
elements  are  defined  as: 

J3i0'i,j2,j3,mi,m2,m3);KG0'3,j2,ji,  -m3,m2,  -mi)/(2ji+l)0'5.  (82) 

10.  Genealogical  coefficients  are: 

Gen:=matrix([[-  (8/15)0’5,  -  (7/15)0>5, 0,0,0], 

[ -l/5°’5,2/50’5, 0,0,0], 

[(7/30)0’5,  -2/ 1 5°’5,0, l/2°’5,0], 

[  -l/2*(3/5)0,5,  -l/2*(7/5)0,5,2/150’5,  -1/2/30’5,  -l/2*(3/5)0’5],  (83) 

[  - 1  /2  *  (7/5 )0’5 , 1  /2  *  (3/5 )0’5 , 0 , 3/2/70’5 ,  -l/2*(5/7)0’5], 

[(2/5)0’5, 1/1 0°’5,0,  -1/70’5,  -  (5/14)0’5], 

[0, l/2°’5,0,  -  (5/2 1)0,5,(  1 1/42)0'5], 

[0,  -l/2°’5,0,0,l/2°'5]]). 

1 1 .  The  matrix  elements  of  electrostatic  interaction  in  the  d  configuration  are: 

Vqq(n1,n2,k)=3*(Gen[n1,l]*Gen[n2,l]*(  -1)1*63(2,2,1 ,2,2,k) 
+Gen[ni,2]*Gen[n2,2]*(-l)3*6J(2,2,3,2,2,k)+ 

Gen[m,3]*Gen[n2,3]*(  -l)°*6J(2,2,0,2,2,k)+  (84) 

Gen[ni,4]*Gen[n2,4]*(  -l)2*6J(2,2,2,2,2,k)+Gen[n,,5]* 

Gen[n2,5]*(  -l)4*  6J(2,2,4,2,2,k))*25*(3Jo(2,k,2,0,0,0))2. 

12.  The  matrix  elements  of  the  electrostatic  interaction  with  parameters  are: 

q(ni,n2)=Vqq(ni,n2,0)*Fo+Vqq(ni,n2,2)*F2+Vqq(ni,n2,4)*F4.  (85) 

2 

13.  The  matrix  elements  of  the  spin-orbit  interaction  in  the  d  configuration  are: 

W11(ni,n2,Li,L2,Si,S2);3*(-l)(L1+S1_1/2+1)* 

[3/2*(2Li+1)(2L2+1)(2Si+1)(2S2+1)]0  5*(Gen[ni,l]*Gen[n2,l]* 

6J(  1 ,2, Lj,  1  ,L2,2)*6J(  1 , 1/2,S  i ,  1  ,S2,  l/2)+Gen[m,2]*Gen[n2,2]  * 
6J(3,2,L1,l,L2,2)*6J(l,l/2,Si,l,S2,l/2)+Gen[n1,3]*Gen[n2,3]*  (86) 

6J(0,2,Li,  1  ,L2,2)*6J(0,  l/2,Si ,  1  ,S2,  l/2)+Gen[m,4]*Gen[n2,4]  * 

6J(2,2,Li,  1  ,L2,2)*6J(0, 1/2,S  i ,  1  ,S2,  l/2)+Gen[m,5]*Gen[n2,5]* 

6J(4,2,Lh  1  ,L2,2)*6J(0, 1/2,S ! ,  1  ,S2, 1/2)); 
J(n1,n2,L1,L2,S1,S2,J)=(2*3Hs5)0'5*(-l)J+LI+S2  *6J(L1,L2,1,S2,S1,J)* 
W11(n1,n2,L1,L2,Si,S2). 

14.  The  matrix  elements  of  the  tensor  operators  Uk  are: 

Uk(ni,n2,S,L1,L2,k);3*(-l)L1*[(2Li+l)5|s(2L2+l)]0'5*(Gen[n1,l]* 
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(87) 


Gen[n2,l]*(  -l)1*  6J(2,Li,l,L2,2,k)+Gen[ni,2]*Gen[n2,2]*(  -l)3* 
6J(2,Li,3,L2,2,k)+Gen[n1,3]*Gen[n2,3]*(-l)°*  6J(2,L1,0,L2,2,k)+ 
Gen[nb4]*Gen[n2,4]*(  -l)2*  6J(2,Li,2,L2,2,k)+Gen[m,5]* 
Gen[n2,5]  *(  -1  )4*6J  (2,L1,4,L2,2,k)). 

15.  The  matrix  elements  of  the  crystal  field  potential  are: 

Crd3(n1,n2,S,L1,L2,JI,J2,k)=(-l)(S+L2+J1+k)* 

[(2J1+l)*(2J2+l)]°'5*6J(Ji,J2,k,L2,Li,S)*Uk(n1,n2,S,Li,L2,k). 

16.  The  matrix  elements  of  the  angular  momentum  L  are: 

MuL(S,Li,L2,J i,J2);(  _i)(L1+s+J2)*  [(L1*(Li+1)*(2Li+1)*(2J1+1)*(2J2+1)]°  5 

Hs6J(L1,J1,S,J2,L2,1). 

17.  The  matrix  elements  of  spin  angular  momentum  S  are: 

MuS(S,Li,L2,J1,J2)=if  Li=L2  then  2(-l)(L1+s+J1)* 
[S*(S+l)*(2S+l)*(2Ji+l)*(2J2+l)]°  5*6J(S,Ji,Li,J2,S,1)  else  0. 

18.  The  matrix  elements  of  the  electro-dipole  operator  are: 

>  Qi(11,12,S,L1,L2,J1,J2)=(  -1)(11+12+s+j2)* 
[(2L1+1)*(2L2+1)*(2J1+ 1)*(2J2+1)]°  5*6J(12,L1,11,L2,12,1)*6J(L1,J1,S,J2,L2,1) 
*  [(21i+l)*(212-l- 1)]°  5*(-  l)((11+12+I)/2)*3  J(l2, 1  ,li,0,0,0). 

19.  Finally,  the  matrix  elements  of  the  g-factor  are: 

g(S ,L, J)=  1  +( J( J+ 1 )  -L(L+ 1  )+S(S+ 1  ))/2/J/(J+ 1 ). 


(88) 

(89) 

(90) 

(91) 

(92) 


A.3  Program  for  calculation  of  the  eigenvectors  and  eigenvalues  of  the  total  matrix  A 

The  total  matrix  A  is  a  summation  of  the  matrices  of  the: 

1 .  Electrostatic  interaction  with  parameters  F22  and  F44,  which  correspond  to  the  Slater 
parameters  as  Fk=Fkk,  F22=F2/49,  F44=F4/441; 

2.  Spin-orbit  interaction  with  the  parameter  ^34; 

3.  Crystal  field  with  the  second-order  parameters  A20,A2i,A2.i,  A22,A2_2,  and  the  fourth- 
order  parameters  A40,A4i,A4_i,A42,A4_2,A43,A4_3,A44,A4_4. 

Maple  calculates  the  eigensystem  as: 

>  Eigenvals(A,vecs):vecs:  M:=(vecs(%)):Ml:=evalm(l/M):  Ksi:=[eigenvals(A)]:  (93) 
And  the  energy  of  the  transition  between  the  first  and  nth  levels  is  given  by: 

E(n)=Ksi[n]  — Ksi[  1  ] .  (94) 


A.4  Program  for  the  calculation  of  the  matrix  elements  of  the  second  and  fourth  order 
tensor  operators  U2  and  U4 

The  matrices  of  the  tensor  operators  U2  and  U4  of  the  second  and  fourth  order,  respectively,  are 
transformed  with  the  matrix  M  of  eigenvectors  of  the  total  matrix  A: 
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U22=M'1&*U2&*M;  (95) 

U44=M'1&*U4&*M; 

Operator  values  of  the  n  to  m  level  transition  with  energy  Ksi[m]-Ksi[l]  are  given  by: 

>  FiU2(m,n)=  print(m,n,Ksi[n]-Ksi[m],U22[m,n],U44[m,n]).  (96) 

The  effective  cross-section  of  the  n  to  m  state  transition  with  the  use  of  the  experimental 
parameter  a2  is  given  by: 

>  Q(a2,n,m)=a2(U22[n,m])2+a2(U44[n,m])2; 

Qt0=sum(sum((U22[i  ,m])2+(U44|j  ,m])2,j=l .  .4),m=5  ..16);  (97) 

Q  i  (a2,m)=sum(Q(a2,j  ,m),j= 1 .  .4). 

The  calculation  of  experimental  parameter  a2  from  the  lifetime  to  is  the  following: 

a2=4 . 5  2  *  s/ t0/E(6)3/ (2  *  Qt0(ji)+  Qto(c)).  (98) 

Oscillator  strengths  of  the  transitions  of  the  absorption  spectrum  from  levels  1-4  to  levels  5-120 
were  calculated.  The  absorption  spectrum  fabs  vs  hv  is  given  by: 

>  fabs(a2,s2,hv)=sum(Qi(a2,p)*(Ksi[p]-Ksi[l])exp(-(hv-Ksi[p]+Ksi[l])2/(s2)2), 

p  =  5..120)/104*hv.  (99) 

The  oscillator  strength  of  the  transition  from  the  ground  state  to  the  level  m  is  given  by: 

f(a2,m)=Qi(a2,m)*E(m).  (100) 

The  absorption  spectrum  oscillator  strength  fnm  vs  hv  for  level  n  to  m  transition  is: 

fnm(a2,S3,n,m,hv)=Q(a2,n,m)*(Ksi[n]-Ksi[m])*exp(-(hv-Ksi[n]+  (101) 

Ksi[m])2/(s3)2)/104*hv. 

A  table  of  the  initial  STS[1]  and  the  mth  level  STS[m],  the  transition  energies  Ksi[m]-Ksi[l]  and 
the  oscillator  strengths  107f(a2,m)  are: 

>  for  m  from  5  to  120  do  print(m,STS[l],STS[m],(Ksi[m]-Ksi[l]),  (102) 

107f(a2,m))  end  do. 

The  absorption  spectrum  diagram  is  given  by: 

plot(fabs(a2,s2,hv),hv=22000..2000,colour=black).  (103) 

The  oscillation  force  of  transitions  from  an  excited  state  ji  to  the  levels  j2=l..ji-l  is: 

>  fji(a2,S3,ji,hv)=sum(fnm(a2,s3jij2,hv)j2=l..ji-l).  (104) 

The  program  was  developed  to:  1.  calculate  the  oscillator  strengths  and  lifetimes  of  any 


transition,  2.  to  model  the  absorption  and  luminescence  spectra,  and  3.  to  model  the  absorption 
transitions  starting  from  an  excited  state. 

To  use  the  program,  one  needs  to  know  the  lifetime  of  a  level  and  the  experimental  absorption 
spectrum.  This  data  is  used  to  choose  the  parameters  of  the  transition  energy,  intensity,  and 
average  transition  bandwidths  s2,  S3,  to  enable  maximal  experimental  -  modeled  spectra 
agreement.  Finally,  we  used  the  program  to  calculate  the  spectral  characteristics  of  Co“  ,  Fe“  , 
and  Ni2+  in  ZnSe. 
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LIST  OF  SYMBOLS,  ABBREVIATIONS,  AND  ACRONYMS 


IR 

Infrared 

TM 

Transition  metal 

RT 

Room  temperature 

UAB 

University  of  Alabama  Birmingham 

GPI 

Prokhorov  General  Physics  Institute 

PL 

Photoluminescence 

AFM 

Atomic  Force  Microscopy 

CF 

Crystal  Field 

TS 

Tanabe-Sugano 

SO 

Spin-orbit 

MO 

Molecular  orbital 

PMMO 

Phenomenological  method  of  molecular  orbital 

ED 

Electric  dipole 

CFP 

Crystal  field  potential 

cc 

Configuration  curves 

CHV 

Coefficients  of  the  harmonic  vibrations 

cw 

Continuous  wave 

STS 

States 
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